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PREFACE 


The 1988 Johnson Space Center (JSC) National Aeronautics and Space Administration (NASA) / 
American Society for Engineering Education (ASEE) Summer Faculty Fellowship Program was 
conducted by the University of Houston and JSC. The 10-week program was operated under the auspices 
of the ASEE. The program at JSC, as well as the programs at other NASA Centers, was funded by the 
Office of University Affairs, NASA Headquarters, Washington, D.C. The objectives of the program, 
which began in 1965 at JSC and in 1964 nationally, are 

1. To further the professional knowledge of qualified engineering and science faculty members 

2. To stimulate an exchange of ideas between participants and NASA 

3. To enrich and refresh the research and teaching activities of participant’s institutions 

4. To contribute to the research objectives of the NASA Centers 

Each faculty fellow spent at least 10 weeks at JSC engaged in a research project commensurate with 
his/her interests and background and worked in collaboration with a NASA / JSC colleague. This 
document is a compilation of the final reports on the research projects done by the faculty fellows during 
the summer of 1988. Volume 1 contains reports 1 through 14, and volume 2 contains reports 
15 through 26. 
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ABSTRACT 


This report compares the Soviet Space Food and Nutrition programs 
with those of the U.S. The Soviets established the first Space Food programs 
in 1961, when one of the Soviet Cosmonauts experienced eating in zero 
gravity. 

A Soviet scientist recently developed "trophology" - the study of living 
systems which includes assimilation of nutrients. Trophology is expected to 
permit nutritional advances beyond the "Balanced Diet" concept of satisfying 
ongoing metabolic needs. The concept expands and enriches the old con- 
cepts with the new findings in space nutrition. The Soviet scientists have 
conducted a number of studies regarding the concepts of nutrition assimila- 
tion and increased nutritional needs in long duration missions. 

U.S. Space Food and Nutrition programs have been developed over the 
past twenty-five years. From the early days of Mercury and Gemini to future 
Space Station requirements, the U.S. Space Food and Nutrition programs 
have progressively improved. 

This study indicates that some major differences exist between the two 
space food and nutrition programs regarding dietary habits. The major 
differences are in recommended nutrient intake and dietary patterns 
between the cosmonauts and astronauts. The intake of protein, 
carbohydrates and fats are significantly higher in cosmonaut diets compared 
to astronauts. Certain mineral elements such as phosphorus, sodium and iron 
are also significantly higher in the cosmonauts' diets. Cosmonauts also 
experience intake of certain unconventional food and plant extract to resist 
stress and to increase stamina. 
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COMPARISON OF SOVIET AND U.S. SPACE FOOD 
AND NUTRITION PROGRAMS 


INTRODUCTION 


The successful conquest of space was facilitated considerably by the 
creation of a life support system for cosmonauts and astronauts. One of the 
more important elements of life support is nourishment for the spacecraft 
crews. Collective efforts of technologists, nutritionists, engineers, physicians 
and microbiologists have developed a solution for these problems. Foods not 
only provide astronauts or cosmonauts nutrition for normal vital activities, 
but also provide psychological fulfillment. 

For more than thirty years, Russian food industry technologists, scien- 
tists and specialists have studied and developed various forms, types and 
groups of food products. This resulted in a specialized food production in- 
dustry to develop food for the cosmonauts. Nine food products with very 
limited storage lives were developed and tested for Yuriy Gagarin in 1961. 
One of the first tasks of the first man in space was to attempt eating in zero 
gravity. A food ration was then proposed for flights lasting from one to four 
days. The products were canned meats, the first dinner foods, and fruit juices 
in aluminum tubes, and bite-sized pieces of bread. All were considered to 
have very limited shelf lives. (20, 21) 

The primary purpose of the U.S. Space Food and Nutrition Program is to 
provide nutritionally balanced diets to the astronauts during space flights 
and to gather physiological information and conduct nutritional studies de- 
signed to assess the effects of space flights on nutrient metabolism and its 
effects on crew performance. Over the past twenty-five years, the designs of 
food^systems for spacecraft have been developed from the experience of past 
food systems. The foods used on U.S. space flights have been comprised of a 
wide variety of foods which have been specially processed and/or packaged 
to adapt them to zero gravity consumption. Mercury food was experimental 
and transient. Aluminum tubes were used for the semi-solid foods. The first 
manned flight of the Gemini program, Gemini 3, lasted less than five hours, 
but four experimental meals were aboard to test. The meals consisted of de- 
hydrated foods in cubes, dry mix and freeze-dried products providing up to 
2900 Kcal/crewmember. The Shuttle food system has over 150 items for 
astronauts to choose from when selecting a menu. The foods are classified by 
the method of preservation, namely, rehydratable, thermostabilized, natural 
form, and intermediate moisture. (19) 

OBJECTIVES 

1. To study the historical background and progression of Soviet Space 
Food and Nutrition programs. 

2. To acquire knowledge about cosmonauts" daily dietary intake and 
nutritional supplements during longerduration missions. 

3. To study their physiological information and nutritional research and 
assess the effects of space flight on nutrient metabolism on crew 
performance. 

4. To compare the Soviet Food and Nutrition program to the U.S. 
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REVIEW OF THE LITERATURE 


In the Soviet space programs, nutrition is given an important place in 
the life support systems for spacecraft crews. The "Soyuz-9" spacecraft, 
which was launched on July 1, 1970, with a crew of Andriyan Nikolayes and 
Vitaliy Sevastyanov, completed an 18-day flight. Their mission carried a de- 
vice for heating food products in aluminum tubes. This made it possible to 
prepare hot foods onboard the spacecraft. In 1974, aboard the "Salyut-4" 
space station, a device was used to simultaneously heat products in aluminum 
tubes, and bread, which significantly improved the food intake of the space- 
craft crews. (1,20) 

In the Salyut-6 missions, the food technologists, scientists, and other 
specialists expanded the variety of food products, increased their quality, and 
extended their shelf lives. A six-day ration was developed which provided a 
daily four-meal course (first and second breakfast, lunch and dinner). The 
Salyut-6 mission ration in 1977 had over 70 different products including 
twenty types of meat, fish, poultry, and cheese in lOOgm cans, together with 
pureed first and second courses (10 types) in aluminum tubes. There were 
fourteen types of dehydrated (freeze dried) first and second courses and 
garnishes for the canned meats, which could be rapidly reconstituted by 
addinq hot water, Bakery products, and desserts made from twelve types of 
semi-dehydrated fruits were also available. The drink selection included fruit 
juices, tea, coffee, milk, and fermented milk products. 

In Soviet Space Food and Nutrition programs, the make-up of the diets 
are developed on the basis of the flight duration, the complexity of the pro- 
gram, the anticipated energy use, the storage, heating and dispensing 
equipment, the water reclamation systems, ana other equipment. In long 
ternrrspace flights, the relative monotony concerning diet does become 
annoying to some degree. In order to prevent the monotony of the diets and 
uplift the psychological aspects of life, some specially prepared fruits and 
vegetables were transported on an experimental basis. The shipment of fresh 
apples, onions, and garlic also was transported at the request of cosmonauts 
G. Crechko and Yu Romenenko during their Soyuz mission. The use of fresh 
fruits and vegetables was expanded. Oranges, lemons, melons, honey, cran- 
berries, and even fresh cherries, were sent at the request of V. Ryumin during 
his Soyuz mission. The cargo spacecraft were used to send fresh fruits and 
vegetables to Vladmir Lyakhov and Aleksander Aleksandrov, who stayed 
three months in space during Soyuz T-9 mission in 1983. 

The food service provisions aboard the Salyut-6 mission consisted of the 
following elements: Container for stowing and storing foods, tables to make 
meals, an electric warmer, place settings, a device to recycle, measure, and 
dispense hot and cold water, and containers for disposal of packages and left 
over foods. (15,9) 
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TABLE I - MENU EXAMPLE (COSMONAUTS) (4) 


MENU EXAMPLE (COSMONAUTS) 
List of foods for a one-day menu on Salyut-6: 


First Breakfast : 

~ Chicken with Prunes lOOg 

Bread 45g 

Candy 50g 

Coffee with Milk 1 50g 

Lunch : 

Cottage Cheese with Pureed Black Currants 165g 

Honeycake 45g 

Black Currant Juice with Pulp 50g 

Dinner : 

Sauerkraut Soup 165g 

Roast Beef with Mashed Potatoes 57.5g 

Bread 45g 

Prunes with Nuts 60g 

Candied Fruit 50g 

Coffee with Sugar 24g 

Supper : 

Chicken in Tomato Sauce 165g 

Bread 45g 

Cheese lOOg 

Tea with Sugar 23g 


V. P. Bychkov, et. al, 1986, researched the diets of three Salyut-7 prime 
crews. The Salyut-7 cosmonauts experienced the best nutritional systems of 
the Soviet Space Nutrition programs. The rations of the Salyut-7 crews had a 
caloric value of 3150 Kcalories and were well balanced, containing all the 
nutritional requirements. The food system included cosmonauts rations, con- 
tainers for serving and storing food, a dining table, an electric food heater, 
utensils and dishes, devices for regenerating water and measuring hot and 
cold water to packages of freeze-dried food, containers for disposing of 
waste, and a refrigerator for storing fruit and vegetables. The food had 
equivalent food value of diets of Salyut-6, but consisted primarily (65%) of 
freeze-dried products, rehydrated with hot or cold water. (5, 17, 20) 

In a study by V. P. Bychkov, et. al., 1981, concerning the adequacy of the 
protein supplied in the diet of crews of Salyut-6, it was found that the physi- 
cal status of the cosmonauts and the parameters of nitrogen metabolism 
after the missions were indicative of adequate protein intake. Examination 
of nitrogen metabolism in the crew of the first main expedition aboard 
Salyut-6 (96-day mission) revealed that, in spite of individual fluctuations in 
the excretion of the end products of nitrogen metabolism in the post-flight 
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period, their range was within the pre-flight levels. Urinary excretion of the 
total nitrogen was 11.1 gm per day post-flight and 11.6 gm per day pre-flight. 
On the 3rd and 4th post-flight days, it dropped to 9.0 gm per day, then re- 
verted to the base level on the 5th and 7th days. (8) 

In the second expedition aboard Salyut-6 (140-day mission), the post- 
flight excretion of the end product of nitrogen metabolism was close to pre- 
flight levels. On the 3rd post-flight day, there was a 40% decrease in excre- 
tion of total nitrogen, 21% decrease for urea, and 44% decrease for ammon- 
ia, with retention of normal proportion. There was a tendency toward nor- 
malization of these changes on subsequent post-flight days. (8) 

I. G. Popov, et. al, 1982, conducted a study to determine the effect of a 
48-day flight on the blood amino acids content in the crew of Salyut-5. It was 
found that changes in essential and non-essential amino acids ratios, e.g., 
methionine and cysteine, are of major concern. They observed that the 
changes in amino acids metabolism could be due to several nutritional fac- 
tors. The transition from the pre-flight diet to the inflight diet could lead to a 
reduction in the essential amino acids content. A decreased intake of dairy 
products, eggs, leguminous products, freshly prepared food, and an excess 
intake of foods which had gone through severe heat processing and drying. 
They suggested that this offsets the amount of amino acids in foods and their 
accessibility to digestion. They also suggested that changes in living con- 
ditions may change the synthesis of amino acids. Similar studies were 
conducted by I. G. Popov, et. al., 1983, to measure the amino acids levels in 
blood of cosmonauts during the 185-day Salyut-6 flight. It was found that 
appreciable decreases in the concentration of most amino acids occurred in 
both cosmonauts. The threonine and cysteine levels in blood plasma of both 
cosmonauts were below the bottom of the normal range cited in the Russian 
Medfcal Encyclopedia. Alanine and histidine for tne commander, and 
methinine, isoleucine, and arginine for the flight engineer, were also below 
the normal range. The demonstrated changes in plasma amino acids were 
attributed to similar factors cited by Popov in 1982. Recommended counter- 
measures are also similar, such as organizing nutrition rehabilitation for 
cosmonauts by increasing all amino acid intake, particularly methionine, 
cysteine, valine, histidine, tyrosine, arginine, aspartic, and glutamic acids. (23, 
24) 

K. V. Smirnov, et. al., 1982, studied the state of the digestive system 
following long-term space flights. Those studies included: gastric, pancreatic 
and intestinal enzymes. These studies demonstrated consistent changes in 
the digestive system function. The depth and severity were related to dura- 
tion and the condition of the space flights. It was found that weightlessness 
plays an effective role in these changes. There was a noticeable change of 
increased activity of gastric pepsinogen and pancreatic lipase after the first 
day of the 96-day mission. They were normalized by the 25th day of the 
readaptation period for the 96- and 140-day flights, and at the 43rd day after 
the 175-day space flights. They suggested inflight use of preventive measures 
and conditioning of cosmonauts to develop changes in enzyme secretion in 
the gastrointestinal tract during the readaptation period. The marked 
changes were not noted on the first day after the 185-day space flight, in 
comparison with the 175-day flight, due to the preventive measures of ade- 
quate nutrition used during the readaptation period. (27) 
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V. P. Bychkov, et. al., 1982, studied the effects of hypokinesia on man's 
nutritional status. They found that during and after the bed rest studies 
adequate nutrition and certain unidentified nutrients can serve as efficient 
countermeasures against metabolic changes, such as weight loss, negative 
nitrogen, phosphorus, potassium, and sodium balance. The nutrient re- 
habilitation led to a positive nitrogen balance, increase in the utilization of 
vitamin C, B1, B6, and faster termination of compensatory nitrogen re- 
tention, as compared to subjects in the control group. Physical exercise with 
adequate nutrition was found to be the most effective means of preventing 
changes observed under hypokinetic conditions. (3) 

I .A. Radayeva, et. al., 1982, conducted research regarding the biological 
value and shelf life of cultured dairy products in the diet of cosmonauts. 
Freeze-dried cultured milk products, yogurt with sugar, and yogurt with 
fruits and berries were examined for shelf life. It was determined that the 
realistic shelf life was 12 months at 20° ± 5° C and 18 months at 1-4° C. They 
considered these foods to have a high biologic value with respect to their 
protein. The leucine/isolucine ratio is 2:1 for acidophilus paste, 2:3 for sweet- 
ened yogurt, 2:4 for fruit and berry yogurt. It was reported that freeze-dried 
cultured dairy products may be quite beneficial in the diet of cosmonauts and 
were recommended to be incorporated into space diets. (26) 

N. G. Bogdanov, et. al., 1986, studied vitamin levels in cosmonauts 
during pre-flight training and after completion of short-term space flights. 
They reported a statistically significant decrease in the excretion of a number 
of vitamins during the post-flight period. They suggested that it was due to 
an increased vitamin metabolism which leads to an increased need for them 
during post-flight. (2) 

DISCUSSION 


Soviet specialists believe that both a regular meal schedule and carefully 
selected diet are important for the maintenance of overall conditioning in 
space. Nutrition is thought to have a synergestic: effect with other counter- 
measures on the control of adaptive changes, such as musculoskeletal 
strength and mass losses. It has been observed that cosmonauts can actually 
gain weight in space if exercise is combined with vitamins, calcium supple- 
ments, appealing foods, and appetite stimulators, such as onions, garlic and 
spices. (4, 16) 

USSR SPACE FOOD AND NUTRITION CONCEPTS 
FOOD 


In the 30th minute of the 1961 mission, Yuri Gagarin ate and drank. 
This became the first evidence of the possibility of eating, chewing, and swal- 
lowing liquid and solid food in weightlessness. Daily ration of the cosmo- 
nauts on Vostok and Vostok-2 contained about 2800 Kcal, including 100 gm 
of protein, 118 gm of fat and 308 gm of carbohydrates. Foods which were 
used in the flights were packaged in dispenser tubes and included soups, cot- 
tage cheese, as well as drinks: coffee, cocoa, juices. The rations had limited 
shelf lives without refrigeration (up to 5 to 6 days). Meat products in packets 
had to be prepared directly prior to flight. (20, 21) 
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In the Soyuz missions, foods with long shelf life, such as bread, were 
baked in the form of small "one-bite" rolls to prevent crumbs. Meat products 
such as ham, steak, and veal, were also included in the diet. The daily dietary 
intake for cosmonauts in Soyuz had a 3-day menu cycle with four meals per 
day. The Soyuz mission also included dehydrated boiled meat, and Soyuz-9 
cosmonauts were the first to heat meals at 60° to 70° C and drink from the 
dispenser tubes. (20) 

Cosmonauts P. Popvich and Yu Artyukhin on Salyut-3 were the first to 
test dehydrated products, rehydrated with recovered water. These experi- 
ments were continued in Salyut-4. Dehydrated products now comprise up to 
20% of the food rations. The second crew of the Salyut-4 took additional 
food products of limited shelf life, as well as bread, coffee, and tea for the 
first time in the transport ship. The unmanned "Progress" transport ship is in 
widespread use to deliver fresh food rations to the orbital station. (20) 

The water supply system on board the spacecraft functions in conjunc- 
tion with the food systems. The supply of drinking water in the orbital 
station is produced at a rate of up to two liters per man per 24 hours. (17) 

The most recent and satisfactory cosmonaut nutritional system on 
Soviet manned space flights was the one developed for Salyut-7. There is a 
buffet table on Salyut for eating food, as well as a set of table accessories, a 
food heater, facilities for sanitary cleaning of the table accessories, and bags 
for leftovers and packaging. The calorie content of the daily diet was in- 
creased up to 3200 Kcal to combat the negative consequences of weight- 
lessness and physical training exercises during the missions. (5, 1 6) 

NUTRITION 


-Dr. Oleg Gazenko, 1987, in an address to the U. N. committee on the 
peaceful uses of outer space, discussed the concepts regarding digestive 
physiology. He described "a recently developed branch of medicine - 
trophology" as the study of general principles of fundamental vital processes 
of living systems, such as ingestion, processing, and assimilation of nutrients. 
Several discoveries in this area were outlined: New information regarding 
the immunological defense of the small intestine, importance of dietary fiber 
in digestive functioning and overall health, and microbiology of the digestive 
tract in medical support of long-term space flight. (12) 

Dr. Gazenko suggests that trophology will advance the nutritionists' 
knowledge beyond the accepted concept of "Balanced Diet", which means 
satisfying on-going metabolic needs, to the new concepts of "Adequate 
Nutrition." The new concept does not really replace the previous concepts 
but expands and enriches them with new findings in the field of space 
nutrition. (12) 

In research findings, decreases observed in essential amino acids in cos- 
monauts after the 211-day mission led researchers to the conclusion that the 
pre-flight diet should be supplemented with methionine and aspartic acid 
and inflight and post-flight diets with seven essential amino acids plus 
cysteine, arginine, proline, and aspartic acid. (6, 11) 

Soviet scientists believe that increased exercise regimens require an 
additional intake of calories to maintain proper energy balance. They also 
believe that supplements consisting of vitamins, amino acids and minerals 
promote the retention of fluids and electrolytes. They have shown that pro- 
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per dietary combinations can also help to regulate the digestive and en- 
zymatic changes associated with stresses of space flights. (22, 25) 

The Soviet scientists also believe that the plant extract (Eleu- 
therococcus), exercise and high calcium diet are three major ways by which 
reduction in calcium loss is possible in the spacecraft crews. They also use the 
plant extract to resist stress and to increase stamina among the cosmonauts. 
00 ) 

PACKAGING AND PRESERVATION 


Snack products, bread products and appetizers, as well as sweet pastry 
products and fruits, were packaged for the Soviet space missions in film pack- 
ets made of viscotene, a clear plastic film material similar to polyethylene. 
Some of them were vacuum packed. The ration of cosmonauts in the Soyuz 
ship included new shelf stable foods with a long shelf life, and pureed and 
liquid products in dispenser tubes. Meat products - ham, steak, and veal - 
were prepared in the form of meals preserved in metal cans. The sweet 
products included chocolate candy, prunes with nuts, and honey ginger 
bread, all in film packages. Some of the products, in the form of briquettes, 
were covered in edible film. The dehydrated boiled meat included in the 
menu was vacuum packed in a film. (9, 17, 20, 25) 

U.S. SPACE FOOD AND NUTRITION CONCEPTS 


The food systems planned for the U.S. Space Station are detailed in 
reference (28). 

SUMMARY AND CONCLUSION 


The research study shows that the successful conquest of space, either 
for astronauts or cosmonauts, depends on the collective efforts of technolo- 
gists, engineers, physicians, microbiologists, nutritionists and psychologists. 

According to the Congressional Research Service Report prepared by 
Hon. Ernest F. Hollings, May 1988, the Soviets have continued to make steady 
strides toward their goals of having a permanently occupied space station in 
Earth orbit. The Soviets hold a commanding lead in the operational use of 
crews. They have introduced two new launch vehicles and continue to devel- 
op a space shuttle and space plane. (13, 14) 

A comparison time line of U.S. and Soviet space missions is shown in the 
Appendix. 

Cosmonaut Romanenko, after his 326-day flight in space in 1987, stated 
that Mars is getting nearer and nearer to Earth. The Soviets also conducted a 
variety of studies dealing with new concepts regarding nutritional assimi- 
lation and increased needs in long duration missions. However, when com- 
paring the research data, these reports are equivocal regarding counter- 
measures taken, the increased needs, and types of tests used in the analyses. 
(14) 

The present research study indicates that food for U.S. space flight has 
improved steadily throughout the space programs. From the early days of 
Mercury and Gemini to future Space Station requirements, the U.S. Space 
Food and Nutrition program has progressively improved. The U.S. Space Food 
and Nutrition programs are more advanced in terms of their food packaging, 
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preservation techniques and flexible menu patterns for long duration mis- 
sions. (28) However, some major differences exist between the two space 
nutrition programs regarding the respective dietary habits. These differences 
are shown in Tables 2 and 3. 

TABLE 2 - DIFFERENCES IN RECOMMENDED NUTRIENT INTAKE OF 
COSMONAU rS AND ASTRONAUTS 


NUTRIENTS — 

Kilocalories 

Protein 

Fat 

Carbohydrates 

Phosphorus 

Sodium 

Iron 


COSMONAUTS 
3200 Kcal 


1.5g/kg Bwt. (140g) 

1.4g/kg Bwt. (lOOg) 

4.5g/kg Bwt. (395g) 

1.7g 

4.5g 

50mg 


ASTRONAUTS 

2300-3WKcal 
8g/kg Bwt. 

1 .3g/Kg Bwt. (93g) 
4.8g/kg Bwt. (350g) 


9 

18mg 


TABLE 3 - OTHER DIFFERENCES IN FOOD INTAKE OF COSMONAUTS AND 
ASTRONAUTS 


FOOD 

Plant Extract (Eleu- 
therococcus) 500mg/ 
day or I.OOg every 
other day 

COSMONAUTS 

Used to increase 
stamina to resist 
stress 

ASTRONAUTS 

None 

-Garlic 

As food seasoning 

None 

Vodka 

Small amount 

None 

Brandy 

Small amount 

None 

Fresh Fruits and 
Vegetables 

Supplied by Progress 
ship 

Stored 16 hours before 
launch 

Onion, Dill, Parsley 

Cultivated in on- 
board vegetable 
garden (18) 

None 

Multivitamin 

Supplement--twice/ 

day 

Optional (Shuttle) 

Undevit 

Vitamin--twice/day 

None 

Aerovit 

Vitamin--twice/day 

None 

Essential Amino 
Acids (Methionine) 

Supplements/ 
increased amount 

None 

Glutamic Acid 
Decamerit 

Supplement 

None 


1-10 




The cosmonauts' nutrient intake is probably higher than the 
astronauts'. The mean daily inflight nutrient consumption per person during 
Shuttle STS-1 through STS-61C indicate higher intake of nutrients by the 
astornauts compared to the RDA (Appendix). The difficulty exists in 
estimating the nutrient intake since Shuttle crews are not required to 
maintain a food intake log. (29) 

FUTURE RECOMMENDATIONS FOR U.S. SPACE NUTRITION PROGRAMS 

1 . Determine nutrient requirements for longer duration missions. 

2. Understand the increased essential amino acids requirements during 
pre-flight, inflight, and post-flight diets. 

3. Determine the effects of plant extract as a stimulant. 

4. Determine the beneficial effects of alcohol for longer duration missions 
in space. 

5. Study the regulatory effects of proper dietary combinations on 
digestive and enzymatic changes which are associated with stresses of 
long duration flights. 

6. Study the accepted concept of "Balanced Diet", which means satisfying 
only ongoing metabolic needs, to the new concept of "Adequate 
Nutrition", which does not really replace the previous concepts, but 
expands and enriches them with new findings in the field of space 
nutrition. 
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ABSTRACT 


A preliminary study of SOAR, a general intelligent architecture 
for automated problem solving and learning, is presented. The 
underlying principles of universal subgoaling and chunking were 
applied to a simple, yet representative, problem in artificial 
intelligence. A number of problem space representations were 
examined and compared. It is concluded that learning is an inherent 
and beneficial aspect of problem solving. Additional studies are 
suggested in domains relevant to mission planning, as well as, in 
aspects related to SOAR itself. 
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available, the subgoal terminates and pops from the goal stack. The 
goal stack also serves as the anchor to information in working 
memory (WM). Each working memory element (WME) is connected to 
some goal in the stack and can be accessed only by specifying the 
connection from the goal to the WME via augmentations. An example 
of one such connection is: 

(goal <g> A state <s>) 

(state <s> A binding <b>) 

(binding <b> A cell ell A tile <t1>) 

where state and binding are the goal augmentations that provide the 
connection between the goal <g> and the value of the action attribute. 
A chunking mechanism is provided to summarize the system 
behavior in terms of subgoals, and also enables the system to learn 
aspects of problem solving related to subgoals. The overall 
architecture is presented in Figure 1. 

Some of the important characteristics of SOAR are: (1) 
separation between architecture level and the knowledge level; (2) 
problem-spaces for representing knowledge; (3) universal subgoaling 
to resolve impasses ; and (4) production system representation that 
serves as access paths to information in long term memory. In 
addition, there is no conflict resolution mechanism in SOAR. All 
matched productions are fired "in parallel" and add one type of WMEs 
to the working memory. The process of collecting available 
information is called the elaboration phase. The second type of WME 
is the preference element that is used by the architecture in the 
decision cycle to determine the next goal context. The decision 
procedure, sketched in Figure 2, controls the elaboration and decision 
cycles. More detailed descriptions can be found in elsewhere (1-3). 

PROBLEM DOMAIN 

SOAR's capabilities are illustrated below in an output trace of 
the problem-solving process for the "eight-puzzle" problem. The 
board in the eight-puzzle is represented by nine cells occupied by 
tiles numbered one through nine with one blank cell. The numbers on 
the beginning of each line are the decision cycle numbers; 
elaboration cycles are not shown. The Build:P notation signifies a 
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INTODUCTION 


SOAR is a production system architecture for a system capable 
of exhibiting general intelligence. SOAR has three principal 
characteristics separating it from other architectures. These are: (1) 
SOAR can be used to solve a range of problems, from routine problems 
to open-ended problems; (2) SOAR applies a wide range of problem- 
solving methods required for these tasks; and (3) SOAR learns about 
aspects of the problem-solving process and is capable of reporting 
about its performance. 

This document summarizes the work performed in implementing 
simple, yet representative, problems in SOAR. One purpose of this 
exercise was to be acquainted with the SOAR architecture and 
implementation. In the course of this work, some general issues were 
raised and found to coincide with current research topics in SOAR. 

This paper is organized as follows. First, a brief description of 
SOAR is presented. Next, some "toy" Al problems are briefly described 
and their implementations in SOAR are presented. Penultimately, a 
comparison of some problems decompositions is examined and the 
affect of problem presentation on learning and performance is 
discussed. Finally, future studies that can be performed are 
recommended. 

SOAR: AN OVERVIEW 

SOAR is an architecture for exhibiting general intelligent 
behavior. SOAR has evolved from a series of production system 
architectures (1,2). SOAR is embedded in about 255 kilobytes of LISP 
code, and extends 100 kilobytes of modified OPS5.code. 

In SOAR, each task is represented in problem-spaces. The 
problem solving process begins from an initial state, working through 
the state to subsequent states by applying operators. Stages in the 
problem solving process are characterized by a goal context, which 
consists of the current goal, problem-space, state and operator. 
When one stage in the problem-solving process does not have enough 
information, it creates a subgoal (hence the name "universal 
subgoaling") to collect the needed information. The new subgoal is 
then added to the goal stack that keeps track of the goals that were 
created for this problem-space. When the needed information is 
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Soar Architecture 



Figure 1 
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new chunk (rule) being created during the trace. The letters G, P, O 
and S represent the current goal, problem-space, operator and state, 
respectively. 

Learn status: on all-goals print trace 

0 G:G1 

1 P: P3 EIGHT-PUZZLE 

2 S: S4 

3 O: 034 MOVE-TILE(C22) 

C22(S4) ~> S35 


I 2 | 8 | 3 | 


I — I — I- 

I 7 | 6 | 5 | 


4 S: S35 

5 -=>G: G185 (UNDECIDED OPERATOR TIE) 

6 P: PI 86 SELECTION 

7 S: S42 

8 O: (047 046 045) 

9 -- >G: G190 (EVALUATE-OBJECT (M0VE-TILE(C21)) OPERATOR 
NO-CHANGE) 

10 P: P3 EIGHT-PUZZLE 

1 1 S" S35 

12 O: 040 M0VE-TILE(C12) 


C21 (S35) -> S48 



C12(S35) --> S53 
I 2 | 8 | 3 | 


l-l-l-l 

I 0 | 1 | 4 | 

l-l-l-l 

I 7 | 6 | 5 | 


13 S: S53 
Duplicate chunk 
Build:P200 
Duplicate chunk 
Build:P202 
Build:P204 
Build:P205 

14 O: 041 MOVE-TILE(C21 ) 


C21(S35) --> S59 


I 2 | 0 | 3 | 

- l-l-l-l 
I 1 I S | 4 | 

l-l-l-l 

I 7 | 6 | 5 | 


15 S: S59 

16 O: 074 MOVE-TILE(CII) 


C11(S59) -> S64 


I 0 | 2 | 3 | 

l-l-l-l 

I 1 I 8 | 4 | 

|— |— j— | 


I 7 I 6 | 5 | 



17 S: S64 

18 O: 076 M0VE-TILE(C12) 

C12(S64) -> S70 

I 1 I 2 | 3 | 

I— I— I— I 

I 0 | 8 | 4 | 

l-l-l-l 

1716151 

19 S: S70 

20 O: 080 MOVE-TILE(C22) 

C22(S70) --> S77 

I 1 I 2 | 3 | 
l-l-l-l 

I S | 0 | 4 | 

l-l-l-l 

I 7 | 6 | 5 | 

21 S: S77 

goal SOLVE-EIGHT-PUZZLE achieved 
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> (p ri nt-stats) 


Soar 4.4 (external release: created April 19, 1987) 

Run statistics on July 26, 1988 
81 productions (1345 / 5703 nodes) 

229.41667 seconds elapsed (43.133335 seconds chunking overhead) 

21 decision cycles (10924.604 ms per cycle) 

46 elaboration cycles (4987.319 ms per cycle) 

(2.1904762 e cycles/d cycle) 

211 production firings (1087.2828 ms per firing) 

(4.5869565 productions in parallel) 

717 RHS actions after initialization (319.96747 ms per action) 

226 mean working memory size (339 maximum, 296 current) 

2890 mean token memory size (10696 maximum, 3893 current) 

22363 total number of tokens added 

18470 total number of tokens removed 

40833 token changes (5.6184134 ms per change) 

(56.7125 changes/action) 

The particular implementation of this problem in SOAR requires 
twelve productions, which can be divided as follows: (1) four rules for 
setting-up the problem space; (2) one rule to create an operator 
instantiation; (3) two rules for applying an instantiated operator; (4) 
one rule for search control; (5) one rule for monitoring states; and (6) 
three rules for evaluating operators. While the rules described are 
specific to this problem, encoding other tasks typically involve a 
similar division of a problem and hence a similar division of rules. 

LEARNING 

The chunk below is an example of a rule that becomes part of 
production memory during the session outlined above. 


i 
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(SP P218 ELABORATE 

(GOAL <G1 > A OPERATOR { <> UNDECIDED <02> }) 

(OPERATOR <02> 

A NAME EVALUATE-OBJECT A ROLE OPERATOR 

A SUPERPROBLEM-SPACE <P1 > A OBJECT <01 > A SUPERSTATE <S1 > 

A DESIRED <D1 > EVALUATION <E1 >) 

(PROBLEM-SPACE <P1> A NAME EIGHT-PUZZLE) 

(OPERATOR <01> 

A NAME MOVE-TILE A BLANK-CELL <C2> A TILE-CELL <C1>) 

(STATE <S1 > A BLANK-CELL <C2> A BINDING <B1 > { <> <B1 > <B2> }) 
(BINDING <B1> A CELL <C2>) 

(BINDING <B2> A CELL <C1 > A TILE <T2>) 

(DESIRED <D1> A BINDING <D2>) 

(BINDING { <> <B1> <D2> } A CELL <C2>) 

(BINDING <D2> A TILE <T2>) 

— > 

(DISPLACED (MAKE EVALUATION <E1> NUMERIC-VALUE 1) 
A (NLAM-MAKE (QUOTE (EVALUATION <E1> NUMERIC-VALUE 1))))) 

In more transparent form the rule states: 

IF any tile is being moved into its final position 

THEN the next state will yield an evaluation of 1, 

i.e., the next state will be evaluated favorably 

This particular rule is one that is created because of the current 
implementation; the search control rule creates a worst preference 
for those operators not placing tiles in their desired locations. Other 
search control productions are responsible for the creation of chunks 
more in line with their method of search control (4). 

The search control rules are of importance in the outcome of the 
problem-solving process. To explore the effect of these rules on the 
parameters to measure performance, sample runs were performed 
using three different search control strategies. The results are 
presented in Table 1. The third column refers to a search control 
procedure that rejects an operator that will move back a tile into its 
previous location. The fourth column refers to a search control 
mechanism that creates a worst preference to a tile that is to be 
moved out of its desired location. The fifth column corresponds to no 
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*1 

(f. 

82 

268 

39 

73 

383 

339 

72 

281 (114) 
62 
94 
572 
453 

82 

48 (0) 
12 
29 
116 
291 
7 

a 

« c 
c 

» a 

o & 

£ a 

a 

74 

53 

21 

43 

181 

301 

77 

54 (9) 
21 
43 
182 
301 

77 

28 (0) 
12 
25 
104 
279 
3 

Reject * 

nnrin 

73 

177 

48 

84 

442 

276 

81 

112 (24) 
30 
55 
303 
323 

81 

51 (0) 
12 
29 
120 
295 
6 

03 

05 

Q 

number of productions 
elapsed time 
decision cycles 
elaboration cycles 
production firings 
working memory elements 

number of productions 
elapsed time 
decision cycles 
elaboration cycles 
production firings 
working memory elements 

number of productions 
elapsed time 
decision cycles 
elaboration cycles 
production firings 
working memory elements 
chunks used 

Type 

nf run 

No 

Learning 

Learning 

Learned 
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Table 1: Eight-Puzzle Results 




prespecified search control strategy. Each row contains data for each 
type of run. The data contain the number of productions in the system, 
the elapsed time to complete the task, the number of decision cycles, 
the number of elaboration cycles, the number of production firings and 
the maximum number of working memory elements. The last row, 
corresponding to runs that have already "learned," indicates the 
number of chunks that were fired during a problem solving session. 
The numbers in parentheses indicate the amount of time required for 
chunking. 

As expected, chunking added productions to each system. 
However, a learned system did not chunk again. Also as expected, a 
"stronger" search control mechanism (column Four) corresponded to a 
better overall performance: fewer cycles, fewer rule firings, more 
efficient memory and more efficient chunking The reason for this is a 
priori control of the preferences in the "worst*preference" strategy. 
Applying specific preference orientations in the search control 
strategy appears to be an effective way to structure a problem 
implementation (1,4). Finally, it is noted that even without a 
specified search control strategy (column Five), SOAR's default 
strategy is sufficient to solve this problem prior to chunking. 

CONCLUSIONS 

A preliminary examination of SOAR has been performed and, 
from the above-mentioned results, the following conclusions are 
drawn: 

• Learning is a beneficial aspect of automated problem. solving in 
that code can be made more efficient and previously 
unrecognized knowledge, in the form of chunks, can be created. 

• Problem decomposition is the key to an efficient (or even 
successful implementation). The search control 
mechanism(s) used in a specific implementation strongly 
influences the problem-solving process and the learning 
process. 


2-10 



• Measuring the difficulty of a problem is a nontrivial task. The 
performance criteria measured by SOAR need to be 
scrutinized before a definitive measure can be ascertained. 

It is clear that additional studies are needed, with more practical 
problems (such as those presented in (4)), to see if the SOAR 
architecture can be useful in NASA-related applications. 

Acknowledgment 
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microEXPLORER™ and the Lisp environment. 


2-1 1 



References 


1. Laird, J.E., "SOAR User's Manual," Version 4, Xerox Palo Alto 
Research Center, Palo Alto, CA 1986. 

2. Laird, J..E., A. Newell and P.S. Rosenbloom, "SOAR: An 
Architecture for General Intelligence," Artificial Intelligence. 33 
(1987)1-64. 

3. Laird, J., P. Rosenbloom and A.Newell, "Universal Subgoaling and 
Chunking, Kluwer Academic Press, 1986. 

4. Reich, Y. and S.J. Fenves, "Floor-system Design in SOAR: A Case 
Study in Learning to Learn," Technical Report EDRC-1 2-26-88, 
Carnegie-Mellon University, Pittsburgh, PA 1988. 


2-12 



N89 - 20061 


DEVELOPMENT OF A SOFTWARE INTERFACE FOR 
OPTICAL DISK ARCHIVAL STORAGE FOR A NEW 
LIFE SCIENCES FLIGHT EXPERIMENTS COMPUTER 


FINAL REPORT 

NASA/ASEE Summer Faculty Fellowship Projjram — 1988 
Johnson Space Center 


Prepared by: 

Academic Rank: 

University and Department: 


NASA/JSC 

Directorate: 

Division: 

Branch: 

JSC Colleague: 
Date Submitted: 
Contract Number: 


Peter N. Bartram, Ph.D., P.E. 

Associate Professor 

Norwich University 
Computer Science Engineering 
Division of Engineering 
Northfield, Vermont 05663 


Space and Life Sciences 
Life Sciences Project Division 
Project Engineering Branch 
Donald Stilwell 
3 August 1988 
NGT 44-005-803 


3-1 


ABSTRACT 


The current Life Sciences Laboratory Equipment (LSLE) microcomputer for 
life sciences experiment data acquisition is now obsolete. Among the 
weaknesses of the current microcomputer are small memory size, relatively 
slow analog data sampling rates, and the lack of a bulk data storage device. 
While life science investigators normally prefer data to be transmitted to 
Earth as it is taken, this is not always possible. No down-link exists for 
experiments performed in the Shuttle middeck region. One important aspect 
of a replacement microcomputer is provision for in-flight storage of 
experimental data. 


The Write Once, Read Many (WORM) optical disk was studied because of its 
high storage density, data integrity, and the availability of a space-qualified 
unit. In keeping with the goals for a replacement microcomputer based 
upon commercially available components and standard interfaces, the 
system studied includes a Small Computer System Interface (SCSI) for 
interfacing the WORM drive. The system itself is designed around the STD 
bus, using readily available boards. Configurations examined were (1) 
master processor board and slave processor board with the SCSI interface, 
(2) master processor with SCSI interface, (3) master processor with SCSI and 
Direct Memory Access (DMA), (4) master processor controlling a separate 
STD bus SCSI board, and (5) master processor controlling a separate STD bus 
SCSI board with DMA. 


Storage times for 512-byte disk sectors ranged from 56 to 22.5 milliseconds 
without DMA and 53 to 8.9 milliseconds with DMA. The faster times resulted 
from using large blocks of data per disk write request. With the faster times, 
the storage rates significantly exceed anticipated data acquisition rates. 
While DMA is especially attractive, configuring workable systems from 
available boards without hardware modification is difficult. 
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INTRODUCTION 


The success of manned space exploration depends in part on an 
understanding of the physiological effects of microgravity. A major 
concern of life sciences research in space is to guarantee the health and 

safety of astronauts for missions of both long and short duration. 
Additionally, astronauts have experienced many symptoms during exposure 
to microgravity. These latter symptoms have been subsumed under the 
rubric of "Space Adaptation Syndrome." It is part of the charter of the Life 
Sciences Flight Experiments Program (LSFEP) to perform experiments in 

space that will contribute to the understanding of the physiological effects 

of space travel on humans and, where possible, to develop strategies for 

reducing the debilitating effects. 


Most flight investigations require the use of a microcomputer for high 
speed collection, storage, and downlinking of in-flight physiological data. 
The Life Sciences Laboratory Equipment (LSLE) Microcomputer has served 
these needs for life sciences investigations in the shuttle Spacelab module. 
This computer has provided the link between Spacelab experiments and the 
Science Monitoring Area (SMA) in the Life Sciences Project Division 
building at the Johnson Space Center (JSC). To perform this function the 
LSLE microcomputer sends formatted real-time data, in a digital serial form, 
to a High Rate Multiplexer (HRM) in the Spacelab module, which in turn 
transmits the data to the ground via a Tracking and Data Relay Satellite 
System (TDRSS) satellite. 


While the old LSLE microcomputer has served its purposes well, it has 
become obsolescent and must be replaced by a newer design. Shortcomings 
of the old LSLE microcomputer include a lack of replacement units and parts, 
an absence of on-board mass storage or real-time data display, an inability to 
be programmed in a high level language, and poor general specifications 
when compared to today’s standards. Because of these deficiencies, the 
development of a new LSLE microcomputer is now a high priority item. The 
requirements for a replacement LSLE microcomputer include all of the data 
acquisition, HRM downlink, and experiment control functions of the 
current machine. Additionally, the new microcomputer must have better 
and more extensive displays, faster data sampling rates, extensive archival 
storage capabilities for use in middeck experiments where no HRM downlink 
is available, greater processing speed, more memory, compatibility with 
popular commercially available systems such as the IBM PC, and the ability 
to be programmed by the principal investigator using common high-level 
languages. The design should be based to the maximum extent possible on 
commercially available boards and should avoid custom hardware as much as 
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possible. It is hoped that the HRM interface is the only custom design 
required. 


A multiple processor system has several advantages over single processor 
systems. In the multiple processor system, one processor (called the 
"master") is reserved to the maximum extent possible for the principle 
investigator’s use, while additional processors (called "slaves") are dedicated 
to each input/output operation. For example, separate slave processors are 
used for analog to digital conversion (A/D), digital to analog conversion 
(D/A), Parallel I/O, Serial I/O, the HRM interface, optical disk control, etc. In 
this way, work is shared among the many processors. Once this partitioning 
is defined, separate engineering efforts can attack each task and 
development work can proceed in parallel. Expansion would have minimal 

impact on the original parts of the system, as additional processors would 
handle the additional data acquisition load. This modularization facilitates 

software and hardware maintenance. It allows on-going replacement of 
obsolete components at the board level, without significantly disturbing 
other parts of the system. A chief advantage of this concept is that the 
modular nature of the system extends into the hardware itself, thereby 
making development and reconfiguration much simpler. 


One of the most important elements of the new LSLE microcomputer is the 
optical disk. The Write Once Read Many (WORM) optical drive would allow 
the LSLE microcomputer to be used in the shuttle middeck where no 

downlink capabilities are available. This will allow the LSLE microcomputer 
to be used on a much larger number of flights (including those not flying 
the Spacelab), thereby greatly expanding opportunities for the Life 

Sciences Flight Experiments Program. In most cases, the investigator will 

prefer real-time data downlink, but as flight opportunities are limited, he 
may choose to fly his experiment in the middeck rather than to risk losing 

the opportunity altogether. In other cases, the investigator may be given a 
late opportunity to add his experiment to an existing fully developed 

payload. In this case, time may not permit the development of ground 

software to handle the HRM data stream. Again, the use of an optical disk 
will reduce data handling costs and development times. It has been 

estimated that the cost of ground software to handle the HRM data stream for 
life sciences experiments is $250,000 per mission, regardless of the number 
of experiments served. The use of the optical disk instead of the HRM would 
remove this cost entirely. In summary, the use of an optical disk will allow 
single experiments to be inserted into missions on relatively short notice 
and will permit full use of the middeck by life science investigators. Other 

bulk storage systems, such as streaming tape and other magnetic devices, do 
not offer the robustness and high data density of the optical disk. 


The objective of my work this summer was to develop prototype software 
needed to store collected data on an optical disk as well as to recover it. Also, 
timing studies were to be performed to determine if data storage rates would 
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be sufficient for the data collection sample frequency required of the new 
system. 


EQUIPMENT USED 


The commonly used STD bus was chosen for the system backplane. This bus 
is well supported, with several hundred manufacturers offering several 
thousand boards. Thus, prospects are excellent for future upgrades in 
keeping with advances in technology. 


The STD bus processor boards selected for the system are based on the Intel 
8088 and 80188 processors. Advantages of these boards include processor 
compatibility with the familiar IBM PC family and the availability of a wide 
range of software and hardware support products. In the system under 

study, the master processor board is the Ziatech Corporation ZT8815, 
containing an 8 MHz 80188. This board has computing power similar to the 
IBM PC/AT. Of particular interest, the 80188 processor chip contains Direct 
Memory Access (DMA), which provides fast memory -- peripheral data 
transfer in parallel with other processor functions. Additional processor 
boards, slave processors, are Ziatech ZT8830 boards. These contain 8 MHz 
8088 . processors, with computing power approximating that of the IBM PC. 

Both master and slave boards contain a standard SBX interface and connector 

for "piggyback" peripheral modules. Hundreds of modules are available 

from dozens of vendors, providing real-time input and output hardware. 
Thus self-contained data acquisition subsystems can be formed using, for 
example, an SBX analog to digital converter attached to a slave processor 
board. 


The optical disk interface is the Small Computer System Interface (SCSI) 
standard. One configuration uses a slave processor with the Zendex 
Corporation model ZBX-280 SCSI SBX module. The ZBX-280 is based on the NCR 
5380 SCSI controller chip. The same SBX module attached to the master 
processor is an alternative configuration, which allows the use of DMA, not 
available with the slave processor boards. Another alternative SCSI 
interface is the Ziatech STD bus peripheral board, the ZT8850, which 
provides a chip set for SCSI as well as other disk control devices. The ZT8850 
includes its own DMA. Compared to the Zendex ZBX-280, the ZT8850 is simpler 
to program, with more SCSI signals generated by hardware rather than 
under program control. The master processor board controls the ZT8850 as a 
bus peripheral, unlike the ZBX-280 which is directly interfaced and does not 
reside on the STD bus. With the ZT8815 processor and the ZT8850 SCSI boards, 
DMA may be implemented using either the 80188 DMA of the ZT8815, or the 
ZT8850 on-board DMA circuits. 
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The optical disk drive is the Optotech model 5984. This drive was selected 
because Mountain Optech (Boulder, CO) makes a space qualified version of 
the drive (model 200SES) under contract to Goddard Space Flight Center. The 
model 5984 is identical in electrical and software characteristics to the model 
200SES, but is less expensive. Both are 200 megabyte (per side) capacity 
optical disk drives, using 5 1/4 inch removable media cartridges. 


A summary of the equipment used and location of principal vendors is given 
in Table 1. 


SOFTWARE DEVELOPED 


Test programs were written to assure that the optical disk could be both 
written and read in the various configurations considered, and to determine 
read and write times. Best case conditions were used in that no other 
competing activity was required of the processor chip or, where used, the 
STD bus. In all cases, a main program in C calls a collection of assembly 
language functions. For support of the Zendex SCSI piggyback (NCR 5380 
chip) interface with either master or slave processor board, the assembly 
language routines are revisions of routines provided by Mountain Optech. 
Each principal SCSI phase is handled by a separate call: reset, selection 

(including arbitration), status, command, message in, data output, and data 
input. (The message out phase has not been implemented.) The data 
transfer is byte by byte, with polling and handshake signal generation by 
software. The calling program can thus check for proper progression of 
phases and handle error conditions. Additional functions were written for 
DMA use for the data input and data output phases with the Zendex SCSI 
interface used with the ZT8815 master processor board. 


A separate set of assembly language functions, matching the names and 
organization of the above routines, has been prepared for use with the 
master ZT8815 board and ZT8850 peripheral SCSI board. These require no 
significant change in the C calling programs written for other 
configurations. While the example code available from Ziatech for the 
ZT8850 support was influential, these routines followed the organization of 
the functions for the Zendex interface, rather than the task oriented Ziatech 
routines. The routines for the ZT8850 are generally simpler than those for 
the NCR 5380 based Zendex interface. The ZT8850 does not support 

arbitration. The handshake signals required for data transfer are hardware 
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TABLE 1. - EQUIPMENT USED. 


HARDWARE 

Ziatech 8862 Card Cage and Power Supply (STD bus) 

Ziatech ZT8830 Intelligent I/O Control Processor (slave) 

Ziatech ZT8815 80188 based CPU card (master) 

Zendex ZBX-H280 SCSI controller multimodule (iSBX) 

IBM PC compatible personal computer for downloading programs to STD bus 
boards 

Digital storage oscilloscope for obtaining timing data 
Mountain Optech Model 5984 Optical disk drive 


SOFTWARE 

Ziatech 8830 Debug software (and ROMs) to load and debug programs 

Ziatech 8815 Debug software 

Microsoft C Language, version 4.00 

Microsoft Macro Assembler, version 4.00 

Microsoft Linker (loader), version 3.51 

Microsoft MS-DOS operating system, version 3.10 (for program development) 


VENDORS 


Ziatech Corporation 
3433 Roberto Court 
San Luis Obispo, CA 
93401 

(805) 541-0488 


Zendex Corporation 
6700 Sierra Lane 
Dublin, CA 

94568 

(415) 828-3000 


Mountain Optech 
2830 Wilderness Place 
Suite F 
Boulder, CO 

80301 

(303) 444-2851 


Microsoft Corporation 
10700 Northrup Way 
Bellevue, WA 

98004 

(206) 882-8089 
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generated. The data transfer is still byte by byte, with handshake signal 
polling to avoid data loss. Some DMA support routines for the ZT8850 have 
been developed. 


Copies of software developed for this project may be obtained from Peter N. 
Bartram, Division of Engineering, Norwich University, Northfield, Vermont 
05663 (telephone 802/485-2263), or Donald. Stilwell, NASA/Johnson Space 
Center, Mail Code: SE3, Houston, Texas 77058 (telephone 713/483-7308). 


RESULTS 


Polled data transfer can be performed in all configurations studied (ZT8830 + 
ZBX280, ZT8815 + ZBX280, and ZT8815 + ZT8850). DMA writing to the optical 
disk may be performed using the ZT8815 with ZBX280 configuration, 
provided the SBX signal TDMA, not implemented by the ZT8815, is grounded 
(requiring the addition of a wire on either the ZT8815 or the ZBX280). The 
ZBX280 uses this signal for one alternative for ending DMA transfer. With it 
left floating, as on the ZT8815, the ZBX280 attempts to halt DMA prematurely. 
Even with this change, with which DMA writing to disk works well, DMA 
reading of the disk returns incorrect values. The source of difficulty has 

not been determined with confidence. The ZT8850 DMA controller functions 
correctly, provided the memory used is not on the ZT8815 processor card 
controlling it. (A separate memory board was used.) Also at the time of this 
writing, code for using the ZT8815 processor DMA with the ZT8850 is under 
preparation. 


For polled data transfers, timing results were similar for both reading and 
writing. For the LSLE replacement microcomputer, disk writing in real-time 
is critical, whereas reading will be performed at a later time, when speed is 
of lesser importance. With the ZBX280 SCSI interface, writing a single 512 
byte sector each call for data output, the average write time was less than 56 
milliseconds per block, regardless of processor. With each write request 
specifying a ten-sector block (5120 bytes), the average write time for ten 
sectors was 250.4 milliseconds (25.04 milliseconds per sector). Some ten 
sector blocks required 225 milliseconds, others 276 milliseconds. These 
times were identical for both the ZT8830 and ZT8815 processor boards. For 
write request block sizes of 25 sectors (12800 bytes), the average time 
required was 562.5 milliseconds per block (22.5 milliseconds per 512 byte 
sector). This is over 22 kilobytes per second. Since each real-time 
measurement results in a two-byte value, data storage in excess of 11000 
samples per second is possible. 
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For the ZT8850 SCSI with no DMA, write times for single sector requests 
averaged 52.7 milliseconds per sector. For ten sector blocks, the average 
time per block was 197 milliseconds (19.7 per 512 byte sector). With 25 sector 
blocks, the time averaged 369 milliseconds (14.8 milliseconds per 512 byte 
sector). With 25 sector (12800 byte) blocks, over 33.8 kilobytes per second 
can be stored, or over 16,900 measurement samples per second. 


With DMA using the ZT8815 with ZBX280, single sector write times averaged 
52.7 milliseconds per sector. For DMA using blocks of 10 sectors per write 
request, the write time per block fluctuated between 124 and 176 

milliseconds, averaging less than 150 milliseconds per block (15 
milliseconds per 512 byte sector). The time required to write blocks of 25 
sectors fluctuated between 210 and 260 milliseconds, averaging 221 
milliseconds per 12800 byte block, or 8.84 milliseconds per 512 byte sector. 
Thus with DMA, storage rates in excess of 56 kilobytes per second (28000 
measurement samples per second) are possible. 


Using the ZT8850 STD bus SCSI board with DMA (as provided on the ZT8850 
board), single sector write times also averaged 52.7 milliseconds per sector. 
Average write times using blocks of ten seconds per write request also was 
less than 150 milliseconds per block, 15 milliseconds per sector, and the same 
fluctuations in time were observed. However, the times for writing blocks of 
25 sectors were longer: fluctuating between 260 and 312 milliseconds, 

averaging about 269 milliseconds per block. With this, about 46.4 kilobytes 
per second average storage rate results. 


CONCLUSIONS AND RECOMMENDATIONS 


Without DMA, data storage rates are marginally adequate for anticipated data 
collection rates. With DMA, storage rates exceed any envisioned 
requirements. DMA is particularly attractive because during data transfer 
to the disk, the processor is freed for other tasks. 


At this point, it appears better to use the master processor for control of the 
SCSI interface to the optical disk. At the time of this writing, no slave 

processor supporting DMA is available with an SBX piggyback connector 
(for an SCSI interface), though these are under development. In order to 
gain the benefits of DMA, the master must be used to control the SCSI 
interface. Even if the slower transfer rates of currently available slave 
processor boards supporting the SBX piggyback modules (such as the ZBX- 
280 SCSI) are accepted, significant master processor usage would still be 
involved. This is because one slave cannot access memory of others. Thus 
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the master would have to transfer data from data acquisition slaves to the 
SCSI slave memory. With the master controlling the SCSI interface, no 
memory to memory copy is required. The master can send data to the SCSI 
directly from data acquisition slave board memory. 


Block size is the most important factor for write speed. This is probably 
because with large numbers of contiguous sectors written with one write 
command, fewer disk seeks are required to properly position the write 
heads. DMA only marginally improved single sector block write times, 
where seek time seems to be limiting. With larger block sizes, not only is 
speed improved with both methods, but DMA gives markedly improved 
performance. With large amounts of data being written with each output 
command, data transfer time becomes limiting, as the number of seeks is 
reduced. The uneven write times can probably be attributed to the 
occasional write head repositioning required for two contiguously addressed 
sectors located on adjacent disk tracks. 


However, DMA is not without difficulty, especially in mixed-vendor 
configurations. The TDMA signal incompatibility has been mentioned 
earlier. It is suspected that other problems of non-standard standards may 
be involved in using DMA in other configurations. The system designer 

must have a good understanding of both hardware and systems 
programming. 


Several recommendations now can be made. (1) Further work is needed to 
test alternative DMA configurations. (2) As the block sizes were here chosen 
arbitrarily, optimization of the number of sectors written per write request 
in light of memory requirements is required. Along with this, the need to 
interleave data from several data acquisition slaves needs to be considered in 
determining block size. (3) It is important that the timing studies be 
repeated in a prototype system with all parts functioning. With other 
activity on the STD bus and more demands on the master processor, it should 
be verified that data storage rates will be maintained at a level exceeding the 
requirements of the data acquisition slave processors. (4) When slave 
processor boards supporting DMA and the SBX piggyback modules become 
available, the wisdom of placing the SCSI under control of the master 
processor should be reconsidered. 
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ABSTRACT 


NASA programs for manned space flight are in their 27th 
year. Scientists and engineers who worked continuously on 
the development of aerospace technology during that period 
are approaching retirement. The resulting loss to the 
organization will be considerable. Although this problem is 
general to the NASA community, the problem was explored in 
terms of the institutional memory and technical expertise of 
a single individual in the Man-Systems Division. The main 
domain of the expert was spacecraft lighting, which became 
the subject area for analysis in these studies. The report 
starts with an analysis of the cumulative expertise and 
institutional memory of technical employees of organizations 
such as NASA. A set of solutions to this problem are 
examined and found inadequate. Two solutions were 
investigated at length: hypertext and expert systems. 
Illustrative examples were provided of hypertext and expert 
system representation of spacecraft lighting. These 
computer technologies can be used to ameliorate the problem 
of the loss of invaluable personnel. 
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DESCRIPTION OF THE PROBLEM* 


Personnel represent a major share of the 
resources of any organization. As these individuals 
interact with others in and outside of their 
organization, as problems are confronted and solutions 
identified, the organization’s personnel increase in 
value to the organization. The solutions, the 
interactions, the actions of the personnel of the 
organization during their association with it contribute 
to the institutional memory of the organization. Each 
employee is privy to part of this memory to the extent 
the employee has participated in the 
organization's activities. Much of this institutional 
memory resides only within the employees. As such it is 
at risk with the employee. When any employee resigns, 
retires or dies, their unique part of the institutional 
memory is lost. 

This paper examines this problem as it 
confronts NASA. It describes the early stages of an 
effort - using the tools of expert system development 
and hypertext data management schemes - at NASA-Johnson 
Space Center/Life Sciences Directorate to capture some of 
this institutional resource before it is lost. While NASA 
is a unique organization in many respects, its efforts to 
confront and resolve this problem may provide some 
guidance for other organizations which are encountering 
similar problems. 

Loss of Institutional Memory 

Institutional memory is the history of the 
activities and interactions which have occurred within 
an organization. Some of this history is recorded in the 
form of memos, papers, research studies and other 
archival material, but a large share is resident in the 
ephemeral memories of employees of the organization. Much 
of the information is common to more than one individual, 
but each individual has a unique combination of 
elements of the larger whole. They can also claim 
some unique elements of the institutional memory. As 
these individuals are separated from the organization - 
through transfers, retirements or death - the unique 
components of the institutional memory which they possess 
are lost to the organization. 

Not all individuals within the organization can be 
considered to possess equally valuable components of the 
institutional memory. The corporate technical experts 

*The order of authorship is alphabetic: C. Bell, specialized 
in the hypertext solution, R. Lachman in expert systems. 


4-3 


will be more highly valued than the corporate chauffeurs. 
Ideally, a process for the identification of the 
critically important experts and the capture of their 
critical experiences can be in place early in the expert's 
association with the organization. Such a process would 
greatly ease the problems of retaining within the 
organization the expert's uniquely valuable knowledge. 
The failure to identify these individuals prior to their 
leaving the organization will ensure loss of their 

unique components of the institutional memory. 

NASA has a workforce which is rapidly approaching 
retirement age. Due to the unfortunate conjunction of 
similar hiring dates - associated with the start of the 
space program - and the number of individuals rolling 

over their military retirement into the federal 

system, some divisions ( e.g.; Man-Systems Division ) have 
over 50% of their technical and managerial personnel 
eligible for retirement. This is an aging population as 
well, subject to the impact of degenerative disease 

processes and common mortality risks. 

Senior experts serve other roles in organizations 
beside those related to the application of their expertise 
and as a repository of institutional memory. They 
provide mentor services to junior personnel, 
acquainting new hires with successful and unsuccessful 
experimental and behavioral procedures. With the loss 
of the senior individuals the efficiency of the organization 
declines as non-productive techniques are retried, as 

mistakes are repeated by the less experienced new 
personnel . 

Approaches to Solution of Problem 

Several techniques have been tried for managing 
the problems associated with the loss of senior 

personnel . These can be grouped into those designed to 
reduce the rate of personnel loss and those designed to 
compensate for this loss. The organization can attempt to 
reduce the number of personnel lost to retirement by 
providing compensation - increased pay, more recognition, 
etc. - to these individuals. Appeals to corporate 

loyalty or patriotism may be tried to retain vital 
personnel. Legal adjustments - such as removal of 

mandatory retirement laws - have reduced arbitrary 
barriers. If retirement is not preventable, 

efforts can be made to reduce the impact of the 
expert's leaving. These efforts might take the form of long 
term post-retirement consulting assignments. A phased 
retirement package made of gradually extended vacation time 
allotments during the employee's last few years of service 
might help to reduce the impact to the organization of 


4-4 



sudden retirement. 


Several approaches are available! to compensate for 
the loss of personnel. To try to retain some elements 
of expertise at risk, an organization may assign 
apprentices to senior experts. Th€‘.se apprentices can 
reduce the impact of the loss of the expert by 
providing a partial replacement during a transitional 
period. New archival technologies eire being introduced 
which will provide computer assistance to the problem. 
The organization can manage the problem in the traditional 
fashion of accepting the loss of the expert and their 
unique institutional memory components. 

None of the retirement-rate reduction programs 
described compensate for the loss of experts due to 
age related or accidental processes. Consulting 
programs or phased retirement programs are weakened by 
two conflicting issues. To the extent that the expert 
is out of the normal flow of the organization's information 
and duties - referred to as the "loop" - the value of the 
expert's unique contribution is diminished. The more 
attempts are made to keep the expert in the "loop'' the less 
others in the organization are able to accept and begin 
to compensate for the loss of the expert. There is a 
morale problem associated with excessively delayed 
retirement programs as well. The senior expert should 
be able to enjoy a healthy retirement period. Junior 
employees may begin to resent the fact that advancement 
routes are being blocked by the continued presence of the 
senior experts . 

The use of apprentices, trainees or proteges would seem 
to give a semblance of continuity to activities 
requiring high levels of expertise. The success of this 
approach is dependent upon several factors. It 
assumes that there are interested trainees available. 
It assumes that the expert is sufficiently 
introspective to identify the critical features of 
themselves that make them valuable to the organization and 
is able to determine effective techniques to transfer 
these features to another. The expert needs to be able 
to verbalize his expertise and identify situations where 
it may be safely practiced by the trainee. These are not 
common skills nor are they easily learned. Finally it is 
unlikely that the expert can pass along the authority 
and knowledge associated with many years experience in a 
short time period. 

Simple acceptance of continuing losses is the most 
commonly chosen alternative. Organizations such as NASA 
will find the cost of these losses very onerous. This 
leaves NASA with a strong motivation to explore 
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alternate archive generating and organizing procedures or 
other technological fixes. 

An active archival program with a professional staff 
to gather, cross-reference and maintain significant 
documents is an expensive option. A difficulty in using 
such an approach for the problem of the loss of 
institutional memory through personnel attrition is that 
the material which is retained is the material which 
is retainable - printed text - and not the critical 
material of attempts successful or not, of judgements 
validated or disabused. These are retained by the expert 
and make up a crucial portion of the expert's unique 
component of the institutional memory. While these 
experiential items may be vital to the expert and their 
handling of new problems, they are not nearly as likely to 
appear in archival storage as memos or minutes of 
informational meetings. Had the expert maintained an 
active report or diary of problem solutions and attempted 
strategies, of backgrounds to decisions made and 
processes observed, the problems associated with both the 
archival and the apprenticeship approaches would be 
reduced. If the retiring expert can be induced to create a 
set of documents detailing the significant issues used to 
identify and solve problems, the thought processes and 
attack strategies needed, then the problems outlined in 
the preceding paragraphs will recede. It is unfortunate 
that most experts find it extremely difficult to 
accurately verbalize the techniques they use to identify 
solutions. Most are not even conscious of the procedures 
they are applying to the solution. Memory that may be 
critical to the decision making process may not be 
consciously available until needed for the solution of 
a particular problem. The reminiscence approach has a 
further negative associated with it. If a particular domain 
or pattern of problems has not been addressed recently, 
the information may not be brought to conscious memory in 
time to be added to the verbal record. 

Even with complete memoirs of retirees, problems of 
organization, selection, access and structure remain. 
Access procedures must be simple enough to facilitate 
their use by remaining employees. They need to be robust 
enough to survive and be able to produce useful 
information for users who may not have the time to become 
facile in using a sophisticated data search package. 
Cross references between memoirs and other archival 
material need to be maintained, and a data dictionary 
of terms and synonyms established. The standard to 
work toward is an on-line system available to users within 
their desk environment or through simple telephone 
connections - equivalent in ease to contacting the 
former employee. It can be anticipated that a significant 
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portion of the material provided will be redundant, but a 
redundancy of information is to be preferred to a 
shortage of information. A winnowing process can control 
for information duplicated by other retirees or available 
from other sources. 

Experts serve as more than information repositories. 
They achieve their status through the appropriate 
application of their knowledge to the situations they 
confront - they make decisions. While present software and 
theoretical activities point toward efforts to 
duplicate the broad decision making processes of 
experts, the most successful applications thus far have been 
those where the domain of interest is very narrow. Expert 
system tools have been most successfully applied in 
highly bounded situations, where the range of variables 
of interest can be anticipated and their interactions 
predicted and encoded. Unbound situations provide the 
greatest challenge to the extension of this 
technology. The problem of the loss of institutional 
memory and consultative resources is a very unbounded 
situation. 

The replication of the tasks of the expert may be 
approached as a problem in artificial intelligence (i.e. as 
an expert system problem) , or as a data base/archival 
problem. Given the specific needs of the person seeking 
to use the expert as a consultant, a traditional query 
language/key word - synonym archival approach will require 
substantial knowledge of the information domain prior to 
effective use. Many of the individuals contacting an 
expert express great unwillingness to become more 
familiar with the knowledge domain of the expert than 
they already were. Too often keyword searches result in 
too many "hits" to be useful to a simple informational 
search or result in too few "hits" because of a lack of 
overlap between the user's vocabulary and that of the 
expert. Scrolling is a very inefficient way to search a 
data base. 

Given the difficulties outlined above for using 
standard archival approaches, the need for alternative 
approaches seems evident. We will examine two: hypertext 
and expert systems. 

HYPERTEXT 

If sophisticated data storage systems could be 
combined with easy to use linking or searching schemas, 
many of the difficulties of recovering information from a 
text base could be reduced. Such a combination of linked 
knowledge stores and an effective front end interface is 
available in present software shells under the generic 
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term hypertext or hypermedia. 

The concept of hyper-text is quite simple: Windows on 

the screen are associated with objects in a (knowledge) 

base, and links are provided between these objects. 

(Conklin, 1987) . 

General Description 

The term hypertext is used in this paper to describe 
a set of techniques for creating and linking text- 

oriented nodes of information which may be accessed 
dynamically during any given trial. Each node 

represents an item of interest in the hypertext knowledge 
base. Links are established to other nodes containing 
information associated with the node. Links are 
conceptual pathways that are used to move from one node 

to another. A user browses in the hypertext knowledge base 
seeking information. This description of hypertext is 

meant to be general. Each experimental or commercial 
shell or developed package will differ depending upon the 
creator's personal view of the relative importance of 
various features. For an excellent review of the 
specific implementations of hypertext systems, see 
Conklin (1987) . 

Jumps between related nodes occur at the user's 
discretion. The user has the option of pursuing any line 
of inquiry by following links to nodes that may contain 
relevant information. While the links need to be pre- 
established, the exact path chosen through the hypertext 
knowledge base is at the user's discretion and thus not 
rigid. Links may be followed in both directions. 

Any given node may serve as the target for one or more 
links and may have no, one or several links leading from 
it. Frequently, a graph or map of the hypertext 

knowledge base is maintained by the system to serve the 
user as a graphic memory aid. Whether the underlying 

linkage structure should be hierarchical or not seems 
to be at the developers' discretion. Conklin (1987) 
presents arguments for both approaches. Although most 
users inherently favor a hierarchical structure some 
need for the ability to jump between structures seems to 
be supported. Users wishing to jump from one node to 

another may do so with minimal keystrokes. The 
software keeps track of the user's location and 
efficiently permits the user to quickly leave a node to 
move to a higher level (assuming a hierarchical structure) 
or to another node with ease. The simplest analogy to a 
large non- structured hypertext system is to visualize 
interlaced cobwebs with nodes of information at each of 
the cross points. The user is like the spider running 
along the web, each strand representing a link. This is in 
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contrast to traditional data base schemas where the 
pattern of linkage and of searching is rigidly- 
structured (as in hierarchical and network designs) or 
designed for pattern matching (as in relational data base 
organizations) . 

Extraction of nodes 

The creation of the hypertext oriented system 
requires the efficient extraction of nodes from the 
background material. Nodes are edited to represent the 
information they contain in simple useful text. Text 
quantity may correspond to a sentence, a paragraph 
or a screen. Given the desire to produce useful and 
frequently used reference materials nodes should be edited 
with a bias favoring the simple and the short over 
the long and convoluted. As the nodes are extracted, 
they are linked to related nodes. These links are 
established to provide possible routes of access for 
users. While it is impossible to anticipate all 
possible routes of user interest, it is suggested that the 
original source of the material might serve very 
nicely as an initial approximation of the pattern 
of relationships. Cross links are then made to relevant 
topics in nodes extracted from other documents. A 
cognitively manageable structure needs to be 
retained. Commonly, clues to links to other nodes 
have been in the form of embedded references within text 
or menu lists (see Fig. 1 and Fig. 2). Either the 
first line of the node is given or a node title is 
provided to aid the user in deciding whether to use a 
provided link to the node. Both approaches have 
support in commercial systems. 

THE HYPERTEXT PROJECT 

The role of the expert is both as a consultant and as a 
decision- maker. As a consultant the expert is able to 
bring together information from previous work which 
may be relevant to the project at hand. The decision- 
making aspect involves using this information to reach a 
conclusion about the problem. The information 
management aspect of the expert can be 
approximated using hyper-text oriented hypertext 
knowledge base procedures. Information from the 
expert's descriptions of previous projects can be 
merged with data from other sources to approximate the 
knowledge of the expert, and then linked to provide 
a user with access to the information of the expert. 
Suggested sources of material to put into the hypertext 
knowledge base would include retirees' observations, 
technical reports and other collateral material. This 
hypertext knowledge base would be tied to a user 
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Figure 1. - Menu oriented hypertext screen. 


to users the actual methods that Mill lead to a particular decision. 
With hypertext, users directly participate in each and every decision 
that leads to a particular expertise. This process of openly 
displaying the structure CFILE75 STRUCTURE> and uses of knowledge 
readily transmits it to users of hypertext systems. ' * 

SENSITIVITY ANALYSIS Experts machines usually pro^idm a single answer 

supported perhaps by a confidence factor such as 82 percent certainty 
(whatever that means). Hypertext systems allow users to rapidly test 
alternative paths <FILE56 BR0WSE> to see how sensitive the advice may 
be to changes in the initial assumptions. 

In my opinion, given the years of A1 promises and the relative lack of 
success in creation and delivery of workable experts systems, the most 
important pocketbooks are rapidly closing against further use of computer 
processing of rule-based systems as a method of vending expertise from a 
disk. 

Generally, over the last 30 years, technologists haven't made operations 
research practical. Consequently, if you can't build machines that are 
effective in quantitative reasoning, how can you build machines that are 
effective in subjective reasoning. Making subjective decision machines work 
is at least several magnitudes more difficult. <FJLE27 REASONS > 

— More — Page 4 Columns 0-79 Rmargint 77 FILENAME: file21 


Figure 2. - Embedded link hypertext screen. 
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interface designed for easy, effective interactions. 
With the passage of time the organization will progress 
beyond the areas of the technical competency of the 
retired expert. The number of archival inquiries can be 
expected to drop. The nature of those inquiries can be 
expected to change from specific information seeking 
to general viewpoint or attack strategy issues. This 
represents the best presently available solution to the 
problems associated with the need for a partial 
replacement for retiring experts. 

The expert's interpretation of the information 
cannot be duplicated, nor can the expert's judgement as 
to the relative importance to attach to the various 
nodes. The user still must make a judgement as to the 
weight to put on the information so presented. The 
advantage to the user is that the information is available 
quickly. Associated information is rapidly accessed and 
relative importance to the project at hand can be 
judged. The user is integrated with the computer 
hypertext knowledge base, rapidly pursuing leads and 
identifying interesting nodes and rejecting those of 
little interest. While the expert is not available, the 
user is able to make approximation of the expert by 
accessing the same information, linked in the same manner 
as the expert would have linked it. 

NASA has a critical need to establish protocols for 
capturing the expertise of retiring expert personnel. To 
this end an individual was identified who has served as a 
consultant in a variety of areas including lighting and 
window design. A brief sampling of the telephone calls 
received by this individual supported the claim of 
being consulted over fifty times per month on issues 
related to areas of expertise. This particular individual 
is also very verbal and it was hoped that this would also 
help, given the prototype nature of the project. 

Procedure for Gathering Information 

Initial data was gathered during interviews with the 
subject. The interviewers quickly discovered that it was 
necessary to isolate the subject from his normal work 
environment as the demands on the subject for 
consultation were frequent. Interviews were transcribed 
and the transcription files served as the core for the 
hypertext knowledge base. 

Initial interviews were minimally directed. The 
subject was asked to discuss his experiences with 
lighting issues during the development of various 
manned space efforts. Review of the interviews 
identified lacunae which were used as prompting points in 
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subsequent interview sessions. There were two objectives 
being met during the interview series. One was to 
provide a core of material for the development of the 
hypertext knowledge base. The interviews also served as a 
means to isolate some examples of decision making that 
could serve as models during the development of the 
proto-type expert system. Both interviewers were pres- 
ent during the sessions. Some conflicts developed as 
to directions for follow-up questions or deciding the 
point at which a narrative needed to be redirected. Suffi- 
cient information was gathered for both efforts. 

Development of Hypertext Knowledge Base 

Given the goal of creating a useful and used product, the 
hypertext knowledge base uses the narrative material as one 
of several important sources. The second major source is 
a NASA document titled Man-Systems Integration Standard. 
This document is a product of the Man-Systems Division and 
serves as a summary of human factors related information 
for NASA. Other sources useful to this project are 
published reference material dealing with vision and light. 
Material from the latter sources can be updated as new 
interpretations of data cause changes to be made in the 
scientific publications relevant to the topic area. 

The hypertext document being produced for NASA uses the 
commercial hypertext software shell HOUDINI. Efforts 
are * being directed toward examining issues related to 
determining the most effective interface for NASA 
personnel and other technically oriented users. 

Questions to be addressed include: 

How much freedom of movement - cross-linkages as opposed to 
hierarchical links - is needed by the casual user in 
contrast to the experienced user? 

Should the user be able to annotate the hypertext knowledge 
base? 

Should the user be able to add new links or nodes? 

How does training effect user format preferences? 

These will be discussed in future papers by the authors. 

Users of hypertext systems have noted some operational 
problems. There is a significant cognitive overhead 
involved in keeping track of the information sought 
and in identifying nodes of immediately useful information 
giving a large number of irrelevant or partially relevant 
nodes. User disorientation - lost in hyperspace - 
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occurs when the user loses their cognitive map of the 
relationships between nodes. There is the risk of 
a combinatorial explosion of links if nodes are 
thoroughly cross-linked. Each of these areas will be 
actively addressed. Hypertext is not the final solution 
to the problems associated with the loss of 
institutional memory. It does preisent a methodology 
for reducing the effect to the organization of the loss of 
institutional memory through mortality of corporate 
experts. This is an expensive loss. Whether modern 
organizations can acknowledge the threat of this loss 
and make the required investment of time and effort to 
begin to compensate for it, is another question. 

EXPERT SYSTEM APPROACHES 
Project Scope and the Expert System Life Cycle 

In the development of expert systems, the most 
important judgement is the one that determines the 
suitability of a task for expert system support. A positive 
decision in this regard represents the first in a series of 
sequential steps that are recommended for the development of 
an expert system prototype prior to the construction of a 
full-scale knowledge-based system (Waterman, 1986) . 
Prototype development is one phase of the expert system life 
cycle. The life cycle consists of five major phases: system 
conceptualization and formulation of requirements, selection 
of software tools, developing the prototype, constructing 
the end product, and maintaining the product. 

The summer project, according to its initial 
conception, revolved around the professional activities of 
the Man-Systems Division internal specialist in spacecraft 
lighting and windows, CW. Since CW is approaching 
retirement, it appeared prudent to examine the options 
available to NASA in any effort designed to preserve the 
operational experience he has accumulated during a period of 
about 30 years service to the organization. The expert 
system approach was among the two options that were selected 
for detailed examination in consultation with NASA 
colleagues. The scope of the project was conceived as more 
in the nature of a problem analysis and overview than the 
construction of an expert system prototype. Consequently, 
this part of the project is best characterized as a pre- 
prototype problem evaluation and the initial phase of 
prototype construction. 

Typical knowledge-engineering methodology was used. The 
flexibility of the method is an asset but the concomitant 
absence of standards can be a serious impediment to the 
development of a knowledge-based system. Logically, the 
first step in the methodology consists of an evaluation of 


4-13 


the expertise of the domain expert selected for the project. 
It is essential to initially determine if the candidate 
consultant is in fact a genuine expert with respect to the 
subject domain. However, a necessary prerequisite for that 
judgement, as well as subsequent steps in prototype 
development, is an adequate level of familiarity with the 
domain of knowledge. Consequently several weeks were spent 
examining primary sources in lighting engineering and design 
(e.g. Kaufman and Christensen, 1984; 1987). The initial 
acquisition of domain vocabulary and conceptual structure 
was later evaluated, corrected, and extended by questions 
posed to CW during formal knowledge engineering sessions. 
During the period devoted to the study of the literature on 
illumination engineering, a log was kept of the 
consultations that the expert provided to other NASA 
professionals, the meetings that he attended, and the tests 
that he conducted. 

The consultations are best characterized as problem 
solving sessions and covered a surprisingly broad range of 
topics. The expert devoted considerable time to analyzing 
competing designs for the cupola windows in the space 
station. In addition to forecasting general visual effects 
of window design modifications, he constructed AITOFF Equal 
Area Projection of the Sphere showing the kinds of visual 
function (foveal vision, color vision, binocular, etc) that 
could be lost as a result of several contemplated changes in 
window design. He also gave technical advice on attributes 
of space station windows such as size, number of panes, 
shape, replaceability , construction materials and selection 
criteria. The second area consisted of advice on space 
station vibroacoustics, including issues of subliminal 
vibration. Third, he developed TV overlays for viewing the 
grapler on the Canadian arm for handling payload during 
STS-24 . Fourth, he conducted ambient lighting experiments 
that compared the relative efficacy of LCD and gas plasma 
displays for the on board laptop computers that display the 
location of the Orb iter. Fifth, he presented arguments 
dealing with the design of central lighting for workstation 
areas on the space station. (Someone wanted to change the 
overhead lighting specifications from 4800K to 3200K 
ostensibly to achieve superior resolution for TV monitoring 
of workstation activity. CW had photos in his files from 
previous tests that demonstrated the color temperature 
levels that produce better color resolution and rendering 
both to eye and camera.) Sixth, he collaborated on a study 
dealing with mercury containment in fluorescent lighting. 
The seventh and eight were consultations on color rendering 
in still photography and reflectivity of materials in the 
cargo bay of the Orbiter. The ninth problem that he worked 
on during the period of observation dealt with the design of 
an EVA helmet with LCD projection on its viewing surface. 
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The observations of consultative collaboration and 
advisory meetings led us to conclude that the expertise was 
indeed genuine. This conclusion was confirmed in interviews 
with a number of engineers and scientists at JSC who 
occasionally consult with CW. The extent of the expertise, 
however, was unexpected. Expert system technology does not 
work well with broad domains of knowledge and the issue of 
project suitability required that the project be narrowed to 
aspects of the domain of spacecraft lighting design. The 
domain represents one of CW's primary area of expertise. The 
design area may be segmented so that realizable goals may be 
set for pre-prototype problem evaluation and the initial 
steps of prototype construction. 

System Conceptualization 

The specification of system requirements and the 
description of user populations are part of the system 
conceptualization phase of prototype development. Several 
two hour long knowledge engineering sessions were conducted 
with the expert. These were followed by interviews with 
potential users of a lighting design expert system. The 
combination of literature review and interviews lead to the 
formulation of the first cut in the system design: IVA 
(Intravehicular Activity) and EVA (Extravehicular Activity) 
lighting design. Although the lighting principles may be the 
same for the two, the heuristics, rules, and tests required 
appear to be quite different. IVA was selected. The next cut 
was "made in accordance with major areas of spacecraft 
activity: general circulation, habitation, and workstations. 
Emergency and portable lighting was added as a fourth 
category of lighting design. Habitation was subdivided into 
crew quarters, ward room, waste management, personal 
hygiene, and showers. Workstation was divided into 
exercisers, general maintenance, health care, galley, 
windows, and science labs. The development of a prototype is 
typically focused on a restricted subdomain so that the 
potential of a full scale expert system can be demonstrated. 
CRT and other display and control workstation components 
were not included in an overall menu design at this time 
because instrument panel lighting entails unique design 
principles that add considerable additional complexity. The 
problem was deferred for future analysis. Also, an 
enumeration of the full set of components of lighting design 
to be described in a prototype and later incorporated in a 
final system was not attempted. 

The conceptualization phase of development also entails 
the identification of representational and problem-solving 
methodologies that are suitable for the task. These issues 
will be briefly discussed below. 
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The major conclusions of the system conceptualization 
and development analysis are: 

1. The development of a Spacecraft Lighting Advisor is 
a major expert system construction task and accordingly 
requires the methodologies and effort that is typical in 
large-scale knowledge based system development. The design 
experience and advice appropriate for building small 
knowledge systems does not apply and can be seriously 
misleading for the task at hand. 

2. The development of a full prototype system will 
require approximately one man-year provided that the one or 
more members of development staff have the necessary 
training or experience in artificial intelligence 
programming and in the psychology of knowledge engineering. 

3. The finished system may take an additional 1.5 man- 
years or longer depending on the number of features 
included. 

Tool Selection 

The most expensive expert system development packages, 
in general, tend to be the most powerful and possess the 
widest array of representational technologies for 
development. They also include procedures for effective 
execution, as well as the maintenance and modification of 
the -knowledge base. Unfortunately, the most expensive 
packages also require the longest period to master, and full 
mastery only occurs when the tool is used in a real world 
project. There is also the problem of matching the package 
to the problem. For example, Clancey (1985) maintains that 
"when presented with a given 'knowledge engineering tool' 
such as EMYCIN (van Melle, 1979) , we are still hard-pressed 
to say what kinds of problems it can be used to solve." 

Typically, a large-scale problem, such as the one under 
consideration, will require hybrid knowledge representation 
including rules, frames and other object-oriented devices. 
At this stage of the evolution of the technology, however, 
no one can definitively select a knowledge engineering tool 
as appropriate for a given type of problem, let alone the 
one that is optimum for that problem. Yet, effective tool 
selection remains a process that can place boundaries around 
a project and render it computationally and financially 
tractable. 

Early in the project, a decision was made against 
allocating any of the ten week project period to an in-depth 
evaluation of tools, as the required analysis can consume a 
project of short duration leaving little time for anything 
else. Consequently, an expert system shell developed for an 
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advanced graduate course in knowledge-based systems 
(Lachman, 1989) was modified for the task. The first menu of 
the shell contains a selection that jumps the program to a 
module patched into the main program. The patched module was 
used to demonstrate a partial prototype for spacecraft 
lighting. Fig. 3 shows a menu nested at the second level and 
contained in the patched code. The knowledge base developed 
for the pre-prototype system can be readily ported to a 
commercial heuristic programming environment once a tool has 
been selected that was designed for building large-scale 
systems. The final programming environment will require 
demonstrable advantages for maintenance, explanation 
facility, control of reasoning, and knowledge base 
representation of the final product. 


NASA 


Lyndon B. Johnson Space Center 


MAN-SYSTEMS 1 

DIVISION 

ROY LACHMAN 

SP-34 


SPACECRAFT LIGHTING ADVISOR 


SPACECRAFT LIGHTING 
ADVISOR 


Expert System: Select Area For Lighting Design 


GENERAL CIRCULATION 

A PASSAGEWAYS 

HABITATION 

B CREW QUARTERS 
C WARD ROOM 
D WASTE MANAGEMENT 
E PERSONAL HYGIENE 
F SHOWER 


EMERGENCY AND PORTABLE 

G GENERAL 

WORKSTATION 

H EXERCISERS 
I GENERAL MAINTENANCE 
J HEALTH CARE 
K GALLEY 
L WINDOWS 
M SCIENCE LABS 


During Any Procedure: Press <F1> for Main Menu, < H > Exit to DOS 


Figure 3.- A menu from the pre-prototype system. 
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The pre-prototype software includes elements of a 
backward chaining module that, when completed, would 
evaluate the design problem and interactively determine if 
the problem was appropriate for an expert system lighting 
advisor. This is not a trivial part of the project since 
design engineering can involve creative processes that 
cannot be implemented in automata at this time. 

The overall plan for the expert system design is to try 
to duplicate processes generally used by human designers. 
The general theoretical orientation is the Newell and Simon 
(1972) information processing theory of problem solving. The 
problem space, in their theory, consists of a mental 
representation of the initial problem state, mental 
operators that change the problem state, and a termination 
criterion that is an acceptable solution to the problem at 
hand. The problem space consists of all the states or 
situations that, in principle, can be produced from the 
initial state by application of cognitive operators. A 
solution is any sequence of individual state-operator pairs 
that leads to a termination state. A problem-solving 
strategy is a method of search that identifies specific 
paths through the search space as worthy of consideration. 
Engineering design can be generally characterized in the 
search-space formulation. An expert system tries to capture 
from a human domain expert actual search-space strategies 
employed. A prototype should demonstrate the efficacy of the 
heuristic reduction of the search space and selective search 
operators that were copied from the human expert. It also 
should demonstrate the potential value of the final product 
to the sponsoring organization. 

A pre-prototype system, in contrast, displays features 
in a pre-computational fashion that might be included in the 
final system. The goal is to demonstrate the logical 
adequacy of a given set of properties and a given approach. 
In addition, the pre-prototype may include a computational 
implementation of one or two of the features. Major 
properties of engineering design should be captured in the 
software. Engineers, typically, generate and test design 
alternatives for the search space of a given design 
specification. A rough design is produced and backtracking 
through the search space then produces incremental 
improvements in the design. 

The simplest situation of problem solving in 
engineering design occurs when a previous design is 
implemented without modification; specifications are tested 
and a single unidirectional path through the search space is 
produced. Because of the project's time constraints, that 
type of solution was selected for the computational 
demonstration along with an uncomplicated area of lighting 
design on a spacecraft, the crew quarters. A sample 
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87-30-1988 18s 16s 19 


VER-S14 RULE-BASE-SLA2 . RB 


Try To establish one of the following hypotheses froa the rule base: 

1 =HO-DESI6H-CHAM6E:USE-STAIiDARD-15IIATT-LUHIHAIRE:THO-PARALLEL-FIXTURES 

2 =ADAPTI VE-DES I6N : CHAN6E-L0CAT I OM-OF-F I ITURES 

3 =VAR I ANT-DES 1 6N : CHARGE- NUHBER- OF -F I XTURES-AND-LOCAT I ON 

4 =AM-0RI6INAL-DESI8N-REQUIRED: REDESI6N-LUHINA1RE 

ENTER HYPOTHESIS NO. — >1 

The Systea is trying to prove goal N0-DESI6N-CHAN€iE:USE-STAHDARD-15WATT- 
LUHINAIRE: TWO-PARALLEL-FIXTURES 

The Systea is trying to prove goal STANDARD-CREU-QUART. DIMENSIONS 

The Systea is trying to prove goal STANDARD-D1 

The Systea is trying to prove goal DlMENSIONl (HEI6HT)~8BIN. 

All relevant rules have fired, can't prove DIMENSION! (HEIGHT)=8BIN, , which 

is not in STH-DB nor the consequent of a Rule. 

IS THIS TRUE: DIMENSION (HEI6HT)=80IN. <Y!N!D!H) — >Y 

RULE 5 Deduces STANDARD-D1 

Short-tera aeaory (DATA BASE) now contains: 

DIMENSION (HEI6HT)=8IIN. STANDARD-D1 

The Systea is trying to prove goal STANDARD-02 

The Systea is trying to prove goal DIHENSI0N2(NDTH)=41.5IN. 

All relevant rules have fired, can't prove DIHENSI0N2(HIDTH)=41.5IN. 

IS THIS TRUE: DIMENSI0N2(NDTH)=41.5IN. <Y!N!D!H) — >Y 

RULE 8 Deduces STANDARD-D2 

Short-tera aeaory (DATA BASE) now contains: 

DIMENSION! (HEI6HT) s 8BIN. STANDARD-D1 DIHENS!0N2(HIDTH)M1.SIN. STANDARD-02 

The Systea is trying to prove goal STANDARD-03 

The Systea is trying to prove goal DIHENS I0N3 (DEPTH) =39. SIN. 

All relevant rules have fired, can't prove DINENSI0N3(DEPTH)=39.5IN. 

IS THIS TRUE: DIHENSI0N3(DEPTH)=39.5IN. <Y!N!D!H) — >Y 


Figure 4. - A saaple run of the pre-prototype expert systea. 
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IS THIS TRUE: STANDARD-WEIGHT-CONSTRAINTS-APPLY (YIND.'W) — >Y 


The Systea is trying to prove goal STANDARD-HEAT-CONSTRAINTS-APPLY 

All relevant rules have fired, can't prove STAND ARO-HEAT-CONSTRA I NTS-APPLY , 

which is not in STM-DB nor the consequent of a Rule. 

IS THIS TRUE: ST ANDARD-HEAT-CONSTRA INTS-APPLY <Y!N!DiW) — >Y 

The Systea is trying to prove goal ST ANDARD-POWER-CONSTRA I NTS-APPLY 

All relevant rules have fired, can't prove STANDARD-POWER-CONSTRAINTS-APPLY 

IS THIS TRUE: ST ANDARD-POWER-CONSTRA INTS-APPLY <Y!N!D!W) — >Y 

RULE 9 Deduces ST AND ARD -HEIGHT : HEAT: POWER-CONSTRAINTS 
Short-tera aeaory (DATA BASE) now contains: 

DIMENSION (HEI6HT)=88IN. STANDARD-D1 DIHENSI0N2(WIDTH)=41.5IN. STANDARD-D2 

DIHENSIQN3(DEPTH)=39.5IN. STANDARD-03 STANDARD-CREH-QUART. DIMENSIONS 

STANDARD-HEIGHT -CONST RA INTS-APPLY STANDARD-HEAT-CONSTRA I NTS-APPLY STANDARD- 

POUER-CONS TRAINTS-APPLY STANDARD-HEISHT : HEAT : POWER-CONSTRAINTS 

RULE 1 Deduces NO-DES I 6N-CHAN6E: USE-ST ANDARD- 1 5W ATT-LUH INA I RE : THO-PARALL EL - 

FUTURES 

Short-tera aeaory (DATA BASE) now contains: 

DIMENSION (HEI6HT)=8BIN. STANDARD-D1 DIHENSI0N2(UIDTH)-4I.SIN. STANDARD-D2 
DIHENSI 0N3 (DEPTH) =39. 5IN. STANDARD-03 STANDARD-CREH-QUART. DIMENSIONS 

ST ANDARD-HE 1 6HT-C0NSTRA I NTS-APPLY STANDARD-HEAT-CONSTRA I NTS-APPLY STANDARD- 
POHER-CONSTRAINTS-APPLY STANDARD-HEIGHT: HEAT : POWER-CONSTRAINTS N0-DESI6N- 
CHAN6E : USE-STANDARD- 1 5HA TT-LUM INA I RE i THO-PARALLEL-P I X TURES 

FINAL RESULT: NO-DES I 6N-CHANGE: USE-ST ANDARD- 1 5WATT-LUH INAI RE: TWO-PARALLEL- 
FIXTURES 

87-38-1988 16:16:51 NORMAL TERMINATION 

RULE 7 Deduces STANDARD-03 

Short-tera aeaory (DATA BASE) nou contains: 

DIMENSION (HEIGHT) =88IN. STANDARD-Di DIMENSI0N2(WIDTH)=41.5IN. STANDARD-02 

DIMENSI0N3(DEPTH)=39.5IN. STANDARD-D3 

RULE 8 Deduces STANDARD-CREW-QUART. DIMENSIONS 

Short-tera aeaory (DATA BASE) nou contains: 

DIHENSI0N1 (HEIGHT)=88IN. STANDARD-DI DIHENSI0N2(WIDTH)=41.5IN. STANDARD-02 
DIMENSI0N3(DEPTH)=39.5IN. STANDARD-03 STANDARD-CREH-QUART. DIMENSIONS 
The Systea is trying to prove goal ST ANDARD-HE I GHT : HEAT s POUER-CONSTRA I NTS 
The Systea is trying to prove goal ST ANDARD-UE I GHT -CONSTRA INTS-APPLY 
All relevant rules have fired, can’t prove ST ANDARD-HE 1 6HT-C0NSTRA INTS- 
APPLY, uhich is not in STH-DB nor the consequent of a Rule. 


Figure 4. (continued) - A saaple run of the pre-prototype expert systea. 
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demonstration run is shown in Fig 4. 

It should be possible to develop a spacecraft lighting 
advisor in an incremental fashion starting with the least 
problematic design areas and working up to the most 
difficult. Successful development of a heuristic system that 
evaluates the difficulty and feasibility of design problems 
for solution by an expert system is not guaranteed at this 
stage in the evolution of the technology. 

SUMMARY 

In this paper, we have examined the problems associated 
with the loss of specialized expertise and institutional 
memory through the attrition of personnel. The costs to the 
NASA organization were examined and several possible 
solutions were critiqued. Two solutions were selected for 
detailed examination: hypertext and expert systems. A 
paradigmatic case was selected in the person of an in-house 
expert on spacecraft lighting and windows. Sample 
computational implementations were produced and described. 
These were determined to offer reasonable solutions to at 
least part of the problem. 
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ABSTRACT 


This paper is concerned with the application of 
parallel techniques for electrical power system analysis. 
The Newton-Raphson method of load flow analysis was used 
along with the decomposition-coordination technique to 
perform load flow analysis. The decomposition-coordination 
technique enables tasks to be performed in parallel by 
partitioning the electrical power system into independent 
local problems. Each independent local problem represents a 
portion of the total electrical power system on which a load 
flow analysis can be preformed. The load flow analysis is 
preformed on these partitioned elements by using the Newton- 
Raphson load flow method. These independent local problems 
will produce results for voltage and power which can then be 
passed to the coordinator portion of the solution procedure. 
The coordinator problem uses the results of the local 
problems to determine if any correction is needed on the 
local problems. The coordinator problem is also solved by 
an iterative method much like the local problem. The 
iterative method for the coordination problem will also be 
the Newton-Raphson method. Therefore, each iteration at the 
coordination level will result in new values for the local 
problems. The local problems will have to be solved again 
along with the coordinator problem until some convergence 
conditions are met. 
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INTRODUCTION 


This paper will use the decomposition-coordination 
technique which enables task to be performed in parallel 
when solving sets of nonlinear equations (1, 4, 5, 6). Sets 
of nonlinear equations occur when performing a load flow 
analysis on an electrical power system (3, 7, 8) . Load flow 
is the solution of an electrical network that gives the 
values of currents, voltages, and power flows at every bus 
(node) in the electrical power system (3, 7, 8) . In the load 
flow problem nonlinear relationships between voltage and 
power occur at each bus (3, 7, 8). The values of the voltage 
and power must be solved for at each bus so that the 
response of the electrical power system can be determined. 
With the increase in the size of new proposed space based 
electrical power systems, it will become necessary to have 
very fast simulation (solution) of these systems. Very fast 
simulation of the electrical power system will aid in the 
evaluation of the system performance by decreasing the speed 
of calculations (1, 3, 4, 6, 7, 8). 


LOAD FLOW ANALYSIS 


The Newton-Raphson method for preforming the load flow 
calculation was used (2, 3, 7, 8) . Taylor series expansion 
for a function of two or more variables is the basis of the 
Newton-Raphson method. Partial derivatives of order greater 
than 1 are neglected in the series terms of the Taylor 
series expansion. The Newton-Raphson method was used 
because it calculates corrections while taking into account 
all other interactions. The number of iterations required 
by the Newton-Raphson method using bus admittances is 
practically independent of the number of buses (3, 7, 8). 
For these reasons shorter computer time for a solution of 
the load flow problem could occur when analyzing large 
electrical power systems (3, 7, 8) . 

The solution of the load flow problem is initiated by 
assuming voltage values for all buses execpt the slack 
bus (3, 7, 8). The slack bus is the point at which the 
voltage is specified and remains fixed. The voltage at the 
slack bus is fixed because the net power flow of the system 
cannot be fixed in advance until the load flow study is 
complete (3, 7, 8). The power calculation at the slack bus 
supplies the difference between the specified real power 
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into the system at the other buses and the total system 
output plus losses (3, 7, 8). The Newton-Raphson method for 
load flow analysis will be used to solve the load flow 
problem at the local and coordinator problem levels (2, 3, 7, 
8 ) . 


DECOMPOSITION-COORDINATION METHOD 


The decomposition-coordination method enables tasks to 
be preformed in parallel by partitioning the electrical 
power system into independent local problems (1, 4, 6). The 
power system presented in Figure 1 was solved by ' defining 
three local problems and a coordinator problem(l, 4, 6). 



FIGURE 1. -CIRCUIT DIAGRAM OF POWER SYSTEM TO BE ANALYZED. 


The local problems were solved by the 


procedure (1, 2, 3, 4, 

6, 7, 8) : 


v n+1 = v n - 

3F i (P i ,v i ) 

1 

v i v i 

3v^ 

v i- v ij 


Newton-Raphson 


F. 


n 


• (P ' 


V i^ 


. ( 1 ) 
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Starting with an initial guess for the voltage, equation (1) 
was used successively until convergence was achieved (1, 2, 
3, 4, 6, 7, 8) . Convergence is achievied when the power 
equations are in balance at each bus (3, 7, 8). This is 
represented by the following equation: 

F i (Pi* = o 


The system of equations which describe the three local 
problems are given in Table 1. 


TABLE 1 . -LOCAL PROBLEM EQUATIONS 


F (P ,V ) - -P + 14E -10E E -4E E =* 0 

11 11 11 11 11 12 11 13 


F(P ,V ) = -P -10E E +15E -5E E =* 0 

12 12 12 11 12 12 12 13 


F (P ,V ) = -P -4E E -5E E +9E = 0 

13 13 13 11 13 12 13 13 


F(P , V ) = -P +2E -E E -E E = 0 
13 13 13 13 13 23 13 33 


F (P ,V ) = -P -E E +2E -E E = 0 
23 23 23 12 23 23 23 33 


F (P ,V ) = -P -E E -E E +2E = 0 

33 33 33 13 33 23 33 33 


F<P ,V > * -P +14E -10E E -4E E = 0 

21 21 21 21 21 22 21 23 


F(P , V ) = -P -10E E +15E -5E E = 0 

22 22 22 21 22 22 22 23 


F(P ,V ) - -P -4E E -5E E +9E = 0 

23 23 23 21 23 22 23 23 


The solution for the local problems are usually not 
available in an explicit form. Therefore, the coordinator 
problem must be solved iteratively. This implies that for 
each iteration at the coordinator level, new values for the 
input to the local problem will result (1, 4, 6). 
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The Newton-Raphson procedure was also used for the 
solution of the coordinator problem: 






v k“ v k 


F k^k’ V 


Starting with an initial guess for the voltage (these values 
should come from the last solution of the local problem) , 
equation (3) was used successively until convergence was 
achieved(l, 2, 3, 4, 6, 7, 8) . Convergence is achieved just 
like that of the local problem when the power equations 
balance to zero (or within a preset error) . This is 
represented by the following equations: 

Fv (Pv» v v) = 0 . (4> 


The equations which describe the coordinator problem are 
given in Table 2 . 


TABLE 2. -COORDINATOR PROBLEM EQUATIONS. 


p/p v ) = -P + 14E -10E E -4E E =0 

31 31 31 31 31 32 31 33 


p (p ,v ) = -P -10E E +15E -5E E = 0 

32 32 32 31 32 32 32 33 


F(P ,V ) = -P -4E E -5E E +9E = 0 

33 33 33 31 33 32 33 33 


EXAMPLE SOLUTION AND RESULTS 


To solve the load flow problem of Figure 1 the local 
and coordinator problems of Table 1 and 2 were first solved 
by the following process: 


1. There were three local problems designated and 
slack buses were assigned (given Table 3). 
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2. The local problems were then solved to 
determine what value of voltage and power must 
be present at the buses of each local system 
to meet the load requirement of each local 
problem. 

3. The coordinator problem was then solved by 
assigning a slack bus (given Table 3) and using 
the most recent values of the voltage and 
power from the local problem. 


This process was continued until the local and coordinator 
problems converged on a value of voltage for all buses. The 
results of this process are given in Table 3. 

TABLE 3. -SOLUTION OF THE LOAD FLOW PROBLEM FOR FIGURE 1. 

Local Problem Iteration 1. 

Bus Power Voltaae 3us Power Voltage Bus Power Voltage 

11 120.000 30.588 13 -30.579 30.303 21 120.000 30.738 

12 -150.000 30.228 23 20.017 30.595 22 -100.000 30.440 

13 30.579 30.500* 33 11.213 30.500* 23 - 20.017 30.500* 

Coordinator Problem Iteration 1. 

3us Power Voltage 

31 112.051 30.500* 

32 -100.000 30.212 

33 - 11.213 30.299 


Local Problem Iteration 
Bus Power Voltage Bus Power Voltage 

11 120.000 30.388 13 -31.703 29.797 

12 -150.000 30.026 23 18.874 30.359 

13 31.703 30.303* 33 13.399 30.299* 


2 . 

Bus Power 
21 120.000 
22 - 100.000 
23 - 18.874 


Voltage 

30.828 

30.532 

30.595* 


Coordinator Problem Iteration 2. 
Bus Power Voltaae 

31 114.446 30.500* 

32 -100.000 30.209 

33 - 13.399 30.289 


Local Problem Iteration 3. 

Bus Power Voltage Bus Power Voltage Bus Power Voltage 

11 120.000 29.882 13 -31.763 29.785 21 120.000 30.534 

12 -150.000 29.515 23 18.877 30.348 22 -100.000 30.235 

13 31.763 29.797* 33 13.458 30.289* 23 - 18.877 30.299* 

Coordinator Problem Iteration 3. 

Bus Power voltage 

31 114.507 30.500* 

32 -100.000 30.208 

33 - 13.458 30.288 


Bus Power 
11 120.000 

12 -150.000 

13 31.765 


Local Problem Iteration 4. 


Voltage Bus Power Voltage 

29.871 13 -31.765 29.785 

29.503 23 18.880 30.348 

29.785* 33 13.449 30.289* 


Bus Power 
21 120.000 
22 - 100.000 
23 - 18.889 


Voltage 

30.583 

30.284 

30.348* 


Coordinator Problem Iteration 4. 
Bus Power Voltage 
31 114.498 30.500* 

23 -100.000 30.209 
33 - 13.449 30.289 


Slack Bus 
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As can be seen from the results of Table 3 convergence of 
the values of voltage and power have occured to within an 
error of 0.01. To continue this solution procedure for 1 
more iteration will result in convergence at errors of less 
then 0.001 for the voltage and power. 


CONCLUSION 


This research used the decomposition-coordination 
technique along with the Newton-Raphson method to preform 
load flow analysis. It has been shown that parallel 
techniques can be used to solve standard power system 
problems. The solution procedure required the formulation 
of a set of local problems which could be run in parallel. 
The results of these local problems were then passed to the 
coordinator problem to determine if any correction is needed 
to the local problem. This process was continued until 
convergence conditions were met. 
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ABSTRACT 


The sain objective of this research is to develop an 
algorithm for isolated-word recognition. This research is 
focused on digital signal analysis rather than linguistic 
analysis of speech. Features extraction is carried out by 
applying a Linear Predictive Coding (LPC) algorithm with 
order of 10. Continuous-word and speaker independent 
recognition will be considered in future study after 
accomplishing this isolated word research 

To implement and test the proposed algorithm a 
microcomputer-based data acquisition system has been designed 
and constructed. The system digitizes the voice signal , 
after passing through a 100 c/s-3.8 Kc/s band pass filter, 
with a sample rate of 8 KHz and stores the digitized data 
into a 6 4Kx 10 -dynamic random access memory (DRAM) buffer. A 
squelch circuit consists mainly of comparators (741s85) 
detects the beginning of the spoken word. The end of the 
word is detected by a software algorithm based on comparing 
the speech energy with a precalculated threshold. A flag 
signals the end of sampling and the data is transferred from 
the buffer to an IBM-PC where it is segmented into frames 
each, 30 millisecond long, and the LPC coefficients are 
calculated. 

To examine the similarity between the reference and the 
training sets, two approaches are explored. The first is 
implementing traditional pattern recognition techniques where 
a dynamic time warping algorithm is applied to align the two 
sets and calculate the probability of matching by measuring 
the Euclidean distance between the two sets. The second is 
implementing a backpropagation artificial neural net model 
with three layers as the pattern classifier. The adaptation 
rule implemented in this network is the generalized least 
mean square (LMS) rule. 

The first approach has been accomplished. A vocabulary 
of 50 words was selected and tested, the accuracy of the 
algorithm was found to be around 85%. The second approach is 
in progress at the present time. The topology of the 
backpropagation model consists of three layers: input, 
hidden, and output. The actual output of each node is 
calculated using a sigmoid nonlinearity function of the inner 
product of the weight and the input; the weights are adapted 
by using the formula W^j new = W i j° ld + u EjX* where u is the 
gain factor, Ej is the'' error, and X^ is J the input. The 
network is being simulated on a PC. 
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INTRODUCTION 

For more than a decade the United States Government, 
foreign countries especially Japan, private corporations, and 
universities have been engaged in extensive research on 
human-machine interaction by voice. The benefits of this 
interaction is especially noteworthy in situations when the 
individual is engaged in such hands/eyes-busy task, or in low 
light or darkness, or when tactile contact is 
impractical/ impossible. These benefits make voice control a 
very effective tool for space-related tasks. Some of the 
voice control applications that have been studied in NASA-JSC 
are: VCS Flight experiments, payload bay cameras, EVA heads 
up display, mission control center display units, and voice 
command robot. A special benefit of voice control is in zero 
gravity condition where voice is a very suitable tool in 
controlling space vehicle equipment. 

Automatic speech recognition is carried out mostly by 
extracting features from the speech signal and storing them 
in reference templates in the computer. These features carry 
the signature of the speech signal. These reference 
templates contain the features of a phoneme, word, or a 
sentence, depending on the structure of the recognizer. If a 
voice interaction with the computer takes place, the computer 
extracts features from this voice signal and compares it with 
the reference templates. If a match is found, the computer 
executes a programmable task such as moving the camera up or 
down. 


Several digital signal processing algorithms are availa- 
ble for speech feature extraction. The efficiency of the 
current algorithms is limited by: hardware restriction, exec- 
ution time, and easiness of use. Some of these algorithms 
are: Linear Predictive Coding (LPC), Short-time Fourier Anal- 
ysis, and Cepstrum analysis. Among these algorithms, LPC is 
the most widely used since it is easy to use, has short exec- 
ution time, and do not require large memory storage. 
However, this algorithm has several limitations due to the 
assumptions upon which it is based upon. 

Current speech recognition technology is not sufficiently 
advanced to achieve high performance on continuous spoken 
input with large vocabularies and/or arbitrary speakers. A 
major obstacle in achieving such high performance is the 
limited capability of the traditional pattern recognition 
(classifier) algorithms that are currently implemented. Due 
to this limited capability, and considering the fact that 
humans have a fascinating capability of recognizing the 
spoken words, researchers have started to explore the 
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possibility of implementing human-like models, or what is 
known as artificial neural networks, as the pattern 
classifiers. 

Neural net models have the greatest potential in area 
such as speech and image recognition where many hypotheses 
are pursued in parallel, high computation rates are required, 
ability to learn is desired, and the current best systems are 1 
far from equaling human performance. Most neural net 
algorithms adapt connection weights in time to improve 
performance based on current results. Adaptation or learning 
is a major focus of neural net research. The ability to adapt 
and continue learning is essential in area such as speech 
recognition where training data is limited and new talkers, 
new words, new dialects, new phrases, and new environments 
are continuously encountered. 

LINEAR PREDICTIVE CODING (LPC) 


Feature extraction in our research is carried out 
through a Linear Predictive Coding algorithm. The signal 
corresponding to a spoken word is segmented into frames each 
30 milliseconds long and the algorithm is applied to replace 
each frame with 10 coefficients. In the following, we 
briefly review the LPC algorithm. Details of this algorithm 
can be found elsewhere [1-9]. This algorithm is built on the 
fact' that there is a high correlation between adjacent 
samples of the speech signal in the time domain. This fact 
means that an nth sample of speech signal can be predicted 
from previous samples. The correlation can be put in a 
linear relationship as: 


Y n = a l^n-l +a 2 Y n-2 + * * * +a p Y n-p 


( 1 ) 


where p is the order of analysis. Usually p ranges from 8 to 
12. y n is the predicted value of speech at time n and a*s 
are the linear predictive coefficients. The prediction error 
B n that resulted from the above linear relationship is: 



( 2 ) 


To find the predictive coefficients which give least mean 
square error, the above equation is squared, partially 
differentiated with respect to a's, and time average term by 
term. The result is p equations in p unknowns as shown 
below: 
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(4) 


where rj is a correlation coefficient of waveform (y n ) and 
r_j = rj by the assumptions of stationary state of y n The 
coefficients a^’s exist only if the matrix in equation 3 is a 
positive definite. To ensure that this condition is 
satisfied, y n is multiplexed by a time window W . This 
multiplexing makes y n exist in a finite interval from 0 to 
N-l, where N is the interval of the Window; a stable solution 
for equation 3 is always obtained. Accordingly, r j can be 
written, as: 


N-j-1 

r j •* ^ < 5 > 

In our study, a Hamming window is implemented. 
Calculation of the correlation coefficients by window 
multiplexing is called the correlation method. a i’ 8 
correspond to the resonance frequencies of the signal, and if 
p, the order of the analysis is selected correctly, these 
ai's represent the formants, frequencies at which peaks of 
the power spectrum of the speech signal occur. A block 
diagram representing an algorithm for voice recognition based 
on LPC analysis is shown in Figure 1 . 
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FIGURE 1. CALCULATION OF THE LPC COEFFICIENTS [5] 
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EXPERIMENTAL DESIGN 

To calculate the LPC coefficients and test the 
performance of the proposed algorithm, a microprocessor-based 
data acquisition system has been designed and constructed. 
Figure 2 shows a block diagram of the system. The system is 



FIGURE 2. A BLOCK DIAGRAM OF THE DATA ACQUISITION SYSTEM. 

designed to lay the foundation for further expansion and 
enhancement for more sophisticated microprocessor-based 
speech identification/ recognition research. The system 
receives the voice signal through a microphone coupled with 
an audio amplifier. The voice signal passes through an active 
multiple feedback bandpass filter with a 3db bandwidth of 3.5 
KHz approximately. The output signal of the filter is 
applied to the analog-to-digital converter where it is 
digitized with a sampling rate of 8 KHz. A squeloh circuit 
is constructed to detect the beginning of the utterance and 
accordingly activates a temporary storage buffer to store the 
digitized data. The buffer consists of DRAMS with maximum 
capacity of 64 Kbyte. A hardware flag (the output bit of a 
flip flop) signals the end of sampling and the data is 
transferred to the microprocessor where digital signal 
analysis and pattern recognition algorithms are applied. The 
digital-to-analog converter, power amplifier and loud speaker 
are used to verify the storage. If the output of this 
circuit matches the original signal, then the storage is 
successful. The system has been tested successfully by the 
aid of a function generator. Details of the hardware of the 
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system can be found elsewhere [9]. 

PATTERN RECOGNITION 

To recognize the spoken word, an algorithm that compares 
between the reference and the training patterns to see 
whether they match or not should be developed. Two 
approaches will be discussed here. The first is a 
traditional pattern classifier approach where training 
pattern is aligned in time with the reference pattern and the 
Euclidean distance between the two patterns is calculated and 
taken as the probability of matching. The second is based on 
implementing a backpropagation artificial neural network as 
the pattern classifier. In the following, we discuss briefly 
the two approaches . 

A. Dynamic Time Warping (DTW) 

The dynamic time warping (DTW) algorithm finds the 
"optimal" (least cost) warping path w(n) which minimizes the 
accumulated distance, D, between training and reference 
patterns, subject to a set of path and endpoint, constraints. 
The dynamic programming algorithm is based upon the fact that 
the optimal path to point (i,j) in the two dimensional matrix 
illustrated in Figure 3 must pass through either the point 
(i-l,j), or (i-l,j-l), or (i,j-l). The minimum accumulated 
distance to point (i,j) is then given by: 

D( i, j)=Dist( i , j )+Min(D( i-1 , j ) ,D( i-1 , j-1 ) ,D( i , j- 
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FIGURE 3 * DYNAMIC TIME ALIGNMENT 

where Dist (i,j) is the distance between the reference and 
training pattern at tine j. The algorithm recursively 
computes this distance column by column to determine the 
minimum accumulated distance to the point (M,N), where M is 
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the number of frames in the reference and N is the number of 
frames in the unknown (training). This path results in a time 
alignment in which the reference word has the maximum 
acoustic similarity with the input. 

B. Artificial Neural Network [10-18] 

Artificial neural network is a non-algor i thmic 
information processing structure based on the architecture of 
our biological nervous system. The structure is composed of 
a massive number of processing elements operating in a 
predetermined parallel operation. The processing elements 
are connected by links with variable weight. The topology of 
each network determines the way each processor is connected 
to the other. The link can be excitory, inhibitory, or has 
no effect on the activity of the processing element. See 
Figure 4. The primary processing at each element consists of 

Type I Type R 



the calculation of weighted sums of the form f(W i{ ,X|) and 
weight changes of the form W ij new= ,X A ,Xj . . . ) . the 
function f is usually a nonlinear function, Wjj is the weight 
of the link from element i to element j and X^ is the input 
to element i. See Figure 5. Figure 6 shows some of the most 
popular neural 
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net models that can be used as classifiers with the classical 
algorithms that are most similar to the neural net models are 
listed along the bottom! 10]. As shown in this figure, these 
nets are first divided according to whether the input is 
binary or cont inuous-value . Second, they are divided 
according to whether they need supervision during training or 
not. Since a multi-layer perceptron backpropagation model is 
implemented in our study, we discuss in the following section 
the topology and the learning rule of that model. 

The Backpropagation Model [10,14,15] 

Figure 7 shows the topology of the backpropagation 

Output Patterns 



Input Patterns 

FIGURE 7. BACKPROPAGATION NETWORK. 


6-10 


model. The model consists of 3 layers (input, hidden, and 
output). each node in the input layer is connected to every 
node in the hidden layer also each node in the hidden layer 
is connected to every node in the output layer. The input of 
the model is a continuous valued vector x_, Xp .Xu ^ 
representing either the LPC coefficients or tne frequencies 
of the Formants which are calculated from the coefficients. 
Investigations, mainly experimental, will be carried out to 
see which input is most suitable for word recognition* The 
actual output y Q , yj,...y M _i is calculated as: 

y = f( ^ w ij x j ~ ®) where 9 is a predetermined threshold 

and the function f is the sigmoid nonlinearity: 


f (a) = 


1 + e 


-(a-9) 


Training of the network is carried out by setting the 
output of the model to the desired output vector d Q , 
d^,...d^_^. All elements of the desired output vector are 
set to zero except for that corresponding to the current 
input training word which is set to 1. The weights are 
adjusted recursively from the output nodes to the hidden 
nodes by the formula: 


5 »ij old * » Vi 

where u is the gain factor, is the error, and X’^ is 
either the output of node i or i.8 an input. If node j is an 
output node , then 


6 j = yjd-yjHdj-yj), 

where dj is the desired output of node j and y* is the 
actual output. If node j is an internal hidden node, then 

6 d = xj’U-Xj’) A k w jk , 


where k is over all nodes in the layers above node j 
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ABSTRACT 


An environmental health monitoring plan for Space Station 
will ensure crew health during prolonged habitation. The Space 
Station, "Freedom" will operate for extended periods, 90+ days, 
without resupply. A regenerative, closed loop life support system 
will be utilized in order to minimize resupply logistics and 
costs. Overboard disposal of wastes and venting of gases to space 
will be minimal. All waste materials will be treated and 
recycled. The concentrated wastes will be stabilized and stored 
for ground disposal. The expected useful life of the station 
(decades) and the diversity of materials brought aboard for 
experimental or manufacturing purposes, increases the likelihood 
of cabin contamination. Processes by which cabin contamination 
can occur include: biological waste production, material off- 
gassing, process leakage, accidental containment breach, and 
accumulation due to poor removal efficiencies of the purification 
units. 

An industrial hygiene approach was taken to rationalize 
monitoring needs and to identify the substances likely to be 
present, the amount, and their hazard. This requires a thorough 
knowledge of the onboard processes, their products and by- 
products. Many factors influence the monitoring requirements for 
Space Station: the enclosed space, the recirculation of supply 
air, the experiences of past missions, the unique experimental 
and manufacturing facilities, and the interfacing of other 
modules with the U.S. core modules. Monitor development and 
selection will be complicated due to the many technologies 
competing for the life support systems, and the number of 
experimental payloads under development, each having their own 
unique monitoring requirements. 

Monitoring options include: on-line sensors for process 
control and determining efficacy of the life support regeneration 
and purification units; monitors for specific components and 
contaminants critical for life support; monitors for surrogate 
parameters representative of contaminant groups likely to be 
encountered; broad spectrum analyzers capable of identifying and 
quantifying nearly any contaminant for leak detection and 
remedial action; manual aboard sampling and analysis; and sample 
collection/preservation with ground based analysis. Other 
factors considered in the monitoring plan include: results of on 
ground system tests; consequences of contaminant detection 
failure; instrument parameters such as weight, volume, 
reliability, specificity, detection limit, and maintenance 
requirements; crew time and effort; expendables; and waste 
production. 
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INDUSTRIAL HYGIENE PERSPECTIVE 


It is standard practice in industry to monitor for those 
agents in the workplace that may occur at levels approaching or 
exceeding safe limits. Defining those limits and deciding on 
which agents to monitor, calls for a systematic and thorough 
approach, often employed in the field of industrial hygiene. 

Industrial Hygiene is the science and art devoted to the 
anticipation, recognition, evaluation, and control of 
environmental factors and stresses arising in or from the 
workplace that may cause sickness, impair health and well being, 
discomfort or inefficiency among workers, their families or 
members of the community at large. * 

Work hazards should be anticipated before they exist and 
action must be taken to prevent their occurrence, so no 
individual is placed at undue risk. Recognition of potential 
hazards requires familarity with the process and work operations 
involved, the maintenance of an inventory of agents, a periodic 
review of the job activities, and the effectiveness of control 
measures . 

Evaluation principally involves performing the monitoring to 
determine the exposures, comparing the results with standards, 
and communicating the judgement on the: degree of hazard with the 
individuals affected and those in authority to take corrective 
action. Monitoring is a continuing program of observation and 
judgement. Reasons for monitoring are: to determine exposure 
levels; to determine effectiveness of control measures; to detect 
process changes; to investigate complaints; and to confirm 
compliance with standards. 

In monitoring, part of the environment is sampled and the 
quality of the whole environment is inferred. A number of factors 
must be considered in order to take a representative sample. In 
industry the concern is worker exposure, so air samples are taken 
in the breathing zone of the worker, specified as twelve inches 
from the nose. Data from such personal samples are better 
correlated with exposures than any other type because of large 
concentration gradients and awkward operator/machine interfaces 
that can occur in industry. Sampling general workroom air is 
effective when the emission rates are uniform and there is good 
mixing. However, when relying on area sampling, documentation of 
its correlation with personal samples must be available. Sampling 
at the operation usually yields the highest concentration and 
worst case exposures. Rest areas should be occasionally monitored 
to confirm a clean zone of zero exposure. 

All individuals exposed at or above the action level must be 
monitored. The action level is some fraction of the standard 
below which exposures must remain, usually one-half or one-fourth 
the standard. It is prudent to monitor nearby workers and 
complainers, as they may be sensitive to the agent, or an 
unidentified process leak may be present, or a control device may 
have failed. 
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Sample duration depends on the detection limit of the 
technique used for analysis, the standard, and the estimated air 
concentration. The sample duration should represent an 
identifiable period, consistent with the standard, the work shift 
and the process. Common sample durations are eight hours, for 
time weighted average standards and ten minutes for peak 
concentrations. Continuous monitoring is required if ceiling 
values are likely to be exceeded or when conditions immediately 
dangerous to life and health may occur. 

Monitoring results are compared with standards to determine 
the severity of the exposure and a plan of action. Standards 
which must be met are the Permissible Exposure Limits (PEL'S) 
published by the Occupational Safety and Health Administration 
(OSHA). 2 These are legal standards and represent the minimum 
effort required to maintain a safe and healthful workplace. Other 
standards or guidelines are recommended by professional 
organizations and consensus groups. Most widely used are the 
Threshold Limit Values (TLV's) of the American Conference of 
Governmental Industrial Hygienists (ACGIH) . 3 Most PEL'S were 
adopted from TLV's and most are eight hour time weighted average 
or ceiling values. For each standard there is a criteria document 
which presents data on which the standard is based. 

Control can be achieved by reducing the emission source, 
interrupting the air path, and protecting the receiver. Specific 
control methods are: elimination or substitution with a less 
harmful agent? process change; enclosure of the process or 
source; isolation in time or distance; wet methods to reduce dust 
loadings and to scrub gases; local exhaust ventilation to remove 
the . contaminant at its source; dilution ventilation? 
housekeeping; adequate maintenance; training and education of 
workers; area and personal monitoring; and use of personal 
protective equipment, mainly as backup for harmful or life 
threatening situations. 

Elimination, substitution, and process change are ground 
based decisions that must be made well before flight, strict 
flight requirements are being relaxed to allow the use of off the 
shelf items for Space Station. Isolation, wet methods, and 
dilution are not compatible with space vehicles. During space 
flight, control methods will rely heavily on enclosure, venting 
to the trace contaminant control system, housekeeping, and 
monitoring to assure a healthy cabin environment. 

Supplied Breathing Air 

On ground, when we encounter situations with probabilities 
of oxygen deficiencies or conditions immediately dangerous to 
life and health, breathing air must be supplied. 4 It must be 
Grade D or better, with specifications as follows: oxygen 19-23%, 
carbon monoxide (CO) £20 ppm, carbon dioxide (C0 2 ) <1000 ppm, 
condensed hydrocarbons <5 mg/M 3 , and the water content must be 
stated. The air source is ambient air, supplied by a compressor 
to a delivery system or cylinder for storage. Any one of these 
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can be a source of contamination. A frequently encountered cause 
of morbidity from contaminated air is carbon monoxide, from the 
incomplete oxidation of hydrocarbons of the compressor fuel, 
exhaust, and lubricants. Another source is from over heating of 
charcoal filters used for purifying the intake air. For this 
reason, continuous carbon monoxide and/or temperature monitors 
are required on air supply compressors. The condensed 
hydrocarbons can cause lipoid pneumonia and decrease the gas 
exchange membrane surface in the lung. They may also oxidize to 
CO and C0 2 within the storage cylinders. The water content is 
important: if excessive, it can cause regulator valves to clog or 
freeze, and promote cylinder corrosion; and if the air is too dry 
then irritation of mucous membranes, eyes, nose, and throat can 
occur. The oxygen (0 2 ) content should be checked routinely to 
confirm adequate concentrations. Industrial users check every 
cylinder of air before use, since several fatal incidents have 
occurred due to low oxygen concentrations. 

Confined Spaces 

In confined spaces, an area is enclosed or partially 
enclosed with poor ventilation and mixing with the outside air is 
limited. Before entry, the atmosphere should be tested to ensure 
the 0 2 concentration is 19.5-21%. Toxic, flammable, and oxygen 
displacing gases and vapors should silso be monitored. Common 
contaminants to check are hydrogen sulfide (H 2 S) , CO, C0 2 , and 
methane (CH 4 ) , plus any other material that is likely to be 
present because of prior storage. If conditions immediately 
dangerous to life and health are probable, then continuous 
monitoring is required, since conditions within a confined space 
can change rapidly. 5 

Nuclear Submarines 

Submarines are enclosed environments with controlled 
atmospheres and are capable of remaining submerged for extended 
periods. Each submarine is equipped with a monitoring system 
designed to measure gaseous atmospheric constituents which are 
important in life support, called the Central Atmosphere 
Monitoring System (CAMS-I) . It continuously analyzes air from the 
main fan room for eight substances, with a mass spectrometer (MS) 
and an infrared analyzer (IR) . The mass spectrometer looks at 
specific m/e in the range of 2-300 amu. The average failure time 
was reported to be 3500 hours.® It is; considered a monitor and 
not a trace contaminant detector, since sensitivity was 
restricted to meet stability and dependability requirements. 
CAMS-II is under development. 7 It will use a scanning mass 
spectrometer to measure 12 substances continuously from multiple 
sample ports located throughout the submarine. The MS will 
measure : 
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non-methane aliphatic hydrocarbons 

0-100 

ppm 

aromatic hydrocarbons 

0-10 

ppm 

benzene 

0-10 

ppm 

carbon dioxide 

0-3.3 

% 

hydrogen 

0-5 

% 

nitrogen 

0-80 

% 

oxygen 

0-25 

% 

refrigerant R-ll 

0-1200 

ppm 

refrigerant R-12 

0-1200 

ppm 

trichloroethylene 

0-100 

ppm 

water vapor 

0-4 

% 


Since the MS cannot resolve CO and N 2 , as both have m/e of 28, an 
IR measures for CO. CO is a major concern in submarines because 
of combustion processes, smoking, cooking, and smoldering of 
activated charcoal filters designed to remove odors, 
refrigerants, and hydrocarbons. Refrigerants themselves are of 
little concern, however their decomposition products from 
compressors or fire are corrosive, toxic, and will poison the 
catalytic oxidizer. Also on board is a paramagnetic 0 2 analyzer 
and a Dwyer C0 2 analyzer. A portable photoionization detector is 
used for hydrocarbons, it is referred to as the "trace gas 
analyzer." Earlier CAMS used a gas chromatograph with a flame 
ionization detector as a hydrocarbon detector in which 100 ml 
samples were injected. It was not compatible with submarine 
duty. 8 An assortment of Draeger detector tubes are used for leak 
detection and for CAMS backup. 

MONITORING ON SHUTTLE 

The air revitalization system and many other systems on the 
Shuttle are different from those proposed for Space Station. 9 On 
the Shuttle, C0 2 is removed by LiOH canisters which are changed 
when the C0 2 concentration is 50 mm of Hg, or every 12 hours. 
Activated charcoal filters are used to remove hydrocarbons and 
odors. Carbon monoxide is oxidized to C0 2 in a low temperature 
catalytic oxidizer, located downstream from the humidity and 
thermal control unit. The activated charcoal filters and the low 
temperature catalytic oxidizer are capable of removing most 
contaminants that may occur. Fresh air is supplied from cryogenic 
liquid nitrogen (N 2 ) and (0 2 ) stored on board. Thus on the 
Shuttle, continuous monitoring is done for total pressure, 
temperature, humidity, 0 2 , and C0 2 . Air samples are taken in 
evacuated bottles, and by active collection on various sorbent 
tubes, such as charcoal, tenax, and molecular sieve. Those 
samples are post-flight analyzed in ground based laboratories. 
The charcoal and LiOH air purification filters are also often 
analyzed for contaminants. 
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MONITORING NEEDS 


Experiences of past missions and ground based systems tests 
have identified a number of health concerns that should be 
addressed in a monitoring plan for Space Station. Paramount is 
the flight and post flight health complaints of the crews: 
headache; irritation of the eyes and upper respiratory tract; and 
odor complaints, symptomatic of noxious air .- 1 - 0 Early missions 
had insufficient monitoring data for evaluation, which indicated 
a need for a more comprehensive monitoring system. Analyses of 
activated carbon and lithium hydroxide filters of the atmospheric 
revitalization systems, and the active sampling and analysis for 
air contaminants of later missions have identified over 250 
contaminants in spacecraft air . 11 Most were observed at trace 
levels, well below the Spacecraft Maximum Allowable Concentration 
(SMAC) . Others may have elicited symptoms among crew members, may 
accumulate to harmful levels, or may have potential to poison the 
spacecraft life support system. 

Nitrogen tetroxide (N 2 O 4 ) , hydrazine, and monomethyl 
hydrazine are the main liquid propellants to be used on Space 
Station. Because of the quantities involved and the frequency of 
extra vehicular activity (EVA) , some Space Station contamination 
will occur. An air lock will likely serve as a decontamination 
station and will contain a propellant monitor. If elevated 
propellant concentrations are detected in the air lock, then that 
atmosphere will be dumped to space to prevent contamination of 
the cabin atmosphere. The air revitalization and trace 
contaminant control systems were not designed to handle high 
pollutant loads. Some N 2 0 4 contamineition occurred on Apollo- 
Soyuz . 10 

Halon 1301 is the fire suppressant to be used on Space 
Station. Halon was detected on spacelab mission SL-1 and on 
Shuttle missions STS-3, and STS-4 . The trace contaminant control 
system (TCCS) will only handle modest quantities. Halon 
degradation products are toxic and will poison the catalytic 
oxidizer. If a halon release occurs it will be necessary to vent 
the cabin air to space and repressurize. 

Methane is a metabolic product which usually accumulates as 
each mission progresses. It will likely be the contaminant of 
greatest concentration. The Bosch C0 2 reduction system, a 
candidate for the air revitalization system (ARS) , will produce 
large quantities of methane. A high temperature catalytic 
oxidizer will be required to keep CH 4 concentrations below 1 
ppm. 12 • 13 

CO, a product of incomplete combustion, may be released from 
metabolic processes, smoldering of carbon filters, or fire. The 
Bosch C0 2 reduction system produces CO and the potential for 
rapid accumulation exists, if not removed by the trace 
contaminant control system . 12 • 13 

Ammonia (NH 3 ) , a product of metabolism will be released from 
urine processing, and it is probably a degradation product of the 
solid amine resin proposed for the ARS . 14 Phosphoric acid 
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impregnated charcoal filters can remove NH 3 . 

Hydrogen (H 2 ) will be produced by electrolysis and used in 
C0 2 reduction by both the Bosch and the Sabatier 
processes . 12 • 13 • 15 A pressure gradient will be used to minimize 
the likelihood of explosive mixtures from developing, if a leak 
occurs . 

Toluene was detected on a number of missions. On Shuttle 
mission STS - 2 , toluene approached the SMAC value in one sample. 
Subsequent analyses indicated that for the sample, the additive 
toxicity hazard index for systemic poisons was exceeded by 1.22 
times, with toluene the major constituent . 10 Toluene is also a 
contaminant which off gases from the solid amine resin of the 
ARS . 14 

Trimethylamine is a principal breakdown product of the solid 
amine resin of the ARS. The trimethylamine concentration has 
exceeded safe limits in tests of the ARS . 14 Because of the 
numerous trace organics off gassing from solid amine process a 
post sorbent bed such as phosphoric acid impregnated charcoal 
will be used. 

Glutaraldehyde and silicon escaped containment on Spacelab 
mission SL-D1. Glutaraldehyde is a preservative and disinfectant 
with irritating properties. It may also be used in 
electrophoresis experiments on Space Station. Silicon compounds 
are catalyst poisons and will occur on Space Station. 

Freons have been detected on all Shuttle missions . 16 The 
degradation products are corrosive, irritating, toxic, and 
catalyst poisons. Freon 12 will be on Space Station. 

A computer model developed from Shuttle charcoal canister 
analysis for TCCS contaminant removal studies indicated that five 
contaminants may exceed SMAC values: propenal (acrolein) , an 
irritant; benzene, a systemic poison and carcinogen; o- 
diethylphthalate, an irritant; propyl fluorosilane, an irritant 
and catalyst poison; and 2 -methylhexane, a central nervous system 
depressant . 17 Benzene has also temporarily exceeded SHAC values 
during preflight off gassing tests . 16 

Ethanal (acetaldehyde) , ethanol, dichloromethane , and 
acetone have a high frequency of occurrence on shuttle missions 
and are likely to be present on Space Station . 16 

Oxidation products will be produced in the catalytic 
oxidizer. Post sorbent beds are necessary to prevent the release 
of oxidants and free radicals to the cabin air from the TCCS. 
Also, it has been hypothesized that secondary pollutants are 
important in cabin atmospheres. Trial simulations have indicated 
that spacecraft cabins may develop elevated N0 2 concentrations 
and ozone (0 3 ) concentrations exceeding SMAC values . 18 Oxidation 
products, N0 2 , O 3 , and formaldehyde, were among the contaminants 
suspected of causing irritation on Shuttle flights, although 
particulates from biological sources were the undisputed cause of 
crew discomfort . 19 

Major metabolic products which must be removed by the ARS 
and the TCCS are CO, co 2 , NH 3 , H 2 S, CH4, organic acids, and 
mercaptans . 2 0 
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Foul odors have been observed on a number of missions. Many 
were attributable to metabolic products. An unusual odor and crew 
headaches occurred on Shuttle flight STS -6 . Burnt wire insulation 
from an electrical short was the suspected causal agent. 10 
Electrical fire can produce a number of noxious agents including 
halogenated organics, benzene derivatives, nitriles, and 
cyanates. 21 Space Station design must be able to handle such 
contingencies either through the TCCS or a smoke removal unit, 22 
without having to rely on venting the cabin air to space and 
repressur i z ing . 

For fire safety concerns, Halon 1301 will be used for fire 
suppression, followed by venting cabin air to space and 
repressurizing. Smoke detectors are an integral part of the fire 
detection and suppression system. To protect from toxic 
combustion products, infrared monitors are recommended for CO, 
hydrogen fluoride (HF) , and hydrogen cyanide (HCN) . 22 

Volatiles will be released to the atmosphere from 
electrolysis and from phase change urine processing. Carboxylic 
acids and phenols will be major contaminants. 24 ' 25 Iodination 
products from the water disinfection process may cross the 
air/water interface and permeate the entire life support 
environment. The identity of these products, their expected 
concentrations, and their medical effects are largely unknown. 26 
However, I suspect the byproduct concentrations and effects of 
iodination are less than those resulting from chlorination. 

SPACE STATION CONFIGURATION 

Space Station is designed to operate for extended periods, 
90 plus days, without resupply. A regenerative, nearly closed 
loop life support system will be required to minimize resupply 
logistics and costs. Overboard disposal of wastes and venting of 
gases to space will be minimal. All waste materials will be 
treated and recycled. The concentrated wastes will be stabilized 
and stored for ground disposal. The expected useful life of the 
station is decades and a diversity of materials will be brought 
aboard for experimental or manufacturing purposes. The likelihood 
of cabin contamination is great. Cabin contamination can occur 
from a number of sources: biological waste production, material 
off gassing, process leakage, accidental containment breach, and 
accumulation due to poor removal efficiencies of the purification 
units . 

The Space Station, "Freedom," will have four modules: the 
U.S. Laboratory (USLAB) ; the U.S. Habitation module (USHAB) ; the 
Japanese Experimental Module (JEM) ; and the European Space Agency 
(ESA) module, Columbus. The modules are connected by four 
resource nodes. Two airlocks and a logistics module are connected 
to the resource nodes. Each module will have an independent 
Environmental Control Life Support System (ECLSS) , complete with 
a Trace Contaminant Control System (TCCS). The U.S. modules will 
have four Air Revitalization Systems (ARS) , two in each module. 
Each ARS is designed to support four crew members; One ARS at a 
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time will operate in each module. 

The ARS will provide ventilation to each module and node but 
not to the airlocks. Intramodule circulation will approximate 
near perfect mixing. Intermodule air exchange design is 130 cubic 
feet per minute (CFM) through 4-4.5 inch ducts . 27 The 
ventilation design is primarily based on: heat transfer and 

humidity control to maintain crew comfort; and 0 2 supply and C0 2 
removal requirements to maintain a healthful atmosphere . 28 The 
air exchange rate will be 1-2 years, achieved through air loss 
from leakage and airlock EVA. 

ESA Columbus Module 

ESA Columbus monitoring requirements are based on the fact, 
that the types of contaminants and their buildup characteristics 
are not precisely known, and that safety will require the 
monitoring of all contaminants permanently . 29 Hence, ESA is 
developing, for Columbus, a GC/MS to monitor N 2 , C0 2 , 0 2 , H 2 0, 
plus a lengthy target list of trace contaminants. The contaminant 
list includes 15 alcohols, 5 aldehydes, 12 aromatics, 11 esters, 
3 ethers, 21 halocarbons, 37 hydrocarbons, 8 ketones, and 11 
miscellaneous compounds , such as NH 3 , CO, H 2 S, H 2 , O 3 , S0 2 , NO, 
and N0 2 • 

Japanese Experimental Module 

The 70 M 3 JEM and its 24 M 3 logistics module are designed 
for two crew members . 30 The JEM will have an independent ECLSS 
interfaced with the Space Station core. However, concentrated C0 2 
will be returned to U.S. modules for processing, and wastewater 
will be returned to the U.S. modules for processing and H 2 0 
recovery. Its TCCS design is based on a contaminant load of 15.4 
grams/day . 31 

The JEM will continuously monitor for total pressure, 
temperature, humidity, 0 2 , and C0 2 . The JEM will rely on the U.S. 
modules for trace gas analysis via sample lines. A GC/MS is 
expected. 

U.S. Modules 

The U.S modules will provide facilities for on-orbit 
repair, health maintenance, and a number of material processing 
and biological experiments intended to lead to manufacturing in 
space . 

A maintenance work station will allow on-orbit repair of 
defective or damaged hardware. Processes likely to be required 
are drilling, sawing, welding, soldering, and epoxy gluing. A 
work bench/contaminant control console is envisioned that will 
collect the particulate and gaseous emissions generated in the 
repair process near their source . 32 The rack would be equipped 
with filters and the air recirculated with some venting to the 
TCCS. The work station would be a source of particulates, metal 
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fumes , and gases not encountered on prior missions. 

The health maintenance facility will provide critical care 
for one individual for 28 days and outpatient care for the crew 
complement for the mission duration. The equipment and supply 
list for this facility will be lengthy. 33 It may be an 
additional source of trace contaminants, mainly sterilants. 

The U.S. Laboratory will provide facilities for experiments 
and manufacturing. 34 The on board processes will generate 
biologicals, combustion and oxidation products, acid gases, metal 
and crystal fumes, and assorted lab wastes. Many of these 
materials are capable of adversely affecting the ECLSS subsystems 
by poisoning the catalyst or absorption beds, or they could 
appear in the humidity condensate, the potable water supply. 
Materials will have to be stored, then transported to the point 
of use, and the waste products handled. The lab racks will be 
contained with at least a two failure tolerant design. They will 
be equipped with some type of contaminant control equipment and 
vented to the TCCS. The lab racks should be equipped with 
monitors, specific for the process they contain to detect 
internal leaks. The chemical storage area should be monitored, 
and the cabin atmosphere must be sampled to alert the crew of any 
leak. 

Trace Contaminant Control System 

The technology base for the TCCS is good. Only limited 
system tests have been conducted but they have worked as 
predicted. The TCCS will consist of fixed bed charcoal filters, 
high -efficiency particulate filters, and a high temperature (680 
*C) catalytic oxidizer (palladium/alurainum) with pre and post 
sorbent beds of LiOH. There will be four units, two in each 
module. The air flow through each catalytic oxidizer is 2.5 CFM, 
or 5 CFM for the two U.S. modules. 12 This is only 4 air changes 
per day of what should be considered as fresh air. It may be too 
low. The TCCS will receive cabin air from the temperature and 
humidity control system. Purge gases from the ARS, waste water 
recovery, urine processing, waste reduction and storage systems, 
and lab racks are to be routed to the TCCS for contaminant 
removal. For comparison, the indoor air quality ventilation 
guideline is 15 CFM per person. 35 The guideline is intended to 
keep odors to an acceptable level to 80% of the visitors entering 
the space and it assumes that one third of the occupants are 
smoking at the rate of 2.2 cigarettes per hour. 

TECHNOLOGY ASSESSMENT 

A technology assessment study on monitoring systems was 
performed by Battelle Columbus Division for NASA. 11 They 
recommended: 1) a long path Fourier transform infrared (FTIR) 
analyzer for rapid detection of high risk contamination 
incidents, and 2) a GC/MS with mass selective or ion-trap 
technologies for detailed monitoring of extended crew exposures 
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to low (ppb) contaminations. Priority requirements for rating the 
candidate monitoring systems were: real time output, which is 
particularly important in closed environments and long duration 
missions; ability to detect and quantify a wide range and number 
of volatile compounds; ability to trigger a warning when the SMAC 
is approached; and the ability to monitor several modules and 
airlocks simultaneously. In the assessment, many instrument 
parameters were considered, each weighted equally: weight, 
volume, power requirements, sensitivity, dynamic range, response 
time, selectivity, growth capability, crew time, by-product 
generation, consumables, reliability, and maintenance. 
Instruments and systems considered included: CAMS, long path 
FTIR, matrix isolation FTIR, GC/MS, GC/photoionization detector, 
electrochemical devices, and chemiluminescent monitors. Long path 
FTIR and GC/MS appeared to be the most promising technologies for 
Space Station monitoring. However, several critical questions are 
yet to be answered: what is the target list of compounds to be 
monitored; which must be monitored on a continuous real time 
basis; at what concentration ranges; and in how many locations? 

CONCLUSIONS 

Development of the monitoring plan for Space Station will be 
a continuing process. The monitoring system must be adaptable to 
accommodate new parameters and concentration ranges. All agents 
should be monitored that have a reasonable probability of 
occurrence at or above some action level, such one-half the SMAC. 
The analytical method relied upon must be able to quantify at 
action level concentrations. The basis for monitoring should be 
the contaminants: toxicity, quantities or production rates, 
removal efficiencies of the ECLSS system, and capacity to poison 
the ECLSS system. 

Monitoring for all contaminants at multiple sites and at 
part per billion concentrations is impracticable and should not 
be attempted. Monitoring for a contaminant should not be done 
simply because there is an assigned SMAC. Equipment control 
monitors , and the monitoring of a surrogate parameter as a 
substitute for the etiological agent should not be relied upon 
when making health evaluations. 

A GC/MS should be developed for Space Station. It should be 
a fixed instrument that continuously samples the well mixed cabin 
return air. It should have a sample line to each module and 
airlock for routine comparison of atmospheres from remote 
sections of the spacecraft. It should be capable of monitoring 
for major atmospheric components (e.g. N 2 , 0 2 , C0 2 , H 2 0, H 2 , CO, 
CH 4 , and Halon 1301) and numerous trace gases (e.g. hydrocarbons, 
halocarbons, silanes, and phenols) . 

A LP/FTIR could detect modest levels of many compounds not 
readily amendable to GC/MS analyses (e.g. inorganics, acid and 
basic gases, oxidation products, and CO) . It could also serve as 
a trace hydrocarbon detector. 

A portable hand held hydrocarbon detector should be 
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available for leak detection in remote areas of the spacecraft. A 
TLV Sniffer (has good response for CH4) or H-Nu photo ionization 
type instrument can detect low ppm concentrations of hydrocarbons 
and should be adequate for this purpose. 

Redundant monitors should be present for critical 
parameters: IR for CO; IR or electrochemical for CO 2 ; and 
electrochemical for 0 2 . 

Particulates and fumes could be measured by either 
photometric or piezoelectric aerosol monitors. Sample collection 
and ground based analyses will be necessary for verifying mass 
loadings, particle size distributions, and chemical components 
often present particle bound (metals, semi-volatile and non- 
volatile compounds) . 

Each experiment and manufacturing process must be evaluated 
for possible sources of cabin contamination. The lab racks should 
be equipped with monitoring devices specific to the process being 
contained. There are many miniature photometric, IR, and 
electrochemical devices which could serve this need. 

Finally, sample collection and preservation will have to be 
continued for ground based analyses, to confirm the accuracy and 
reliability of the on board monitoring system. 
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ABSTRACT 


Several radiation models for low density nonequilibrium 
hypersonic flow are studied in this report. It is proposed that 
these models should be tested by the 3-D VRFL code developed at 
NASA/JSC. A modified and optimized radiation, model may be 
obtained from the testing. Then, the current VRFL code could be 
expanded to solve hypersonic flow problems with nonequilibrium 
therma 1 radiation. 

INTRODUCTION 


When hypersonic flight vehicles, such as Aeroassisted Qrbiter 
Transfer Vehicles (AOTVS), traveling in low density env 1 r o nmen t , 
strong bow shocks occur adjacent to the vehicles. Behind the 
shock front, the flow involved could be at thermal and chemical 
nonequilibrium. At the extreme high temperatures, the effect of 
nonequilibrium gas radiation may be prominent. Determining the 
magnitude of radiation in this environment is considered to be one 
of the most important problems associated with AOTVS (1-3). 

There are many research works dealing with thermaL and 
chemical honequ i 1 i br ium hypersonic flow, some of the recent papers 
are given in the references (4-10). Although the nonequilibrium 
radiation is part of the nonequilibrium precesses, it is usually 
assumed as equilibrium radiation or neglected completely in the 
analysis. The relatively few work in nonequilibrium radiation 
research will be mentioned in the later sections. 

The purpose of this report is to give a brief discussion of 
nonequilibrium radiation in low density hypersonic flow, with 
various models for radiative interaction, numerical solution 
methods, and some of the available codes. In the next section, 
the solution methods used for the general nonequilibrium 
hypersonic flow will be discussed, this is followed by the section 
dealing with radiation. 

THERMAL AND CHEMICAL NONEQU I L I BR I UM 
HYPERSONIC FLOW 

There are several numerical methods which can be used for the 
nonequilibrium hypersonic flow (11). The Direct Simulation Monte 
Carlo method (12,13), based on probab i 1 ist ic collision model, 
requires excessive amount of computer time to solve a two - 
dimensional flow problem so it may be not able to solve a three - 
dimensional flow problem. The Viscous shock Layer method (14, 
15), which used the simplified conservation equations in the shock 
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layer, fails to predict the flow for the whole vehicle* To solve 
the complete flow field around a hypersonic vehicle, the full 
Navier stokes equations should be used,, 

Full Navier-Stokes Equations 

In order to calculate radiation from the nonequi 1 ibr ium 
region of the flow, one needs to know the thermal and chemical 
state of th gas in the region- This includes the nonequ i 1 i br i urn 
species concentrations and various temperatures to express the 
different energy modes of the species- These state variables can 
be calculated by simultaneously solving a set of mass, momentum, 
and energy conservation equations, the so— called full Navier - 
Stoke equat ions (2,16,17). The equations consist of one mass 
conservation equation for each species, one overall momentum 
equation, and several energy equations. The number of energy 
equations should equal to the number of temperatures which 
characterize the various energy modes of the species. For 
example, when using a two-temper a ture model of one electron 
temperature and one heavy-par t i c I es temperature, the two energy 
equations could be one electron energy conservation equation and 
one overall energy conservation equation. If the two- temper a ture 
model used involving one electron-vibrational tempera ture and one 
rotational-translational temperature, the two energy equations 
could be one e 1 ec ton-v i brat i ona 1 energy conservation equation and 
one overall energy conservation equation. For a three- tempera ture 
model of an electron temperature, a vibrational temperature for 
all species, and a t r ansi a t i ona 1 -r o t at i ona 1 temperature for all 
species, then there are three energy conservation equations for 
them. In the literature, the thermal nonequ i 1 i br i urn model could 
involve more than three temperatures (18), then the computational 
time has to increase greatly. 

In the overall energy equation, there is a term represents 
the contribution of thermal radiation. This term is expressed by 
a divergence of a radiation heat flux vector. To evaluate this 
radiation heat flux vector, a radiative transport equation is 
needed. The coupling of these equations makes the nonequ i 1 i br i urn 
flow problem very difficult to solve. The radiative transport 
equation is discussed in the later section. 

Rate Equations 

To formulate the full Navier-Stokes equations for 
nonequilibrium flow, various rate expressions for source functions 
should be specified. In the mass conservation equation for each 
species, the source term is the mass production rate of species 
due to chemical reactions. The chemical reaction rate is a 
function of the various temperatures used for species. In the 
energy equation involving translational - vibrational energy 
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exchange, the relaxation time in the rate equation is of the 
Landau - Teller type <2, 17). For* the el.ectron-vibrat ional energy 
exchange, the relaxation time expression is a modification of the 
Landau-Teller expression (2,19). 

NONEQUILIBRIUM GAS RADIATION 


As mentioned in the earlier section, to include a radiation 
energy term in the overall energy equation, an additional equation 
is needed, this is the equation of radiative transfer. The 
coupling between this radiative transfer equation and the energy 
equation in gas dynamics is of a complicated i n tegr a 1 -d i f f er en t i a l 
form. When nonequilibrium requirement is included, it adds more 
complexity to the problem. In the following sections, the basic 
equation, various simplication and solution methods are briefly 
d i scussed . 

The Equation of Radiative Transfer 

The equation of radiative transfer is a continuity equation 
for the number density of photons in a specific solid angle and at 
a specific frequency (17). The equation can be formulated in 
terms of intensity, emission and absorption coefficients. In a 
nonequilibrium environment, a knowledge of the emission and 
absorption coefficients is a requisite to a careful solution of 
the equation of transfer. The discussion of various spectral 
models for emission and absorption coefficients is in the next 
section. since the radiation transport is coupled with flow 
equations, for a formal solution procedure, one must rely on an 
iterative method to solve the system of flow equations with 
radiation tr anspor t . 

Radiation Transport Models 

There are a number of spectral models for radiative 
interaction computations in the literature. These models are 
briefly outlined in the following: 

Line By Line Transport Calculation Park, etc (3,20,21) has 
developed a model to calculate radiative properties for 
nonequilibrium air. the model assumes that the radiation 
transport is not strongly coupled with the flow equations, so the 
heavy-particle translational temperature, vibrational teniperature, 
electron temperature, and species mole fractions are solved 
independently from flow equations. These items will be used as 
inputs to calculate the number densities of internal state of 
atoms and molecules by using the so called quasi steady-state 
assumption. A detailed line-by-line calculation is performed to 
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evaluate emission and absorption characteristics independently, 
using the calculated, nonequilibrium, excited internal state 
density values. This model requires a tremendous amount of 
computer time and memory space for the required computations. It 
becomes even more so when four vibrational temperatures were used 
for the gas (18). It is very difficult to apply this model to a 
general 3-D hypersonic flow problem. 

The Step-Function Absorption Coefficient Model A step function 
model has been developed to compute the radiation absorption 
coefficients by Zoby, Sutton, etc. (22). In general, the model 
includes, free-free, bound-free, and bound-bound transitions of 
various atoms and molecules. The inputs to the model are species 
niimber density and temperature of the flow field. The frequency 
dependence of the absorption coefficient is represented by a set 
of steps, these steps has a fixed, but not necessarily equal 
widths chosen to resolve the detailed spectra adequately. The 
total absorption coefficient of a particular step is a summation 
of all the average absorption coefficient of individual energy 
transitions in that step. The average absorption coefficient of 
an individual transition can be obtained by integration the 
absorption coefficient over frequency in that step. When the step 
width is properly chosen, there should be analytic expressions to 
approximate the absorption coefficient in that step for the 
purpose of integration. Thus the frequency dependence for each 
step- becomes a set value, and the total absorption coefficient for 
that step can now be expressed in terms of species number density 
and temper a ture . When this is done, the radiative transfer 
equation is well expressed in each frequency step. In Reference 
EE, a 58 steps were used to cover the frequency range from 0 to 
1 7 ev . 


Although this model is mainly applied to equilibrium 
radiation, it includes the frequency dependent absorption 
coefficient (22,23,24). In addition, with some mod i f i ca t i ons , 
this method can be applied to nonequilibrium radiation (25). 

There is a similiar but simpler approximation called the band 
approximation model for absorption coefficient (26,27). This 
model is also for equilibrium raidation and usually is a two-band 
approximation model . 

Parametric R a diation Model The nonequilibrium radiation 
environment corresponds to conditions where population densities 
of the energy levels deviate from the equal ibrium Boltzmann 
distributions based on a single temperature. At sufficiently low 
densities, the probability of a radiative transition becomes 
comparable with the probability of a correspond i ng collisional 
transition. Unless the gas is optically thick, the emission of a 
photon is not balanced by its inverse. Consequently, the 


8-5 


population distribution among the energy levels departs from the 
predicted by the Boltzmann equation (25,28). In the formulation 
of the radiative transfer equation, Tiwari, etc. (25,29) devised a 
model that the absorption coefficient of different species is 
calculated by using the 58-step absorption model developed by 
Sutton (22), while the major nonequ i 1 i br i urn effect will enter 
through the source function in the equation. For a multilevel 
energy transitions model, this is done by introducing a parameter 
n which represents the ratio of the collisional relaxation time 
and the radiative lifetime of the first excited state, and the 
higher level energy transitions can be related to this parameter. 
The nonequilibrium radiation become important for conditions when 
n = 0 (1). The key to this model is the determination of 
collisional relaxation time between various particles and 
radiative lifetime of the excited states. Additional comparison 
and analysis are needed for this model. 

Other Radiation Models There are a number of other radiation 
models in the literatures (30-32). Further investigations and 
testings are needed for these models. In the field of 
nonequilibrium hypersonic flow with thermal radiation, it is 
desirable to have more • theoret ical analysis, computational code 
developments, and experimental data collections. 

Integration of Transfer Equation 

.To obtain the radiant heat flux vector in the energy 
equation, the radiative transfer equation needed to be integrated 
spectrally as well as spatially. This is very difficult to 
perform on a flow field over a whole vehicle. For an 
approximation, the tangent slab approximation technique is 
commonly used (20,22,23,29). This approximation implies that the 
radiative energy transfer along the body is negligible in 
comparison to that transferred in the direction normal to the 
body. Other approximation such as optical thin or/and optical 
thick assumptions may also be used to solve the problem (18). In 
general, additional works are needed in the area of fast spatial 
integration of transfer equation. 

CONCLUSIONS 


Several radiation models were found in the literature servey. 

These models need to be tested by some existing 3-D Nov i er -s to kes 
computer codes. It is proposed that the VRFL code is used for 
this purpose (7). From the result of testing, comparisons can be 
made on accuracy of prediction, requirement of computer time and 
memory space, and generality of the model. To save computer time 
and to facilitate easy implementation, these models should be 
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simplified but retained their essential features. In the computer 
tests, the flow field may be limited to the stagnation region of a 
symmetric blunt body for tangent slab approx imat ion, use minimum 
number of temperatures for thermal model, and assume radiation 
heat transfer equation and gas dynamic equations are uncoupled. 
From this testing and additional analysis, a modified and 
optimized radiation model may be obtained. 

After the optimized radiation model is incorporated into the 
VRFL code, the heat transfer rate to a entire hypersonic vehicle 
should be calculated. The radiation heat transfer equation and 
gas dynamic equations are still uncoupled but the tangent slab 
approx imat ion is removed. The purpose of this test is to find a 
numerical model for the spatial integration of the radiative 
tr ansfer equation. 

Finally, the expanded VRFL code can he applied to vehicles 
such as AOTVs. For these computations, the coupling of radiation 
and gas dynamics should be maintained. The purpose of calculation 
is not only to accurately predict the heat transfer to the 
vehicle, tout also to find the effect of radiation to the 
hypersonic flow field. 

New radiation model and numerical solution technique may 
still need to be developed. At present there are still too many 
unanswered questions in the field of nonequilibrium hypersonic 
flow with thermal radiation. 
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ABSTRACT 


f 


It is well established that plasma volume decreases during spaceflight 
and simulated weightlessness (bedrest). The decrement in plasma volume is 
thought to contribute to the orthostatic intolerance that has been 
observed in some crew members following spaceflight. To date, no studies 
have evaluated the effectiveness of fluid countermeasures of varying 
osmolality in the restoration of plasma volume and orthostatic tolerance 
in a controlled study. The overall objectives of this project were to: 1) 
provide a model that would rapidly and safely produce a fluid loss 
comparable to that which occurs during short duration spaceflight and 2) 
design a study that would determine the optimal drink solution to restore 
orthostatic tolerance and describe the mechanism(s) whereby orthostatic 
tolerance is restored. In the first part of the study, we used a diuretic 
(Lasix) as a model for the plasma volume changes that occur during short 
duration spaceflight. Four subjects (3 males and 1 female) participated 
in the study. Blood samples were drawn before Lasix administration (IV) 
and every 30 minutes for 3 hours after Lasix administration. Changes in 
plasma volume and plasma osmolarity were determined for each blood sample. 
Urine flow was followed for 3 hours after Lasix injection. Plasma volume 
decreased by an average of 11.5 £ .78% by hour 2 and then leveled off. 

Most of the loss of plasma occurred in the first 30 minutes after Lasix 
administration (-9.5 £ .45%). Plasma volume was lowest at 2 hours and 
increased slightly between 2 and 3 hours (11.5 £ .85% to 10.35 £ .95%). 
This might be a result of fluid shifting from other fluid compartments 
(interstitial and intracellular) into the plasma. There was a slight 
decrease in plasma osmolality during the three hours following lasix 
(286.3 to 283.6 £.85 mosra), thus the fluid loss was primarily isotonic. 
Urine flows peaked at approximately 40-60 minutes after lasix (14.7 £ 2.4 
ml/min) and returned to pre-lasix values by 3 hours. These data suggest 
that Lasix can be used as a model for fluid loss during short duration 
spaceflight where other models (i.e., bedrest) are not appropriate because 
of the experimental design. In the second part of the study, we have 
designed an experiment to examine the effectiveness of fluid 
countermeasures in the restoration of plasma volume and orthostatic 
tolerance. For each subject, a euhydrated control orthostatic challenge 
(lower body negative pressure - LBNP) will be performed in order to 
establish a baseline response. A second condition will involve 
dehydration by Lasix administration (20 mg IV). Two subsequent 
dehydration experiments will be performed where LBNP will be administered 
following rehydration with a hypertonic saline solution and an isotonic 
saline solution. This will be done to produce a range of plasma volumes 
and osmolarities to compare LBNP responses. Plasma volume, leg volume, 
plasma osmolarity, forearm blood flow, heart rate, blood pressure and 
cardiac dimensions (echocardiography) will be measured during rest, LBNP 
exposure, and recovery. During spaceflight, man is subjected to 
gravitational extremes (weightlessness to high G+ force) . It will be 
important to document the effectiveness of countermeasures that could 
potentially minimize the negative effects of spaceflight. 
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INTRODUCTION 


The cardiovascular system adapts to the weightless environment of 
space. This is sometimes referred to as cardiovascular deconditioning. 
The extent of the deconditioning and the mechanism of deconditioning vary 
with the individual and the duration of spaceflight. The effects of the 
deconditioning are manifested upon return to the earth's atmosphere. 
Crewmembers have experienced orthostatic intolerance when they are re- 
exposed to the +1 G environment. Symptoms include tachycardia, a 
reduction in systolic blood pressure and pulse pressure, and a reduction 
in exercise capacity. In order to maintain crew safety, it is necessary 
to understand the mechanisms that underlie this cardiovascular 
deconditioning and find potential countermeasures to minimize its 
negative effects. 

During spaceflight, there is a redistribution of fluid from the lower 
extremities to the head and upper body. It is thought that this central 
blood volume expansion results in a suppression of antidiuretic hormone 
(ADH) and a potentiated diuresis. This has been reported to decrease 
plasma volume by 8-15% (13). Current techniques to simulate the fluid 
shifts during spaceflight are limited to bedrest and water immersion. 

Both of these are time consuming and expensive and not appropriate for 
those studies that involve repetitive testing. 

As a result of the plasma volume loss, the central hypervolemia that 
occurs initially is compensated for and thus an individual is well suited 
to a weightless environment. However, upon reentry crewmembers are 
exposed to 1.2 to 2 +Gz, and at this point due to the hypovolemia, some 
individuals experience cardiovascular instability. This is characterized 
by decreased blood pressure, decreased brain blood flow and in some cases 
syncope. The potential for these to occur continues and is further 
manifested upon egress from the spacecraft. 

Countermeasures are currently being developed to minimize the 
cardiovascular deconditioning that occurs during spaceflight. These 
include exercise, lower body negative pressure (LBNP) during spaceflight, 
and fluid ingestion. In short duration flight it is thought that the 
primary factor involved in deconditioning is the loss of plasma volume. 

In longer duration flights, there could be venous compliance or 
baroreceptor sensitivity changes that contribute to the cardiovascular 
deconditioning . 

Fluid countermeasures have been used for the past 15 years in both the 
Russian and US space programs. The initial countermeasures included water 
and saline in the form of bouillon to offset the natriuresis and diuresis 
that were thought to occur in response to the headward fluid shift. 

Since that time the makeup of the countermeasure has changed considerably 
including a wide variety of electrolyte supplements. The current US 
countermeasure is salt tablets and water to make an isotonic solution 
(0.9% saline) taken 2 hours before reentry. One recent study by Bungo 
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and Charles (2) has indicated that crewmembers who ingested a salt and 
water solution before reentry had fewer incidences of orthostatic 
intolerance upon re-exposure to gravity than those that did not use the 
countermeasure. This has yet to be validated in well controlled ground 
based research. In addition, we dont know if the mechanism for the 
restoration of orthostatic tolerance is the restoration of plasma volume. 
It is possible that there is a restoration of hydration status in other 
body fluid compartments and that this is restoring orthostatic tolerance. 
In order to evaluate the use of different drink solutions, it is first 
necessary to first develop a model for spaceflight fluid loss that can be 
used in a repeated treatment design. 

Purpose of Research 

The overall objectives of this study were: 1) develop a ground based 
simulation of changes in plasma volume during short duration spaceflight 
and 2) develop a research protocol that would determine the optimal drink 
solution to restore orthostatic responses and describe the mechanism(s) 
whereby orthostatic tolerance is restored. 


HISTORICAL BACKGROUND 

The fluid shift that occurs during spaceflight can be as much as 2 
liters. This shift is thought to increase central venous pressure at least 
transiently, although this remains quite controversial (15). Atrial 
receptors are sensitive to changes in atrial pressure. Increasing stretch 
on these receptors decreases output to the anterior pituitary and thus 
decreases the output of ADH (6). This in turn decreases the reabsorption 
of water in the distal tubule of the kidney. In the non-human primate, 
this mechanism has been questioned (8). There was no suppression of ADH in 
response to volume expansion in the monkey. A study on Sky lab failed to 
observe a diuresis. They report that urine flow was actually suppressed 
during spaceflight (16). 

Recent evidence suggests that the increase in CVP is transient and might 
not be sufficient to stimulate a suppresion of ADH (21). However, it might 
be sufficient to stimulate an increase in the release of ANF. There are 
several resports of increased levels of ANF in response to central volume 
expansion by head-out water immersion (4, 22). There have also has been 
reports of a decrease in fluid intake during flight, which could also 
contribute to a reduction in plasma volume. Whatever the mechanism, there 
is a decline in plasma volume during spaceflight and simulated 
weightlessness. The decrease in plasma volume has been reported to be 
between 8-15% depending on the length of the flight (13). During simulated 
weightlessness (bedrest), the blood volume loss is equivalent to what is 
observed in spaceflight however, the plasma volume losses after horizontal 
bedrest are greater (13). Greenleaf et al. (10) reported an average plasma 
volume decrease of 15.2% following 14 days of bedrest. 
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Other studies have demonstrated that dehydration either as a result of 
heat stress (18), diuretic administration (19), or long duration exercise 
(11) decreases LBNP tolerance. Luft et al. (18) found that leg volume 
increased and arm volume decreased progressively with LBNP, but these 
changes were significantly less after dehydration. In a separate study, 
Luft et al. (19) found that LBNP tolerance was reduced following oral Lasix 
administration. The Lasix produced a 16.8% reduction in plasma volume. 
Hilton et al. (11) found that after dehydration as a result of 60 minutes 
of moderate exercise, LBNP intolerance occurred even after subjects were 
rehydrated . 

To date, there have been no studies that have identified the source of 
the fluid loss. Total body water has been shown to decrease by 3% in short 
duration shuttle flights (17). It has been previously mentioned that plasma 
volume decreases by approximately 10-15% during spaceflight. It is 
unlikely that all of the fluid was coming solely from the plasma. There 
must be equilibration of body fluid compartments in the latter portions of 
the flight. It has yet to be determined how much fluid is coming from the 
interstitial space (ISF) and/or intracellular space (ICF). 

As a result of cardiovascular deconditioning that occurs during 
spaceflight, there is an impairment of orthostatic responses (7,14). 

Several studies have indicated that orthostatic tolerance following bedrest 
is also impaired (10,12). In an attempt to reduce some of the potentially 
debilitating cardiovascular deconditioning of spaceflight and its effect on 
orthostatic tolerance, three potential countermeasures have been developed: 
exercise, LBNP during spaceflight, and fluid ingestion before reentry. The 
effects of exercise on cardiovascular deconditioning have yet to be 
documented. However, anecdotal information collected on the Sky lab missions 
indicates that recovery from spaceflight was faster in those who 
participated in an inflight exercise program (20). Inflight LBNP has been 
used during Sky lab with mixed results (14). However, a ground based 
bedrest study where LBNP was applied daily has demonstrated that LBNP can 
be used to offset some of the effects of weightlessness (12). 

Several countermeasures have been used to offset the naturesis and 
diuresis that was observed during spaceflight. One of the first ingested 
countermeasures was 9-alpha-fluorohydrocortisone which minimized bed rest 
orthostatic intolerance (1). Hyatt and West (12) used one of the first 
fluid countermeasures and found that the combined use of oral saline and 
LBNP was quite effective in returning the heart rate and blood pressure 
responses to LBNP to prebedrest values, but saline ingestion alone was less 
effective. Greenleaf et al.(10) reported that saline consumption improves 
acceleration tolerance after bedrest by restoring plasma volume. Recently, 
Bungo et al. (2) examined the use of fluid countermeasures in several 
shuttle flights. They found that crew members that ingested a salt and 
water solution before reentry had fewer incidences of orthostatic 
intolerance upon re-exposure to gravity than those that did not use the 
countermeasure. In a recent study by Fry et al. (5), different drink 
solutions were ingested in order to determine the solution that produced 
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the greatest elevation in plasma volume and maintained the increase the 
longest. They found that a slightly hypertonic saline solution (1.07%) was 
most effective in increasing and maintaining the higher plasma volume. It 
was superior to a number of drink solutions including an isotonic saline 
solution (the present countermeasure). It has yet to be determined whether 
the hypertonic solution will increase plasma volume after dehydration and 
whether it will improve orthostatic tolerance as well as the isotonic 
saline fluid countermeasure. 

There are various techniques to assess orthostatic responses. These 
include the stand test, tilt table test, human centrifugation, and lower 
body negative pressure (LBNP). The NASA stand test has the subject lie 
supine for 5 minutes and then stand quickly. Heart rate and blood pressure 
are then monitored for 5 minutes. This is a simple test, however, it does 
not rule out the confounding effects of muscle contraction and the 
resulting volume redistribution associated with it. By putting an 
individual on a tilt table and positioning the table to various inclines, 
the gravitational effects on the cardiovasular system can be varied. 
However, this test has limited usefulness in that it can only place a +1 Gz 
stress on the individual and it is difficult to get many physiological 
measurements during this type of testing. Human centrifugation is 
expensive and data collection is limited by the size of the centrifuge. 

Lower body negative pressure allows for the gradation of the orthostatic 
stress while the subject is in the supine position. This technique 
facilitates data collection and allows for high +Gz forces to be applied to 
the subject. Lower body negative pressure results in venous pooling in the 
lower extremities. This in turn decreases cardiac filling pressure, end 
diastolic volume, and ultimately via the Frank-Starling mechanism stroke 
volume. Reflex mechanisms from both the carotid baroreceptors and 
cardiopulmonary receptors are activated which act to increase heart rate 
and total peripheral resistance to maintain arterial pressure. In 
addition, there is an increased sympathetic outflow and catecholamine 
release resulting in a reduction in end-systolic volume, thus maintaining 
stroke volume and cardiac output (24). 

To date, no studies have evaluated the effectiveness of fluid 
countermeasures of varying osmolality in the restoration of plasma volume 
and LBNP tolerance in a controlled study. Studies of this nature have been 
limited by the need to do repetitive spaceflight simulations. For bedrest, 
there is an effect of one bedrest on the next bedrest making it impossible 
to look at independent effects of drink solutions in subsequent bedrests. 
Diuretics will cause a dehydration as a result of increasing fluid loss in 
the urine. Lasix is a commonly used diuretic for the treatment of 
hypertension that uses furosimide as its active ingredient. Furosimide 
acts to inhibit the reabsorption of sodium and chloride in the proximal and 
distal tubules, and the loop of henle. Its action is independent of any 
inhibitory effects of carbonic anhydrase or aldosterone (9). It is quick 
acting and has few side effects. For these reasons, we tried Lasix as a 
model for the fluid losses that occur during spaceflight. 
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METHODS 


Model for Plasma Volume Changes During Spaceflight 

In order to determine the total loss of plasma volume as a result of 
diuretic administration (Lasix), the time course of PV changes, the 
optimal time for fluid ingestion, and whether the loss of fluid with the 
diuretic is isotonic, four subjects were recruited and participated after 
obtaining informed consent (3 males and 1 female). Subjects entered the 
laboratory after a small meal in the euhydrated state. Subjects voided 
their bladder and were weighed. They then sat for 30 minutes to control 
for posture before blood sampling. An intracath was inserted in a 
superficial arm vein. A preliminary blood sample (5 ml) was drawn without 
stasis. Twenty mg of Lasix was injected IV. Urine samples were collected 
for 3 hours after Lasix injection for urine flow determinations. Blood 
samples were obtained every 30 minutes after Lasix. Hematocrit was 
determined by microcentrifugation, hemoglobin was determined by a Coulter 
Counter, osmolality was determined by freezing point depression on all 
samples. From the hematocrit (HCT) and hemoglobin (HG), relative changes in 
plasma volume were determined using the formula of Dill and Costill (3). 

Fluid Countermeasure and Orthostatic Tolerance 

Eight volunteers will be recruited for this part of the study. Subjects 
will be given several practice LBNP tests prior to the start of the 
experimental phase of the study. For each subject, control LBNP test(s) 
will be performed in order to establish a baseline response to LBNP. In a 
second condition, dehydration will be produced by a 20 mg IV injection of 
Lasix. This will produce approximately an 11% reduction in plasma volume 
from euhydrated plasma volume as has been determined from the first part of 
this project. Two subsequent experiments will be performed where each 
subject is rehydrated after dehydration by Lasix with a 14 ml/kg solution 
of a hypertonic saline solution and an isotonic saline solution. For the 
dehydration and two rehydration experiments an LBNP test will be performed. 
This will give us a range of plasma volumes and osmolarities to compare 
LBNP reponses. Each subject will complete all tests. Below is a summary of 
the 4 conditions: 

1. Control LBNP 

2. LBNP following dehydration 

3. LBNP following dehydration and rehydration with hypertonic saline 

4. LBNP following dehydration and rehydration with isotonic saline 

The order of treatments will be randomized in order to eliminate any 
order effects on the LBNP results. All experiments will begin at the same 
time of day for each individual subject to eliminate any circadian effects 
on the findings. Testing will begin 2 hours after the ingestion of a small 
meal. LBNP will be performed 60 minutes after ingestion of fluids. At 
least 72 hours will be allowed between testing. 
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LBNP Protocol 


Lower body negative pressure will be applied using a plywood chamber 
sealed at the waist with a rubber gasket. A vacuum cleaner is then used to 
withdraw air out of the box. The following protocol will be used: 

Pressure (Torr) Time at Stage(min) 


Rest 

0 

10 

Stage 1 

-5 

5 

Stage 2 

-10 

5 

Stage 3 

-20 

5 

Stage 4 

-30 

5 

Stage 5 

-40 

5 

Stage 6 

-50 

5 

Stage 7 

-60 

5 

Recovery 

0 

10 


The test will be terminated if the subject’s systolic blood pressure 
drops suddenly (greater than 25 mmHg in 1 minute), systolic pressure 
reaches 80 mmHg, bradycardia occurs (drop in HR greater than 15 BPM), 
subject distress or subject request. 

Before, during, and immediately after LBNP the parameters described below 
will be measured. 

Leg volume is calculated from leg circumference measurements made every 
minute with a mercury-in-silastic strain gauge. The changes in leg volume 
during LBNP exposure will be used to approximate the amount of venous 
pooling in the legs. 

A blood sample will be taken from an antecubital vein before drink 
ingestion, 30 minutes after drink ingestion, immediately before, and 
immediately after LBNP (a total of seven blood samples). Changes in plasma 
volume will be calculated from hematocrit (microhematocrit technique) and 
hemoglobin (Coulter Counter) ratios using the formula of Dill and Costill 
(3). A direct measurement of plasma volume will be performed once (human 
serum labeled 125-1) in each subject before the first experimental testing. 
This value will be used as the absolute plasma volume and changes will be 
expressed relative to that value. Plasma osmolality will be determined 
using freezing point depression on all blood samples. A 5 ml sample will be 
required to do all of the hematological analysis. Plasma levels of 
antidiuretic hormone (ADH), atrial natriuretic factor (ANF), aldosterone, 
and catecholamines (epinephrine, norepinephrine) will be measured for each 
blood sample. These assays require a blood sample of 15 ml. 

Venous occlusion plethysmography with a mercury-in-silastic strain gauge 
will be used to measure forearm blood flow. The forearm vasoconstrictor 
response as indicated by this technique will be measured at each step of 
LBNP, and plotted as a function of systolic blood pressure. The slope of 
this response will be used to estimate vasoconstrictor responses to 
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baroreceptor (cardiopulmonary and sinoaortic) activation induced by the 
decreasing blood pressure. Heart rate will be determined using a three 
lead EKG. The R-R interval will be determined for each stage of LBNP. The 
change in R-R interval will be plotted as a function of the systolic blood 
pressure. The slope of this relationship will be used as an estimate of 
overall baroreceptor function. Blood pressure will be measured once a 
minute using a standard auscultatory technique. 2-D and M-Mode 
echocardiography (ATL 4000 S/LC ultrasound system, ATL, Botheli, WA.) will 
be performed at each stage of LBNP in order to assess relative changes in 
left ventricular dimensions with LBNP. 

EMG will be used to assess abdominal and leg muscle tension in order to 
control for muscle tensing. The subject will be asked to maintain a resting 
EMG level throughout all LBNP tests. 

Urine will be collected 30 minutes after drinking and immediately after 
LBNP in order to calculate urine flow rates. 

RESULTS 

Model for Changes in Plasma Volume During Short Duration Spaceflight 

Plasma volume decreased in the first 30 minutes after Lasix 
administration (Figure 1). The loss of plasma volume leveled off at 2 
hours (11 ± .78%). By 3 hours, the plasma volume was partially restored, 
most. likely as a result of fluid shifting from the ISF and ICF compartments 
into the plasma. 

There were no major changes in plasma osmolality for three hours after 
Lasix in comparison to pre-lasix injection (Figure 2). There was a slight 
fall in plasma osmolality (286.3 to 283.6 + .85 mosm), indicating a 
slightly hyperosmotic urine. 

Urine flow for each subject is displayed in Figure 3. For all of the 
subjects, urine flow peaked between 30 and 60 minutes, then decreased back 
to control level. By 3 hours, urine flow had returned to the control value. 

In an attempt to determine the source of the fluid loss, the % loss of 
fluid from the plasma was estimated. We assumed that the plasma volume was 
45 ml/kg for men and 35 ml/kg for women. These results are displayed in 
Table 1. Overall, for the 3 male subjects, one third of the fluid loss was 
from the plasma, and thus the other two thirds must have been from the ISF 
or ICF. The level of dehydration (-11% of the plasma) and the time (2-3 
hours) would indicate that most of the fluid was coming from the ISF and 
not the ICF. 

DISCUSSION 

Six degree head down bedrest is currently the standard simulation of 
weightlessness. However, it is expensive and can not be used in repetitive 
designs. Results from the present study indicate that Lasix produces 


9-9 


comparable losses of plasma volume (11%) to that observed following 
spaceflight (10-15%) (13). This is also comparable to oral Lasix 
administration. Luft et al. (19) reported plasma volume decreased 16.8% by 
3 hours after oral Lasix. However, they did not follow the time course of 
the response. Lasix is fast acting as most of the loss of plasma volume 
occured in the present study in the first 45 minutes following IV 
injection. By 2 hours the loss of plasma had peaked. Thus for the second 
part of this study, fluid ingestion should take place approximately 2 hours 
after lasix administration. As is the case during spaceflight the loss of 
fluid as a result of the lasix is primarily isotonic. Although, it is not 
known during spaceflight where the fluid loss is coming from, it is 
probably coming from the ISF first and then the ICF. In the present study, 
we found that approximately 1/3 of the total fluid loss was coming from the 
plasma and the rest from the other two fluid compartments (ICF and ISF). 

There is the possibility that Lasix could influence the control of blood 
pressure independent of the plasma volume changes. Lasix does alter renal 
blood flow. It can increase or decrease renal blood flow depending on the 
experimental conditions (9). Furosimide has also been shown to increase 
atrial natriuretic factor and renin in plasma (23). However, it is 
difficult to separate out the hypovolemic effects of furosimide and its 
direct effect on the kidney and other vascular areas. This needs to be 
taken into account when using Lasix under these conditions. Given these 
limitations, Lasix can be used to simulate the plasma volume losses that 
occur during short duration spaceflight. 

The design presented should determine the optimal drink solution to 
restore orthostatic tolerance and the mechanism whereby orthostatic 
tolerance is restored. We would anticipate, based on the study by Fry et 
al.(5), that the optimal drink solution should be the hypertonic saline 
solution if the mechanism for restoring orthostatic tolerance is the 
restoration of plasma volume. However, it is possible that other 
mechanisms might be operational. It is possible that fluid ingested moves 
through the plasma into the ISF and ultimately the ICF, thus restoring 
hydration in cells. This could in some way restore orthostatic tolerance 
independent of restoring plasma volume. 

CONCLUSIONS 

In summary, Lasix can be used as a way of simulating the plasma volume 
changes that occur during short duration spaceflight. The total loss of 
plasma is comparable to spaceflight. Lasix is fast acting, and has 
relatively few side effects. 

The present design for evaluating the optimal fluid countermeasures will 
have important implications in restoring orthostatic tolerance and function 
in the latter stages of spaceflight when it is essential for safe operation 
of the spacecraft. 
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Figure 1. - Change in plasma volume as a function of time 
after Lasix administration. Means are plotted + 1 SE. 
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Figure 2. - Plasma osmolality as a function of time after 
Lasix administration. Means are plotted ± 1 SE. 
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Figure 3. - Urine flow rate as a function of time after Lasix 
administration. Individual responses are displayed. 
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Table 1. - Total urine output and % loss of fluid from 
plasma for each subject. 


SUBJECT 

TOTAL URINE OUTPUT 
(LITERS) 

% FROM PLASMA 

1 

1.58 

34.5 

2 

1.13 

35.8 

3 

1.40 

31.6 

4 

0.72 

51.7 


CALCULATIONS: 

EST.PV (ML) X % LOSS OF PV 
% LOSS FROM PLASMA = 

TOTAL URINE OUTPUT (L) 
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ABSTRACT 


Some recent problems are no longer formulated in terms of imprecise facts, missing data or 
inadequate measuring devices. Instead, questions pertaining to knowledge and information itself 
arise and can be phrased independently of any particular area of knowledge. The problem considered 
in the present work is how to model a problem solver that is trying to find the answer to some query. 
The problem solver has access to a large number of knowledge systems that specialize in diverse 
features. In this context, feature means an indicator of what the possibilities for the answer are. The 
knowledge systems should not be accessed more than once, in order to have truly independent sources 
of information. Moreover, these systems are allowed to run in parallel. Since access might be 
expensive, it is necessary to construct a management policy for accessing these knowledge systems. 
To help in the access policy, some control knowledge systems are available. Control knowledge 
systems have knowledge about the performance parameters status of the knowledge systems. In 
order to carry out the double goal of estimating what units to access and to answer the given query, 
diverse pieces of evidence must be fused. We use the Dempster-Shafer Theory of Evidence to pool the 
knowledge bases. The present work demonstrates how pooling evidence can determine the necessary 
strategy to solve the problem as well as solving it. The present work deals with the path the problem 
solver could follow to answer the query. 
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INTRODUCTION 


Some recent problems are no longer formulated in terms of imprecise facts, missing data or 
inadequate measuring devices. Many problems today deal with knowledge and information itself and 
can be phrased independently of any particular area of knowledge . Problems in such diverse fields as 
psychology, engineering, artificial intelligence and decision making can exhibit similar formats when 
it comes to dealing with information and uncertainty. It is desirable to be able to handle information, 
not in any particular context, but in the form of general principles valid in many situations. 

In [1] [2] and [8], the authors develop a model for a course of action-feedback loop that is useful 
when considering learning in repeated environments. This model is not applicable to the case where 
one-shot decisions are to be made. For such decisions, the decision maker accumulates information 
and after enough accumulation has taken place, certain alternati ves are ruled out. This process is 
repeated till there is only one remaining alternative. Presumably this last alternative represents the 
decision to be taken. This model is studied in [9]. In this setting, it is essential to be able to update 
the general state of information every time an alternative is ruled out. 

A powerful method in artificial intelligence is to look at certain features of a problem and try to 
put together the evidence coming from these features in order to recognize the corresponding pattern. 
In this context, feature means any indication of what the possibilities of a correct answer are. Special 
cases of such an approach are, for example, taken in computer vision, see [3], [13], and [14], It is clear 
that with such an approach we need to have some kind of mechanism available for evidence pooling. 
We shall use the D-S theory [12] to carry out this pooling of evidence. In fact, the D-S theory will be 
used for two purposes: (a) plan the solution to a query by accessing some pertinent knowledge 
systems (b) define dynamically a policy that would determine what KSs to access at a particular 
time. The problem we consider here is quite general. Given a query Q we look at a set of possible 
answers to Q. We of course would like to determine the precise element of the set that represents the 
correct answer to Q. In order to accomplish this, specialized knowledge systems are consulted. The 
execution is driven by a table that to each goal associates a sequence of control strategies. A control 
strategy is a set of simultaneous constraints on the desired performance characteristics of the 
knowledge systems to be consulted. A number of control strategies weighted by the strength of belief 
in each constraint is given in the table. The weight of each constraint is the degree of belief that a 
specific value, for some performance parameter, is optimal. This will be compounded with the 
uncertainty that some family of knowledge systems has the desired performance parameter value. 
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After the completion of each goal, new data will be read in. A control sequence defines a 
subcycle of the goal execution. After the completion of each subcycle, an update on the information 
status is performed. In order to obtain information about the knowledge systems, special knowledge 
systems, called control knowledge systems, are consulted. Based on the pooled evidence provided 
dynamically by the control knowledge systems, a policy of access is determined. 

There are two common ways of representing uncertain information about propositions. The 
most straightforward (employed in PROSPECTOR) provides a number P(AIE) which is the 
probability of proposition A being true, given current evidence E. One problem with this approach is 
that the precision of P(AIE) is not known. In fact, Garvey et al note "A likelihood represented by a 
point probability value is usually an overstatement of what is actually known.” Another problem 
with this approach is that we do not have the amount of evidence for and against A. The other way to 
proceed is to attach not one but two numbers with each proposition (MYCIN employs this method). 
The two numbers correspond to measures of belief and disbelief. This approach certainly removes one 
of the above criticisms. It is still true that the precision relative to the two numbers is not known. 
Another criticism that may be leveled is that no formal theory indicates how evidence for and against 
might be combined. Also, it is not clear how to detect conflict. Finally, it has been shown that results 
degrade quite fast as uncertainty increases. The D-S theory is used here because it avoids, for the 
most part, the problems outlined above. 

We assume that a large number of knowledge sources is available because we do want to access 
independent sources and therefore do not want to access the same source more than once. In this 
context, "independent sources” means that the errors in these sources are independent. To 
approximate independence, we may force a delay before any knowledge source is reaccessed. When 
all the steps have been taken, a final evaluation is made to determine the answer to the query. 



CONCEPTS AND NOTATIONS 


We start by defining the concept of belief structure. For a complete account of general belief 
structures and related material, the reader is referred to [12], If Xdenotes some set and t(X) denotes 
the power set of X, a belief structure is a function m from t(X) into [0,1 ] satisfying 

( i) m(0) =0 
(ii) V m (A)= 1 

AC X 

Elements A € i(X) for which m(A)> 0 are called focal elements of rrt. 

Now we define the Belief and Plausability measures on t(X) 


If B € v(X) 


Bel(B) = X rn(A) , Pts(B) = £ m(A). 

AC B An«%0 


Suppose we have n experts whose opinions are equally respected. Assume that the opinion of 
the i th expert is that a specified object is to be found in some set A t . Given these n opinions we now ask 
the question: what is the probability that this specified object is in some set B? Assume for ease of 
notation that A;, A$, ..Ak are subsets of B, Ak+i, Ak+ 2 ,-.Ak+( intersect B and its complement and 
that Aic+(+i,...A n , are subsets of the complement of B. A lower bound on the probability that the 
object is in B is Kin, an upper bound is (K+t)/n. We hence see that a possible interpretation of 
BEL (B) and Pts (B) is a "Lower probability and a higher probability” of B. Another way to interpret 
BEL and Pts is to identify the set B with the proposition: "The object is in the set B.” The truth or 
falsity of this proposition is determined by the available evidence, in this case, the opinion of the n 
experts. With this interpretation Bel (B) is interpreted to be the degree of support for proposition B. 
Pts (B) is the degree to which one fails to refute proposition B (Aj, A 2 ,... Ak+i,...Ak+( all intersect B 
and this constitute evidence that fails to refute B). With this interpretation, Pts (~>B) is interpreted 
as the degree to which proposition B is refuted while Pts(B) - Bel(B) is the degree of ignorance about 
proposition B, indeed: 


Pts(B) - Bel(B) = 2 K+t m (A.) 

i=K + 1 ‘ 
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and since the above A* do intersect B and its complement, they fail to refute B and they fail to refute 
~>B. Ignorance on proposition B is the degree to which one fails to refute B and ~>B. An important 
property of belief structures is that they allow to fuse evidence from independent sources of 
information. If mj and m 2 are belief structures defined on the same set X then the combination rule, 
see [12], states that they combine to yield a belief structure m 3 on X. The focal element of m 3 are 
given by taking all possible intersections of focal elements of mi with focal elements of m 2 and 

m 3 (C) = ^ m l (A)m 2 (B)/(l~K) 

A ns = C 

K is the measure of conflict between the two sources and comes from pairs of disjoint focal 
elements. In fact 


K = £ m l (A)m 2 (B) 

A ns = 0 

If m 3 is given by the formula above, we write m 3 = mi® m 2 . The formula above is true if the 
information sources are independent. For a very readable interpretation of the combination formula, 
the reader is referred to the article by L.A. Zadeh [16]. 

In [1], [2], and [8] the authors have defined the concept of an abstract information system. In 
this model, decisions are made sequentially and information is received as a consequence of the 
cource of action taken. That is, each cycle of the decision-information acquisition loop is initiated by 
the decision maker (DM) choosing a cource of action ( COA) and the data generated by this CO A is 
fedback to the DM which completes the cycle. This model is most relevant to situations requiring 
frequent decisions on a regular basis such as ordering equipment or choosing stocks. The DM can 
take the same COA on successive cycles or change to different COA, possibly returning to a COA 
chosen earlier. Uncertainty has been incorporated into this model. Different policies for selecting a 
COA on a given cycle have been discussed in [7]. The model discussed in the works cited above is not 
relevant to "one-shot decisions” where the DM uses all available information to pick a single COA 
which, once chosen, is irreversible. 

In [9] the authors study situations where one-shot decisions are called for. In this context there 
may be many minor decisions, such as which information source to access and how to interpret 
information received. The DM constructs a set of possible CO As and then, based on the information 
received, eliminates the less promising CO As until only one remains. This, of course, is the COA 
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taken. The main idea in [9] is that for each fixed data vector, the DM creates a conditional belief 
structure and then after K cycles, combines these masses conditioned by the data vectors into a single 
mass conditioned by history. From that, the belief and plausability measures conditioned on history 
are formed and whenever the plausability of a COA* falls below the belief of some other COAj, COA l 
is eliminated from consideration. The general problem of gathering evidence through classical 
statistics is that, as already pointed out in the introduction, an enormous amount of data is required. 
In addition, with such approaches, results quickly deteriorate if uncertainty is present. Usually some 
version of Baye’s rule of inference is used as a remedy to the larg e amount of micro-events present, see 
for example [4] and [5]. In [10] J. Pearl makes the point that a complete probability space is not 
required. One needs only to estimate likelihood ratios. However a large number of such likelihood 
ratios would still be required. 

In the next section we will outline an approach that uses the DS theory to collect evidence from 
different cycles as in [9], We will collect evidence by analyzing different features of the problem. 
Breaking a problem down to subsets of features has already been tried in diverse fields of artificial 
intelligence. For example for computer vision, see [3], [13], and [14], In most of the existing 
approaches, information on specific features is relatively static and no update of that information is 
incorporated into the DS approach. In this work, an update mechanism will be built into the decision 
making process. Also the DS theory will not only be used to identify a solution of the problem but also 
will determine how to access the knowledge sources (KSs ) Since access to some KSs might be difficult 
and since we will at the very least delay reaccessing a KS, it is particularly important to determine a 
viable policy for access. We will also allow different KSs to be m ulti-tasked on particular cycles or 
sub-cycles. Finally it is important to have a model that is potentially applicable to many situations 
and for this reason we have made the discussion as general as possible. 


RESULTS 


The setting is as follows: Let Q denote some specific query and let Hq be the corresponding 
frame of discernment, i.e. Hq is a set of possible answers to Q and one of the elements of Hq is the 
correct answer. Of course we do not know a priori which element is the correct answer. The data 
vector x is of the form 
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where s is large and where € Fi(l^i^s). F, is a set of (not necessarily numerical) values 
pertaining to the i th feature. For example for some i, F; could be different values of light intensity. 
These values could range in the set {low, average, fairly high, high}. For another i, F, could be color 
values such as {red, blue, green} etc... Associated with each set F; we have a knowledge system KSi 
that only reports on values in F; . Thus for some values of i, KSi may report only on the value of the 
color. Since we would like in principle, not to access any KSi more than once (in order to truly have 
independent sources of information), we assume that s is large enough to yield appropriate 
information to obtain an answer to Q. Given a feature value € F, we define a corresponding 
subset of Hq by setting 



g € H^H th feature ofg has value/". 1 


A ( ft * ) can be identified with the proposition: "The i tfl feature has value /}*• ”. In fact if G C Hq then 
G can be identified with the proposition: "The answer is in G”. Thus r (Hq) can be identified with the 
set of relevant propositions. 

As already mentioned earlier, the value of the i th feature will be reported by KSj. However, 
there is some uncertainty built into the reporting of KSi. Thus for the corresponding i, KSi might 
report that the intensity is low .6, is average .3 and is fairly high .1 

I 

i 

I 

In general, let a,* be the degree of belief that the i tfl feature value is f t k . Clearly 




We now define a belief structure m, : 

Another way to state this is to say that is the degree of belief that any element of A(f;*) is 
possibly the correct answer. Thus to sum up: in order to get the answer to Q, we look at the values of 
some features. These values are reported by specialized KSs with some degree of uncertainty built 
into the report. It is also clear that n*j, 1<> i ^ s constitute independent sources of information since 
they deal with distinct features. 
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Since there is a large number of KSs to access, it is important to formulate a policy on how to 
access KSis. We base this policy on performance characteristics. Let Pj, P2 , ... P( be sets of 
performance values e.g. Pi could be a set of values corresponding to the cost of access, P2 could 
correspond to response delays etc... To each Pi corresponds a special KS, distinct from the ones 
mentioned above, which knows about performance (as opposed to features). To distinguish these KSs, 
we call them control KSs and denote them by CKSs. Thus CKS t has knowledge of values in Pj. This 
knowledge is much more dynamic than the knowledge possessed by KSis, since performance of the 
KSs strongly depends on time. We would like to have some of the KSi run in parallel. To this end we 
set 

H p Ct { KSi, . . . KS a } 

to be all possible combinations of KSs that may run in parallel. Slome of the combinations are of 
course ruled out as some of the KSs might not be compatible or may have been already used. As 
earlier, for p;*‘ € Pj define 


a (v)«"p 

where each element of A ( pi ki ) is such that it’s ith performance characteristic value is p t ki . For 
example if p**' means the cheapest access and if A ( pfi) = {KS12 , KSig} then KS12 or KSig might be 
the cheapest KS to access. 


In order to find an answer to Q, we will have a sequence of goals to satisfy. The process will be 
table-driven in the sense that each intermediate goal will correspond to a sequence of control 
strategies. A control strategy will be a set of performance constraints to be simultaneously satisfied. 
Ultimate control strategies with their respective degrees of belie f may be given. For example a 
control strategy could be: get the cheapest access KSs (.8) and also the smallest response delays (.7). 
By a subcycle we mean the execution of a particular control strat egy. By a cycle we mean the 
execution of one goal. We use the symbol^ to refer to the k th subcycle of cycle j. 

On each subcycle, CKSi assigns a mass to A( pj*). Thus 
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where 6,* marks the degree of belief that for the i th performance it’s the performance value p t k that is 
needed and c t k is the degree of belief that any element in A ( pt k ) has i th - performance value p,*. 


We have ^ 6* = 1 for all t 

k 

0 5 c k S 1 for all i and k 

t 

and/r* ( KSyKS 2 , ... KS g )= 1— ^ 6* c* . It is assumed that we always have 

a(^)- {»,.«, ...o. 

That is, every p** determines a proper subset of {KSi, ... KSJ. Note that n{ k do satisfy the 
properties of a belief structure that the excess mass, after b/c;* has been assigned, is distributed over 
all KSs. Assume that a particular control strategy dictates that we look at particular values in 
performance spaces Pj,? 2 ,etc. ... Let Xj = {KSi r ... KSiJ be the set of KSs whose first performance 
characteristic matches pi k i, similarly let X2, X3, ... be sets of KSs whose performance characteristic 
matches P2 > ‘ ! ,P3 k 3 , etc... Of course this information is given by CKSi, CKS2, ... at the time the 
driving table is consulted for the appropriate control strategy to be followed. Thus the focal elements 
of up* are of the form r (Xi) X H p X H p X ... The focal elements of n^k are of the form 
H p X v(X2) X Hp X... etc... 


We now use the D-S theory to formulate a policy of selecting a set of KSs. Of course each X t 
reported by CKSi is assumed to be compatible with H p i.e., ifXi = {KS2 , KSg, KS10} and {KS2, KS5} 
$ H p but Xj 4 H p then KSjo will not be considered, but running KS2 and KSg in parallel will be 
considered. For example 




KS^KS C 


X H X H X 

p p 


-) 


will denote the degree of belief that for some value p/*i for access cost (assuming access cost is P;) 
that value p^i is the best to consider compounded with the degree of belief that KS2 and KSg have an 
access cost value of p^i . Now on subcycle jk if we consider, for example, three performance 
characteristics which, for ease of notation, we assume to be Pj, P2, P3 then we will consider the 
performance space 

= « (*,) X « (*,) X « (*,) 
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where € (X() (l^i^ 3 ) are the possible values of Xi (there may be several possible values o fXi due to 
uncertainty wherever p.^i is the best value to consider in Pi). If Yi, Y3 are possible values of Xj, 
X2, and X3 then nfo, for example, will be induced on the set {Yj} X € (X2) X € (X3) by 




X € 


X € 


x 4H‘( 


XH XH X 

p P 


thus, in contrast to Hq, H p jk (which is the set of answers of what KSs to use on subcycle j*) is highly 
dynamic and is generated by the appropriate control sequence. In general H p Jk is a set of elements of 
the form ( Aj , A2 , .. At k ) where t* is the number of performance characteristics considered on a 
particular control strategy on cycle j and each A is an appropriate set of KSs. We now combine the 
masses on H p j k by setting 


n = © n. 
i = 1 1 

The problem is to select the "best combination" V = ( Bi , i?2, ... B tk ). Note again that 
independence of information sources is justified as different CKSs are pooled. To get the "best 
combination ” U we maximize 

Bel(V) -Bel(-'V) 

over V € H p j k, and relative to the pooled belief structure n Jfc . The rationale for the optimization 
described above is to maximize the support of V over the support of its competitors. If several 
elements are tied for this optimization, we select one of the elements at random. If V is an optimal 
element we make an arbitrary selection 6; € Bi,(l^i^tk). Each 6; denotes some KS € B j. Now we 
follow the policy of running the in parallel. We of course keep track of which bi have been run and 
update H p appropriately after each subcycle, so no KS will be accessed twice. 

Now at the beginning of cycle j assume that the belief structure relative to subsets of Hq is 
given by mf. Of course mt itself was obtained by pooling together the belief structures m j where i 
ranges over all features considered prior to cycle j. On the first subcycle of cycle j, assume that the 
policy described above indicates that features ij, t2, ... it should be looked at. Pooling the 
corresponding KSs we obtain the belief structures m (Ji (l^v^t ) and we form 

Am 1 = © t . m. 1 

u = 1 1 
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At the end of cycle j, mf is updated by 


Arn'QAm 2 ... ©A m v 

where v denotes the number of control sequences in cycle j. 

At the completion of each subcycle we check if 


Pls(A(/J))<Be((A^)) forj* k 

We use the current value of raj to do this. Of course that current value is given by 

«— - 1 ©Am 1 ©... ©Am* 

if we have just completed subcycle y*. 

If the inequality above is satisfied we rule out the proposition " the value of the i th feature is 
fl k”. The rationale for doing this is if, for example, the upper probability for value green falls below 
the lower probability of some color value, say blue, then we should rule out that the color is green. Of 
course we could embellish this rule and require for instance that Bel (A ( fiJ) ), at least initially, 
exceeds Pis (A f ,*) ) by some margin which could be relaxed as data accummulates. We then should 
perform a post-elimination update. There are many ways to do this. A straightforward method would 
be to distribute the mass of A f /j*) uniformly over the remaining A(fj). Of course if we have some 
additional information, the distribution may be chosen not to be uniform. We also could use the 
plausability of the remaining values of F; to guide our redistribution of mass. 

At the end of the cycle we have reduced the possible number of values in some of the sets F;. We 
have at this point followed the actions dictated by the current goal. Now we take a new reading of the 
data and proceed to work in a similar manner on the next goal. A special note should be made when 
the control strategy is to reduce ignorance, which was defined in the previous section. Say the largest 
ignorance is for {qj , q$, qrf C Hq. The CKS in charge of ignorance (which is treated here as a 
particular feature) has access to all F; (1 5 1 s s) and the corresponding A (ft 1 ). To reduce the 
ignorance on {qi,qs,q7} one must get hold of those KSs that do not deal with qi, qs, qj simultaneously 
(otherwise the ignorance may continue to propagate). Thus if KSs has information about a set 
containing say qi but not qs or 47, then KSs is a viable candidate to reduce ignorance. 
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Finally when all the actions corresponding to all of the listed goals have been executed we end 
up with some of the sets Fj reduced (e.g. value, green and red may have been eliminated from Fs, if the 
third feature is color value) Now to get the answer to query Q we maximize 

Bel(q) -Bel(^q) ,qtH Q 


relative to the final value of m J . 

The rationale for the optimization defined above is of course the same as earlier. We want to 
maximize the support for a particular answer over its competitors. 

To sum up we have used the D-S theory to pool the evidence given by the CKSs to select a 
policy by which we would use the appropriate KSs to find information relative to query Q. Then we 
pooled the evidence collected from the KSs to rule out certain propositions about feature values. Each 
time we ruled out a proposition we had to redistribute its mass ov er remaining propositions 
pertaining to a specific feature. When all the steps are completed we pick the answer that maximizes 
support over its competitors. Of course, there are other ways to pool beliefs from different sources and 
for a sample of these methods we refer the reactor to [6], [1 17, [IS], and [1 7]. 

Finally, if it is desired to access KSs more than once, some delay must be built into reaccessing 
any KS in order to approximate the independence of information sources. Recently, theories more 
general than the D-S approach have been studied. More flexible rules than the D-S rule of 
combination have been developed. These settings will be investigated in future work. 


10-13 



REFERENCES 


1. Alo, R., Kleyle, R.M. and de Korvin, A. (1985), "Decision Making Mechanisms and Stopping 
Times in a Generalized Information System”, Mathematical Modelling. Vol. 6, pp. 259-271. 

2. Alo, R., de Korvin, A. and Kleyle, R. (1987), "Emergence of a Dominant Course of Action in a 
Generalized Feedback Loop when Goal Uncertainty is Present”, Journal of the American Society for 
Information Science. Vol. 38, pp. 111-117. 

3. Davis, L.S., Rosenfeld, A. " Computer Vision Systems” . Hanson, A.R., and Riseman, E.M. 
eds,.101. Academic Press, New York (1978). 

4. Duda, R.O., Reboh, R., et al (1979) "A Computer Based Consultant For Mineral Exploration” 
SRI Final Report For Project II 6415, Artificial Intelligence Ctr., SRI International, Menlo Park, 
California. 

5. Duda, R.O. , et al (1976) "Subjective Bayesian Methods For Rule-Based Inference Systems” 
Proc. Nat. Comp. Conf. pp. 1075-1082. AFIPS Press, Reston, Va. 

6. Garvey, T.D., (1976) "Perceptual Strategies For Purposive Vision” SRI International Tech. 
Note No. 117, Artificial Intelligence, Ctr, Menlo Park, California. 

7. Kleyle, R.M. and de Korvin, A. (1984), "Switching Mechanisms in a Generalized Information 
System”, Cybernetics and Systems. Vol. 15, pp. 145-167. 

8. Kleyle, R.M. and de Korvin, A. (1986), "A Two Phase Approach to Making Decisions Involving 
Goal Uncertainty”, Journal of Information Science. Vol. 12, pp. 161-171. 

9. Kleyle, R.M., de Korvin, A. "A Belief Function Approach To Information Utilization In 
Decision Making” Submitted. 

10. Pearl, J. "How To Do With Probabilities What People Say Y ou Can’t” The Second Conf. on 
Artificial Intelligence Appl. pp. 6-11, IEEE Comp. Soc. Press Miami, Florida (1985). 

11. Selfridge, P.G (1982) "Reasoning About Success and Failure In Aerial Image Understanding” 
Ph.D. dissertation, Computer Science Dept. Univ. of Rochester. 


10-14 



12. Schafer. G. (1976). A Mathematical Theory of Evidence . Princeton, NJ, Princeton University 

Press. 


13. Tanimoto, S.L., Computer Vision Systems . Hanson, A.R., and Riseman, E.M., eds Academic 
Press, New York (1978). 

14. Waltz. D.L. Computer Vision Systems . Hanson, A.R., and Riseman, E.M., eds Academic Press, 
New York (1978). 

15. Yakimovsky, Y., Feldman, J.A., (1973) in Proc. Third Int. Joint Conf. on Artificial Intelligence. 
580. The International Joint Conf. on Artificial Intelligence. 

16. Zadeh, L.A. (1986), "A Simple View of the Dempster-Shafer Theory of Evidence and its 
Implication for the Rule of Combination”, The AI Magazine, pp. 85-90. 

17. Zucker, S.W. in Pattern Recognition and Artificial Intelligence. Chen, C.H. ed Academic Press, 
New York. 


10-15 



ORIGINAL PAGE IS 
OF POOR QUALITY 


N89-20069 


EFFECTS ON MOTOR UNIT POTENTIATION AND GROUND 
REACTION FORCE FROM TREADMILL EXERCISE 

Final Repori: 

NASA/ASEE Summer Faculty Fellowship P rog ram- - 1 988 
Johnson Space Center 


Prepared By: Reid P. Elam, Ph.D. 

Academic Rank: Assistant Professor 

University & Department: Central Washington University 

Department of Physical Education 
Ellensburg, UA 98926 


NASA/ JSC 
Directorate: 
Division: 
Branch: 

JSC Colleague: 
Date Submitted: 


Space and Life Sciences 

Medical Sciences 

Space Biomedical Research 
Institute 

Victor Schneider, M.D. 
August 12 r 1988 


11-1 



ORIGINAL PACiL IS 

OF POOR QUALITY 

ABSTRACT 

Exercise countermeasures in long duration space flight 
have primarily been focused on two modalities, the treadmill 
and bicycle ergometer. Their use in space flight has been 
shown to offset the deleterious effects of disuse resulting 
from habituation in a weightless environment. Thus, muscle 
atrophy, osteoporosis and neuromuscular inefficiency - - 
physiological consequences from prolonged exposure to 
microgravity - - can be controlled by use of these modalities 
and their proper exercise regimen. Because muscular atrophy 
has been most significant in the posterior calf muscles and 
muscles of the thigh, the use of exercise modalities that 
concentrate on the musculature of the lower body is 
warranted. The treadmill and bicycle ergometer both require 
similar musculature in the activities of running and cycling 
although the sequence of their actions and intensity of their 
recruitment patterns may vary. Also, the reaction between 
the surface of a treadmill and the lower segments in ground 
reaction force (GRF) measures differ from pedal to foot 
reactions in cycling, running G R F 1 s being far greater than 
those measured in cycling. Motor unit potentiation (HUP) as 
measured by electromyography (EMG) is indicative of muscular 
actions and may be greater where stress from increased GRF 
is imposed. Therefore, greater MUP from treadmill exercise 
may be expected due to the greater GRF generated in running 
over cycling. This study was conducted to analyze the 
characteristics of MUP and GRF in treadmill exercise at the 
inclines of 0, 5.5 and 11 X with conjuctive speeds of 7.5, 6, 
and 5 MPH respectively. These speeds and corresponding 
inclines were set to provide equivalient physiological 
workloads at 12.5 METS. EMG recordings were taken from the 
rectus femoris and gastrocnemius of the right leg from 5 
subjects. Simultaneous GRF recordings were obtained from a 
Del mar Avionic treadmill rigged with load cells. Measures 
for MUP and GRF were taken over a period containing* 10 
strides at steady pace. Results of ANOVA revealed no 
significant differences between MUP of the rectus femoris at 
different speed/ incline settings. The same was true for the 
gastrocnemius. Significant differences were found between 
MUP of the rectus femoris and gastrocnemius at all treadmill 
settings (p<.01), with the MUP of the gastrocnemius being 
greater. No significant relationships between MUP of the 
thigh and GRF, and MUP of the calf and GRF were found (r=.097 
and r*-.31 respectively). Significant differences were found 
among GRF with 11X grade higher than OX (p<.05) and 5.5X over 
OX (p<.01). It was concluded that the gastrocnemius was more 
evident in EMG activity in all speed /incline settings over 
the rectus femoris, and that inclines from 5.5-11X produced a 
greater GRF's over OX. Recommendations for future studies 
was made. 
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INTRODUCTION 

Long duration space flight has a profound effect on the 
human musculoskeletal and neuromuscular systems- Reported 
consequences from prolonged habituation in a 
microgravi tational environment include muscular atrophy 
(2,4,5,7,12), bone loss <1,11, 17,1 (I) and decrements in 
neuromuscular performance (2,9,10,12,14,20). 

These physiologic problems are evident particularly in 
the lower body segments and arise from d i suse ( 1 , 1 4 , 1 5 ) - The 
absence of weight-bearing stresses to musculature and their 
supporting structures, as experienced in 1 g, lead to these 
maladies (1,11,15,16). Because the use of lower body 
segments in locomotion is no longer a requisite for 
ambulation, these structures are virtually devoid of 
volitional input. This creates in essence, a neuromuscular 
shunt, and a situation similar to what is observed in 
denervated muscle occurs (3,6). 

Exercise countermeasures have been reported to have 
some degree of efficacy in offsetting the debilitating 
effects of weightlessness (4,11,12,14,20). In particular, 
the exercise modalitites of primary use in manned space 
flights have been the treadmill and bicycle ergo me ter since 
these modalities concentrate their effort-requirements to the 
lower body (11,15,16). 

In previous studies, space flight has been shown to 
have an atrophying effect on muscles of the posterior calf 
(4,19,20) as well as the thighs and buttocks (19). Muscle 
strength deficits from post-flight studies have been observed 
in dynomometric measures of isokinetic force production of 
the lower limbs (14,20). E l ec t r omy og r a ph i c (EMG) studies 
showed that muscular efficiency and fatigue onset was 
increased after 1-14 days in space, and shifts to higher 
frequencies in power spectrum analysis following flights of 
longer than 2 weeks duration, implied muscle deterioration 
(12). Studies like these have supported the postulates 
concerning disuse atrophy, osteoporosis and neuromuscular 
degradat i on. 

Measurements comparing ground reaction or vertical 
forces of cycling and running have shown that ground reaction 
forces are greater in running than in cycling (1,15,16). 
Reasons for this difference center on factors of 
gravitational force, or forces imposed on the body and its 
supporting appendages. Basically, locomotion by running 
involves a great amount of decelerative or negative forces 
which provide for the storage of elastic energy to be 
released in subsequent positive muscular contractions (13). 
This was termed the st retch- shorteni ng cycle by Komi (8) and 
hints at the existence of eccentric loading in phases of 
human locomotion (eg. running) which may be absent in non- 
weight bearing activities (13). Therefore, the greater 
reaction forces between body segments and surface impact from 
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the effects of gravity, the greater the requirements for 
muscular activity. Subsequent ly, it is the increased ground 
reaction forces that provides the implication of greater 
muscular activity involved in running over cycling. 

The purpose of this study was to: 

1. make a basic 1 g analysis of motor unit potentiation and 
ground reaction force during treadmill running at different 
speeds and inclines. 

2. study any relationships between motor unit potentiation 
and ground reaction force during treadmill running at 
different speeds and inclines. 

3. establish a means for assessment of individual 
force/muscle activation characteristics for implementation 
of protocol to other studies, to include those for 
t r a i n i ng/ t es t i ng in a weightless environment. 

REVIEW 


Motor unit potentiation. 

The temperosopatial aspects governing human 
movement are derived from the functioning of the basic 
components of the integrated neural and muscular systems. 
Together these systems comprise the operative element, the 
motor unit. It is the motor unit, when activated via an 
electrochemical sequence that controls muscular contractions 
and subsequent limb movement. The force and speed at which 
limb movements occur are dictated by several factors as 
f o l lows : 

1. the number of motor units recruited. 

2. the frequency of firing of motor units. 

3. the synchrony of motor unit firing. 

4. the pattern of motor unit recruitment. 

5. the type of motor units activated. 

6. the inhibition of antagonist motor units. 

These factors are applicable as determinants of the amount 
muscular activity produced during running and cycling. 

ENG is a means by which the electromechanical dynamics 
of the neuromuscu l ar system can be measured. Activation of 
muscular contractions are measured by ENG as electrical 
current is generated from the difference in potential sources 
of a dipole. The potential differences at each pole reflect 
the electrical activity to ensue, and the term motor unit 
potentiation is used to characterize this. 

Ground reaction force. 

Forces that the body exerts against the earth's surface 
in movement, in line with Newton's third law, can be measured 
by low frequency output from piezoelectric sensors built into 
load cells. Changes in the forces exerted against the 
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running surface of the treadmill are transduced into electric 
signals which can be calibrated to units of G force. 

Often, ground reaction forces are similarly termed 
vertical or transient forces as the line of gravity is 
vertical and transmitted along the longitudinal axes of bone 
structures (transient) in upright ambulation. 

METHODS 

Subjects were selected from among volunteers who were 
screened for CHD, were less than 40 years of age, had a 
completed Air Force Class III physical examination, and were 
familiar with treadmill exercise (M = 5). Following a 
complete orientation, subjects signed an informed consent 
document . 

Preparation of the subjects for E MG recording, entailed 
the use of alcohol prep pads to cleanse and abrade the skin 
for placement of surface electrodes. The skin was rubbed 
until it became e r y t h r oma t ous and has an electrical 
resistance less than 10,000 ohms. 

The rectus femoris and gastrocnemius of the right leg 
were prepared for electrode placement. Prior to placement 
of the electrodes, a measurement was taken from the point of 
origin of each muscle to its insertion with the midpoint 
cited as the point of innervation. Two Ag-AgCl electrodes 
were placed on either side of this point so that their 
centers were 4 cm apart. Leads from these electrodes were 
inserted into a signal conditioner worn at the waist of the 
subject. The input EMG signals were amplified and 
transmitted by infrared telemetry to a signal receiver 
capable of receiving simultaneous signals. These signals 
were transferred via BNC connection to a Compaq 40 MB hard 
disk computer for A/D conversion and storage for later 
analysis. 

Ground reaction forces were recorded from four load 
cell transducers located under the belt of a Delmar Avionics 
motor-driven treadmill. Analog signals were sent via 
hardwire BNC connection to the Compaq computer for A/D 
conversion and storage for future analysis. 

Two signal receivers were placed on the wall on the 
side of the treadmill such that the field of signal 
acceptance would be perpendicular to the emitted signal. 

Each subject was required to perform treadmill running 
in one session at three different speeds and inclines. 
Following a warm-up period including stretching exercises and 
4 minutes of low intensity treadmill running, the treadmill 
speed was randomly (as respects order) increased to settings 
of 5, 6 and 7.5 MPH. These speeds were set conjunctively 
with grades of 11, 5.5 and 0 X respectively. These speeds 
and grades (inclines) were selected to represent a 
physiological equivalency in workload at 12.5 METS (1 MET * 
V02 of 3.5 ml02/kg/min at rest). 
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The subjects were instructed to maintain pace at 1-2 
minute intervals at the prescribed speeds and inclines such 
that EMG and ground reaction forces could be taken for a 
segment of 10 strides. 

After completion of the exercise, a cool-down period of 
5 minutes was administered and consisted of low intensity 
walking at a OX grade. This was followed by a 10 minute 
recovery period with the subject in a seated position. 

Analysis. 

Analysis of the raw ENG was performed by integration, 
and averaged over the ensemble recordings for the 10 strides 
to produce a mean value for motor unit potentiation. 

Measures of ground reaction force was analyzed for the 
l and i ng - suppor t - t a keof f phase of the running stride. An 
ensemble average of the G force was calculated over 10 
strides. 

Averaged values for motor unit potentiation and ground 
reaction forces were compared intersub jectly to note 
differences between variables in prescribed speeds and 
inclines of the treadmill. Correlations for motor unit 
potentiation and ground reaction force were made to analyze 
relationships between the two variables. 


RESULTS 

The results of the study are contained in Table 1, and 
a sample recording of G forces and EMG shown in full wave 
rectified form are shown in Figure 1. 

Results from ANOVA showed that there were no 
significant differences in motor unit potentiation of the 
quadriceps at all speeds and incline settings (p<.05). The 
same was found to be the case in differences between motor 
unit potentiation of the gastrocnemius at the various speed 
and incline settings. There was a significant difference 
between motor unit potentiation of the quadriceps and 
gastrocnemius at all speeds and incline settings (p<.01), the 
mean integrated EMG < 1 EMG ) being greater than quadricep I E MG . 
Significant differences existed in ground reaction forces 
[F(2,12)=12.52;p<.01] . In post-priori comparison t tests, 
ground reaction forces generated at 5.5X grade and 6 MPH 
speed, and 11% grade and 5 MPN speed, were greater than at OX 
grade and 7.5 MPH speed <p<.01 and p<.05 respectively). 
Ground reaction forces between 11X/5 MPH and 5.5X/6 MPH 
showed forces to be slightly higher at 5.5X/6 MPH, however 
the difference was not significant (p<.05). 

The correlation between motor unit potentiation of the 
quadriceps and ground reaction forces was not significant 
(r».097), nor was the correlation between gastrocnemius motor 
unit potentiation and ground reaction force (r»-.31). This 
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TABLE 1. RESULTS OF TREADMILL RUNNING ON MOTOR UNIT 
POTENTIATION (NUP) AND GROUND REACTION FORCE (GRF) FOR 
DIFFERENT TREAOMILL (TM) SPEEDS AND GRADES. VALUES ARE MEANS 
CN * 5). 


TN (MPH/X) 


HUP ( uVs ) 


GRF (G) 

Speed/ | 

incline f 

Rectus 
f emo r i s 

Gastroc- 

nemius 

TOTAL | 

1 

T ouchdown/ 
takeoff 

i 

5/11 | 

1 

166.164 

372.438 

538.602 j 
1 

2.016 

1 

6/5.5 | 

1 

203.838 

302.598 

1 

506.436 | 
1 

2.104 

1 

7.5/0 1 

226.968 

344.258 

1 

571.226 1 

1.816 
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Figure 1. Example from one subject for (a) force recording 
and full wave rectified EMG recording for (b) the 
gast ocnenti us and (c) rectus femoris at 11 X grade and 5 MPH 
treadmill setting. 
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was for all speed and incline settings combined. 

CONCLUSION 

The fact that impact or collision forces exist between 
the surface of the treadmill and the feet of the exerciser, a 
vertical force is transmitted along the longitudinal axes of 
the bone structures. As long as this force is present, the 
resulting transarticular reaction forces, particularly at the 
t a l of i bu l a r/ t i bi a l and femorotibial surfaces, will require 
added supporting contractions of the surrounding musculature 
for joint stabilization. The demands imposed on the 
musculature of the thigh in knee stability and of the shank 
in ankle stability may be increased as the reaction forces 
become greater due to changes in speed and incline. This was 
not found to be the case in this study, as no significant 
differences were found between EMG activities of the rectus 
femoris at all speeds and inclines. The same was true for 
the gastocnemius. Also, no correlations were found between 
ground reaction forces and motor unit potentiation of either 
the rectus femoris or gastrocnemius. Significant differences 
were found between potentiation of the gastrocnemius and 
rectus femoris at all speeds and inclines which showed that 
the calves are more active relative to the quadriceps in all 
cases analyzed in this study. 

Because running requires the individual to be airborne 
in between strides, the resulting G forces upon impact will 
be greater than in walking. That increases in speed require 
a greater pushoff in the takeoff phase, resulting in greater 
airborne time, the ensuing impact may be thought to be 
greater, however the contrary was found to be true from the 
results of this study. Greater G forces were found at 
inclines of 1 1 X and 5.5X over horizontal running at speeds of 
5 and 6 NPH respectively versus the greater speed of 7.5 HP H 
at OX grade. This was attributed to the greater horizontal 
impedence necessary to be exerted against the belt for 
pushoff in the horizontal direction. This detracts from the 
vertical pushoff force that was measured by the load cell 
i ns t rument a t i on . 

Greater EMG activity of the gastrocnemius at 11X grade 
may be attributed to the need for more forefoot running, 
therefore eliciting a greater stretch on the calf muscles and 
resulting in a stronger muscular impulse. This supports 
Koii 's (8) assertion of eccentric factors contributing to the 
stretch-shortening reflex (stronger ensuing positive 
contraction) in running. It was interesting to note that the 
potentiation of the gastrocnemius at 5.5X was lower than both 
the 1 1 X and OX grades. This would show that either increased 
speed or inclines at 1 1 X or higher may be capable of 
eliciting higher EMG responses than at lesser speeds at mid- 
range grades -- perhaps 1-10X. 

Ground, react i on forces were greater at inclines of 1 1 X 
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and 5.5X over OX, however the higher values at 5.5X over 11X 
indicated that there may be an optimal grade for promoting 
higher G forces in running at moderate speeds. 

In summary, the results of this study revealed that 
speed and incline variation in treadmill exercise at 1 g can 
effect ground reaction force such that higher G forces can be 
produced at grades ranging from 5.5 to 11X. Motor unit 
potentiation of the quadriceps and gastrocnemius were found 
to be virtually unaffected by speed and incline variations as 
respects significant difference between the changes in these 
variables, although trends may exist between increased 
speed/ i nc l i ne and increased EMG activity, but this was not 
analyzed in this study. 

Inclines may have no effect in a 0 g environment, 
however, restraint systems to treadmill exercises can be set 
such that forward inclinations in body lean can be made to 
produce resultant forces to bring about greater ground 
reaction forces. Increased G forces imposed by treadmill 
parameters imply greater transient forces to long bones 
involved in locomotion and therefore greater benefit to the 
structural integrity of bone tissue. 

The study also revealed that physiological workload 
equivalencies, as set by MET requirements, are not paralleled 
by biomechanical and neuromuscular factors in actual exercise 
on the treadmill. Hence exercise countermeasures in space 
flight, where muscle, bone and neuromuscular, considerations 
are concerned necessitate monitoring of the effect of the 
exercise regimens and modalities according to the response of 
the human body to the imposed demands of exercise by EMG and 
force transduction. 

As this serves as a pilot study, the results Lead to 
the following recommendations: 

1. A basic 1 g comparison between the effects on motor unit 
potentiation (MUP) and ground reaction force (GRF) in 
treadmill and bicycle ergometer exercise. 

2. The effects of prolonged treadmill running on bone 
density versus the same in bicycle ergometer exercise. 

3. An analysis of the muscular torques and transarticular 
forces about the knee and ankle in treadmill running by means 
of optoelectrical data collection and anthropomorphic 
mode ling. 

4. The effects of treadmill running at various body 
inclinations and speeds on MUP and GRF in 0 gravity. 

5. The differences between total body standing resistive 
exercises and treadmill running on MUP, GRF, muscle torques, 
transarticular reaction force and bone density. 

6. The effects on MUP and GRF in decline treadmill running: 
an analysis of eccentric contractions in locomotion. 

7. The effects of running shoes on MUP and GRF at different 
speeds and inclines in treadmill running. 

8. The re l a t i onsh i ps between increasing speed/ i nc l i ne on HUP 
and GRF. 
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9. The relationships between HUP and physiological measures 
of energy expenditure (ie. V02 and kcal). 
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ABSTRACT 


7 -3 

Particles in the mass range from 10 to 10 grams. contribute 80% of 
the total mass influx of meteoritic material in the 10' 1 '* to 10 b gram mass 
range at Earth (Hughes, 1978) _ On Earth atmospheric entry all but the 
smallest particles in the 10” ' to 10~ J gram mass range, about 60 to 1200 
micrometers in diameter, are heated sufficiently to melt of vaporize. Mars, 
because of its lower escape velocity and larger atmospheric scale height, 
is a much more favorable site for unmelted survival of micrometeorites on 
atmospheric deceleration. We calculate that a significant fraction of 
particles throughout the 60 to 1200 micrometer diameter range will survive 
atmospheric entry unmelted. Thus returned Mars soils may offer a resource 
for sampling micrometeorites in a size range uncollectable in unaltered 
form at Earth. 

The addition of meteoritic material to the Mars soils should perturb 
their chemical composition, as has been detected on for soils on the Moon 
(Anders, et al., 1973). Using the measured mass influx at Earth and 
estimates of the Mars/Earth flux ratio, we estimate a mass influx at Mars 
of between 2,700 and 202,000 metric tons per year. If distributed uniformly 
into a soil with a mean planetary production rate of 1 meter per billion 
years, consistent with radar estimates of the soil depth overlaying a 
bouldered terrain in the Tharsis region (Christensen, 1986), we estimate a 
meteoritic concentration in the soil ranging from 2 to 58%. This range is 
consistent with chemical inferences that Mars soil at the Viking sites is 
fit by a mixture of 60% indiginous rock fragments and 40% meteoritic 
material (Clark and Baird, 1979) . 
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INTRODUCTION 


A number of sources for dust on Mars have been proposed including: 
chemical and physical weathering of indiginous surface rocks , debris from 
impact cratering, volcanism, tectonism (Greely, 1986) , and weathering of 
shock- activated meteoritic projectiles from the period of heavy bombardment 
(Boslough, 1988). However, the contribution from micrometeorites , particles 
sufficiently small to be slowed without melting on atmospheric 
deceleration, has only recently been considered (Flynn and McKay, 1988). 
Particles in the size range which would be expected to survive atmospheric 
entry on Mars without melting constitute a major fraction of the meteoritic 
material presently accreting onto the Earth. On Earth, the rate of 
production of new soil is sufficiently high that unmelted particles are 
difficult to extract from the soils, however the meteoritic component can 
be detected by enrichment in Ir (Kyte and Wasson, 1986) , and the particles 
themselves can be recovered from the stratosphere (Brownlee, 1985) and 
polar ices (Zolensky, et al., 1988), where the terrestrial soil 
contamination is reduced. Partially melted particles can also be extracted 
from sea sediments by magnetic concentration (Brownlee, 1985). Although 
particles in this size range are destroyed on Lunar impact, the Lunar mare 
soils contain between 1 and 2 percent meteoritic material, inferred from 
chemical abundances of the volatile and siderophile elements which are 
abundant in primitive meteorites but rare in the mare rocks (Anders, et 
al., 1973). 

Mars, because of its low surface gravity coupled with an atmosphere of 
sufficient density to decelerate incoming micrometeorites, is one of the 
most favorable sites in the solar system for the unaltered survival of 
micrometeorites on atmospheric entry. The magnitude of the micrometeorite 
component in the Martian soils will depend on the accretion rate of 
meteoritic material onto the planet, as well as the rate at which soil is 
being produced by indiginous processes on Mars. To assess the importance of 
the micrometeorite component in the soil on Mars, we have: 

1) assessed the micrometeorite flux at Mars, 

2) estimated the micrometeorite velocity distribution at Mars, and 
evaluated the survival probability for micrometeorites entering the Martian 
atmosphere , 

3) considerd the chemical perturbations to the indiginous soil produced 
by micrometeorite addition, and suggested measurements which a Mars Lander 
spacecraft could make to determine the fraction of meteoritic material in 
the soils, and, 

4) assessd the possibilities for recovery of micrometeorites from the 
soils returned by a Mars Sample Return mission, and. 

Four types of meteoric and/or meteoritic material should be found on 

Mars: 

1) micrometeorites, many of which will survive atmospheric deceleration 
unmelted, which should fall relatively uniformly over the planet's surface, 

2) ablation products from larger meteors and meteorites which ablate, 
break up and/or burn up in the Mars atmosphere, 

3) debris from large, crater forming objects, which, by analogy to 
terrestrial and lunar impact events, will be concentrated in the crater 
ejecta blankets (except for rare, large events, such as the proposed C-T 
event on earth, which can distribute debris on a planetary scale) , and 

4) debris from the early, intense bombardment, which, in many areas of 
the planet, may now be incorporated into rocks by geologic processes 
subsequent to the intense bombardment era. 


FLUX AT MARS 

To estimate the extent of meteoritic addition to indigenous Martian 
soil, the meteoritic flux at Mars must be known. The meteoritic flux 
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measured at Garth provides a starting point to estimate the flux at Mars. 
For particles in the mass range from 10* grams to 10° grams, the size- 
frequency distribution at Earth has been estimated from satellite, radar 
meteor and visual meteor observations (Hughes, 1978). Fifty percent of the 
total mass influx is in the narrow mass range from 10*° grams to 10" 4 grams 
(Hughes, 1978), and eighty percent is in the range from 10*' to 10*° grams. 
The Earth mass influx per decade of mass is shown in Table I. If the size- 
frequency distribution of meteoritic material at Mars is similar to that at 
Earth, onlv 0.53 au away, then an assessment of the Mars flux for material 
in the 10*' gram to 10'° gram mass range will provide a good estimate of 
the overall influx of meteoritic material on Mars. _ 

,On entering the Earth's atmosphere, most particles in the 10*' gram to 
10" ° gram mass range are volatilized, producing ionized trails detectable 
by radar, while some at the lawest end of this mass range are decelerated 
sufficiently slowly to survive entry unmelted (Brownlee, 1985). Veriani 
(1973) obtained the density distribution of over 5,000 radar meteors in the 
10*° gram to 10*° gram, mass range. More. than 50% of his meteors fell in the 
density range of 10* *° gm/cm° to 10 +U- gm/cm , with a mean of the density 
distribution at 0.8 gm/cm° . Taking 1 gm/cm° as a representative density for 
the particles, the 10*' gram to 10*° gram mass range corresponds to 
particles from 60 to 1200 micrometers in diameter. Particles near the lower 
end of this mass range, up to about 100 micrometers in diameter, have been 
collected from the Earth's stratosphere, after atmospheric deceleration 
(Brownlee, 1985). The properties of these micrometeorites, collected in the 
NASA Cosmic Dust Sampling Program, have been reviewed by Brownlee (1985) . 

Most of the material in the 100 to 3000 micrometer size range incident 
on the top of the Earth's atmosphere would be expected to melt of 
volatilize during atmospheric deceleration (Brownlee, 1985). However Mars 
is a much more favorable site for particle survival on atmospheric 
deceleration since: 

1) the scale height of the atmosphere on Mars is larger than on Earth, 
resulting in a longer deceleration time (and lower peak temperature on Mars 
than on Earth) . and 

2) the lower planetary mass and surface gravity give rise to less 

gravitational infall acceleration, thus particles with similar in- space 
velocities prior to planetary encounter will enter the atmosphere on Mars 
with lower velocity than on Earth. , 

The current mass influx at Earth in the 10* gram to 10° gram mass 
range is estimated, from satellite, radar meteor, and visual meteor 
observations , ,at 16,100 tons per year, with 13,000 tons per year in the 10* 
' gram to 10* ° gram mass range (Hughes, 1978). The measurement of Ir, which 
is believed to be a signature of extraterrestrial material, in atmospheric 
particulate samples at the South Pole suggest a current flux of 
extraterrestrial material of 11,000 tons per year (Tuncel and Zoller, 

1987), consistent with the mass influx obtained by Hughes (1978). 

Measurements of the long term meteoritic influx, by Ir concentrations 
in Pacific ocean sediments, have given values higher than the current flux. 
Kyte and Wasson (1986) infer a relatively constant mass influx (except for 
sharp spikes corresponding to major impact events) of 78,000 tons per year 
over the past 67 million years from Ir in Pacific sedmiments. This is 
consistent with an earlier value of 110,000 tons per year derived from 
Pacific sediments by Shedlowsky and Paisley (1966) . We will perform our 
calculations using both the current Hughes (1978) mass influx at Earth, and 
the Kyte and Wasson (1986) long-term mass influx, which is a factor of five 
higher than the Hughes value. Literature values for the extraterrestrial 
mass influx at Earth, which are tabulated in Tuncel and Zoller (1987), 
exceed even this range. Since the Ir measurements provide no information on 
the incoming size distribution of the particles, we will assume the current 
size distribution from Hughes (1978) is also representative of the size 
distribution of the meteoritic material contributing to the long-term flux 
determined by Ir concentration (Kyte and Wasson, 1986) . Some confirmation 
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of this assumption is provided by the analysis of the size versus frequency 
distribution of microcraters on exposed lunar rock surfaces returned by the 
Apollo missions. For particles larger than 10' y grams, where the effects of 
secondary impacts are no longer believed to be significant, the shape of 
the long-term lunar microcrater size versus frequency distribution is 
generally consistent with the mass distribution inferred from current 
satellite and meteor measurements (Grun, et. al, 1985). 

To extrapolate the mass influx at Earth to the corresponding Mars value 
requires an estimate of the ratio of the Mars flux to the Earth flux. This 
ratio depends on the type of orbital evolution experienced by the 
particles . The orbits of large meteorites are perturbed principly by 
gravitational interactions with the planets. However, for small particles 
solar radiation pressure forces cause significant orbital perturbations on 
time scales comparable to or shorter than the gravitational perturbation 
time scale (Dohnanyi, 1978). 

LARGE OBJECTS 

For large, crater producing, objects whose orbits are dominated by 
gravitational perturbations, the relative crater production rates on the 
terrestrial planets and the Moon have been assessed to establish 
chronologies. Anders and Arnold (1965) have estimated the meteoritic input 
on Mars to be 25 times the Lunar value, while Soderblom, et al. (1974) have 
estimated the input of meteoritic material on Mars to be only twice that on 
the Moon. Shoemaker (1977) estimated that the ratio of impact rates of 
bodies to absolute visual magnitude 18 on Mars and Earth is 2.6. When 
corrected for differences in impact velocity, and scaled for the planetary 
surface gravities, Shoemaker (1977) arrived at a cratering rate ratio of 
1.6. Hartmann et al. (1981) have reviewed the planetary cratering rate 
estimates, and they adopt a best value for the crater production rate on 
Mars of 1.3 that for Earth. Hartmann et al. (1981) indicate there are 
factor of 2 uncertainties in the cratering rate ratio, however, since the 
proportions of objects in various types of Earth and Mars crossing orbits' 
are not well established. Since the Hartmann, et al. (1981) best value of 
the cratering ratio is consistent with the Shoemaker (1977) value, we will 
adopt the Shoemaker (1977) impact rate ratio of 2.6 as indicative of the 
ratio of the mass influx for large, gravitationally perturbed, objects at 
Mars and Earth. 

SMALL OBJECTS 

Objects whose orbits are dominated by F-R drag may, however, have a 
different ratio of the Mars to Earth flux. For particles up to 100 
micrometers in diameter, the dominant radiation effect is a drag force, 
Poynting- Robertson (P-R) drag (Dohnanyi, 1978). For objects starting from 
circular orbits, P-R drag causes them to spiral into the sun. If the 
initial orbit of the particle is elliptical, P-R drag cause a rapid 
decrease in the perihelion and a slower decrease in the aphelion, so that 
the ellipticity of the orbit decreases as the particle falls into the sun. 

The initial orbits of the dust particles will be determined by the 
orbits of the parent objects. The Infrared Astronomy Satellite (IRAS) 
detected two major types of solar system dust sources: the main belt 
asteroids (Low, et al., 1984), and comets (Sykes, et al., 1986). The dust 
bands detected in the main asteroid belt are thought to have been produced 
by low velocity collisions between main-belt asteroids (Sykes and 
Greenberg, 1986). Small particles produced in such collisions would spiral 
in towards the sun under P-R drag. They would pass Mars, providing an 
opportunity for Mars collection, and later pass Earth. Thus the Earth flux 
could be used to estimate the Mars fluut for particles from this source 
region. 

In the case of the particles emitted by comets, those detected by 
Sykes, et al. (1986) all had aphelia outside the orbit of Mars, and thus 
their orbits will also evolve, under P-R drag, through Mars collection and 
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subsequent Earth collection opportunities . Only if significant dust sources 
existed between the Earth and Mars, would there be a catagory of particles 
which, under P-R drag would be collectable at Earth but not at Mars. No 
such sources were reported in the IRAS survey. 

To assess the ratio of the particle flux at Mars to that at Earth, the 
velocity distribution of the particles must be known at both planets. As 
shown by Opik (1951), the effective planetary capture cross-section, <r , 
varies with the velocity of the incident particle, Vp, as: 

T - yRp 2 ’(l + v e 2 /v p 2 ) (Equation 1) 

where v 0 is the planetary escape velocity, and R^ is the planetary radius. 

-The velocity distribution at Earth for radairmeteors , particles in the 
10" 6 to 10" 2 gram mass range, has been determined by Southworth and 
Sekanina (1973) for a set of over 14,000 radar meteors. We have used the 
Zook (1975) approximation to the Southworth and Sekanina (1973) atmospheric 
entry velocity distribution at Earth: 

F(v) - 3.822 x e '°- 2468v (Equation 2) 

where F(v) is the fraction of the particles having a velocity v. This 
distribution is shown in Figure 1. We then followed the same procedures 
used by Morgan, et al. (1988), who calculated the meteoritic velocity 
distribution at Mercury, to calculate the velocity distribution at Mars. 
First, the Earth entry velocity distribution was corrected to an in- space 
distribution at 1 au by removing the near-earth gravitational focusing in 
each velocity increment, and removing the effect of earth infall 
acceleration. Next, the resulting velocity distribution was transformed to 
1.53 au. At this time, the difference in flux at 1.53 au was also accounted 
for, taking the flux fall off with increasing heliocentric distance to vary 
as r" 1 , as determined from zodiacal light observations (Hanner, et al., 

1976; Schuerman, 1980). We then transformed the space velocity distribution 
at 1.53- au to a Mars atmospheric entry velocity distribution taking into 
account both the Mars gravitational focusing effect and the gravitational 
infall acceleration. The resulting atmospheric entry velocity distribution 
is shown in Figure 1. The ratio of the area under the curve of Earth entry 
velocity, which was normalized to 1.0, to the area under the Mars velocity 
distribution curve (0.57) is the flux ratio. This flux ratio must be 
multiplied by the planetary cross-sectional area ratio, equal to 0.29, to 
obtain the mass influx ratio. Thus, we estimate that, for a population 
dominated by P-R drag, the ratio of the mass influx at Mars to that at 
Earth would be 0.17. 

The time for an equivalent amount of orbital evolution under P-R drag 
increases lineraly with particle diameter if the density and reflectivity 
remain constant. Thus, as particles become larger they are more likely to 
be gravitationally perturbed from their orbits or destroyed by catastrophic 
collisions before P-R drag produces a significant change in the orbit. 
Estimates of the catastrophic collision lifetimes for particles in the main 
asteroid belt (Dohnanyi, 1978) and in space at 1 au (Grun, et al., 1985) 
both suggest that particles larger than about 100 micrometers in diameter 
will be destroyed by collisions before their orbits are perturbed into the 
sun. Uncertainties in the flux of particles producing collisional 
destruction might change this size limit. 

FLUX ASSESSMENT 

Thus particles larger than 100 micrometers diameter, which contribute 
75% of the mass influx at earth, are expected to be only slightly larger 
than the size range whose orbital evolution is dominated by P-R drag. The 
actual ratio of the Mars to Earth flux is thus likely to be somewhere 
between that which we have calculated for P-R drag dominated particles, and 
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the value obtained by Shoemaker (1977) for larger objects whose orbital 
evolution is dominated by gravitational perturbations. 

In assessing the meteoritic infall on Mars, we will take as a lower 
estimate the a value of 2700 tons per year obtained using the Hughes (1978) 
earth flux coupled with our P-R drag dominated flux ratio. An upper 
estimate of the Mars infall, of 202,000 tons per year, was obtained using 
the Kyte and Wasson (1986) terrestrial flux and the (Shoemaker, 1977) flux 
ratio. The large range in values reflects the uncertainties in the Earth 
flux and the Mars/Earth flux ratio. 

If distributed uniformly over the planet, the corresponding mass 
accretion rates range from 18 grams per square meter per million years to 
1350 grams per square meter per million years. For density 1 gm/cnr 
material, these correspond to the addition of between 1.8 cm/billion years 
and 1.35 meters/billion years of meteoritic material to the Martian 
surface. Even the lower estimate would be expected to produce some 
detectable perturbations in the soil composition, unless the indiginous 
materials are very similar to meteoritic composition or the soil production 
rate far exceeds that on Earth. 7 

Of this infalling mass, about 80% would be expected to be in the 10* / 
gram to 10 gram mass range. 


MICROMETEORITE SURVIVAL 

Deceleration in the -atmosphere of the Earth heats most of the particles 
in the 10~ b gram to 10' J gram mass range above their melting temperature. 
Some particles in this mass range are recovered as melted spherules from 
the ocean bottom (Brownlee, 1985) and from pools in Greenland (Maurette, et 
al., 1988). However no extraterrestrial particle larger than approximately 
50 micrometers in diameter has been recovered unmelted, and intact in the 
NASA Cpsmic Dust Collection Program (Zolensky, pers. comm.). 

Brownlee (1985) estimates the melting temperature for micrometeor itic 
material at about 1600 K. The distribution of peak temperatures reached by 
a micrometeorite on atmospheric entry can be predicated using an entry 
heating model developed by Whipple (1950) and extended by Fraundorf (1980) . 
This temperature distribution is generally expressed as the fraction of 
particles heated above any given temperature. We show in Figure 2, the 
calculated fraction of incident particles of diameter 10 microns, 100 
microns, and 1000 microns, heated above temperatures from 300 K to 2000 K 
on Earth atmospheric entry. The calculations use the Fraundorf (1980) 
equations, and assume the Southworth and Sekanina (1973) entry velocity 
distribution, a particle density of 1 gm/cm , and a value of 1 for the 
parameters charactarizing drag, kinetic energy transfer, and emmissivity. 
The reasons for this latter assumption are examined in Flynn (1988) . 

The general validity of the Fraundorf (1980) entry heating model can be 
verified by noting that, if 1600 K is taken as the critical temperature for 
unmelted survival, essentially all 10 micron particles, about half of the 
100 micron particles, and essentially none of the 1000 micron particles 
would be expected to survive Earth atmospheric entry unmelted. This is 
consistent with the size cutoff on unmelted particles recovered from the 
stratosphere, and the observation that larger meteoritic material recovered 
from sediments is mostly melted. 

We have applied the Fraundorf (1980) entry heating model to the Mars 
atmosphere. Within 30 km of the surface, the atmospheric scale height 
measured by Viking 1 was 11.70 km, and by Viking 2 was 11.36 km (Seif and 
Kirk, 1977). However spacecraft deceleration data for the atmosphere from 
30 km to 120 km (Seif and Kirk, 1977) give a scale height in this region of 
7.9 km. Since particles in the size range of interest reach their maximum 
deceleration in the upper region, we have used a scale height of 7.9 km. 

The velocity distribution which we inferred for Mars entry (Figure 1) was 
used, and the particle size, density, and interation parameters were the 
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same as for the Earth atmospheric entry calculations. These results are 
also shown in Figure 2. 

Unlike the Earth case, almost 90% of the 100 micron diameter particles 
would be expected to survive Mars atmospheric entry without melting. For 
the 1000 micron diameter particles, 30% would be expected to survive 
atmospheric entry unmelted. Micrometeorites in this size range which have 
survived atmospheric entry are rare in the Earth collections. Since these 
larger particles from Mars orbital distance are likely to sample different 
sources than the smaller micrometeorites collected at 1 a.u. in the cosmic 
dust sampling program on Earth (Flynn, 1988; Zook and McKay 1986), returned 
Mars soils may provide a unique resource for micrometeorite analysis. 

We have used our calculated Mars atmospheric entry velocity 
distribution, and the Fraundorf (1980) entry heating model, to plot the 
predicted fraction of incoming micrometeorites not heated above temperature 
T on Mars atmospheric deceleration for temperatures ranging from 700 K to 
1900 K. These results are shown for particles diameters from 10 to 1000 
micrometers in Figure 3. We have used this survival fraction (those not 
heated above 1600 K at each size), coupled with the size -frequency 
distribution of micrometeorites from Hughes (1978) , and our high and low 
estimates of the total mass influx at Mars to calculate the micrometeorite 
addition rate (particles per square meter of surface per year) to the 
Martian soils. Since each decade of mass spans only a factor of 2.1 in 
particle diameter, and since the particle abundance is a rapidly decreasing 
function of mass, all particles within each mass decade are taken to be at 
the smallest diameter in that decade. Particles are assumed to have a 
density of 1 gm/cm , consistent with the range of 0.7 to 2.2 gm/cnr 
measured for the smaller micrometeorites recovered from the Earth's 
stratosphere (Fraundorf, et al., 1982). These results of these calculations 
are reported in Table II. The expected abundance of micrometeorites in a 
returned Mars soil sample could be estimated from these results if the soil 
production rate on Mars were known. 

Estimates of the thickness of the Martian regolith vary widely, from 
only twice as thick as the lunar regolith (Soderbloom, et al., 1974) to as 
deep as 2 km (Fanale, 1976). However, much of the planetary regolith was 
very likely generated during the intense bombardment era. Depending on the 
mixing depth, the present meteoritic infall may or may not be mixed into 
the soil of that early regolith. 

Various physical properties of the Martian surface can be used to 
constrain the thickness of the current dust deposits. The low thermal 
inertia of the deposits requires a minimum thickness of order 0.1 meters 
(Harmon, et al. , 1982). Harmon, et al. (1982) also suggest that the 
presence of exposed rocks and the degree of visible mantling indicate the 
dust thickness is less than 5 meters. Dual -polarization radar measurements 
in the Tharsis region indicate a rough texture, which suggests that a 
relatively thin dust layer covers near- surface rocks. Based on radar 
penetration properties, Christensen (1986) estimates a dust mantle 
thickness of only 1 to 2 meters. Arvidson (1986) suggests that most of the 
sedimentary debris on Mars was produced relatively early, perhaps in the 
first billion years. In more recent times, the preservation of a large 
number of pristine -looking, small bowl shaped craters at the Viking 1 
lander site suggests a rate of rock breakdown and removal of only meters 
per billions of years (Arvidson, 1986). 

Thus, while the regolith itself could be quite deep (Fanale, 1976), 
much of it is likely to have been produced during the intense bombardment 
of the planet. The surface soil into which the last billion years of 
meteoritic material may be concentrated, could be only a few meters deep. 

DIRECT COLLECTION 

We calculate the expected concentration of micrometeorites in a soil 
whose planetary average production rate is 1 meter per billion years . The 
concentrations obtained can then be scaled to other assumed soil production 
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rates by multiplying by (1 meter/billion years) /(production rate) . Table 
III shows the number of unmelted micrometeorites in each size range 
expected in an average 10 gram Mars soil sample. These results suggest that 
returned Mars samples may offer a new resource for the study of 
micrometeorites. These micrometeorites may sample a different parent 
population than the smaller particles recovered from the Earth's 
stratosphere. The survival lifetime of micrometeorites deposited in the 
harsh Martian enviornment is unknown, however on Earth millimeter size 
fragments of meteoritic material, both unmelted and melted, have been 
recoved from late Pliocene deep sea sediments (Brownlee, et al. , 1982), 
demonstrating that survival is possible on Earth for several million years. 

The soils would also be expected to contain melted micrometeorites in 
the larger sizes. On Earth, such particles can easily be extracted from the 
ocean bottom by magnetic separation, due to the formation of magnetic 
minerals on atmospheric entry (Brownlee, 1985). Depending on the abundance 
of indiginous magnetic minerals in the Mars soil, a magnetic separation may 
permit recovery of the melted particles to determine if the meteoritic 
concentration is sufficient for search and extraction of unmelted 
particles . 

COLLECTING SITES FOR MICROMETEORITES 

Martian surface processes (weathering and wind erosion, transport, and 
deposition) may fractionate the dust by size, density or composition 
providing regions of increased local concentration, suggesting even more 
suitable sites for micrometeorite sampling than the average soil. These 
sites may include placer catch basins or lag surfaces which may accumulate 
high density micrometeorites or their derived and altered minerals. 
Conversely, low density micrometeorites may be wind segregated along with 
finer Martian dust and may constitute a relatively coarse-grained component 
of that dust at its deposition sites. By analogy with Antarctica, 
meteorites of all size ranges may be relatively concentrated in Martian 
polar regions, although the concentration mechanisms may be different. 

MICROMETEORITES AS A TOOL 

Assuming that micrometeorites could be identified in returned soil 
samples, this addition of micrometeorite material to the uppermost Martian 
regolith at a constant rate could conceivably provide a powerful tool for 
tracking rates of erosion, deposition, and weathering. On Mars Sample 
Return missions, an attempt should be made to collect soils from different 
geologic sites (catch basins, lag surfaces, flat high plains, valley 
bottoms, etc.) so as to provide a variety of soils of different sedimentary 
environments. One of the important differences among these environments 
might be the proportion of petrographically or chemically identifiable 
micrometeorites mixed into the soil. 


CHEMICAL SIGNATURES 

The possibility that detectable micrometeorites and their remains can 
be found in the Martian soils depends on the relative rates of infall, 
weathering and alteration, transportation, and mixing. These rates are not 
yet known reliably enough to allow us to predict with certainty whether 
identifiable micrometeorites will be found. While they may be relatively 
quickly destroyed by Martian weathering, the chemical signatures, 
particularly siderophiles and volatiles, may persist in the soils, as they 
have in the lunar regolith (Anders, et al., 1973), and in Earth sediments 
where Ir anomalies are detected (Kyte and Wasson, 1986) . 

All of the meteoritic material collected by the planet, whether it reaches 
the surface unmelted, or if it melts or vaporizes on atmospheric entry, 
should eventually be added to the soils. We have estimated a mass influx on 
Mars ranging from 2,700 to 202,000 tons per year. If spread uniformly over 
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the 5 surface of the planet, this corresponds to an addition rate from 1.8 x 
10* 3 gm/nr’year to 1.4 x 10 gm/n/’year. For material of density 1 gm/cm , 
this corresponds to an accretion rate which ranges from 1.8 centimeters per 
billion years to 1.4 meters per billion years. 

If the soil production rate on Mars is of order 1 meter per billion 
years, these meteoritic accretion rates would give rise to meteoritic 
concentrations ranging from 2% to 58% in the average Martian soil, fioslough 
(1988) and Clark and Baird (1979), applying the subtraction method to the 
Viking chemical data, suggest the Mars regolith can be fit by a mixture of 
40% Cl meteorite and 60% planetary rock fragments. Boslough (1988) suggests 
the meteoritic component is ancient. But it could equally well be the 
micrometeorite component, which dominates the ancient component in lunar 
mare soils. In the lunar case the composition of the non- indigenous 
material was taken as the residual after subtracting rock composition from 
soil composition. Two meteoritic components were detected. In mature soils 
the residual has a trace element composition consistent with the addition 
of 1.5% Cl meteoritic material, attributed to the long term micrometeorite 
infall (Anders et al. , 1973). Less mature highland soils also show a second 
component, characterized by fractionated siderophile content and low 
volatiles, attributed to ancient bombardment (Anders et al., 1973). 

The observations of Boslough (1988) and Clark and Baird (1979) are 
consistent with the Mars soil at the Viking sites containing a substantial 
meteoritic component. The meteoritic concentration they infer is consistent 
with the range of concentrations we calculate from micrometeorite influx, 
provided the rate of production of soil on Mars is no more than a few 
meters per billion years. The possibility exists that the Martian soils 
contain a substantial fraction of meteoritic material. On Earth and on the 
Moon, the Ni/Fe ratio and the Ir abundance have proven to be diagnostic 
indicators of the meteoritic component, since Ir and Ni are enriched in Cl 
meteorites but depleted in crustal materials. Direct spacecraft measurement 
of the Ni and/or Ir abundances in the Mars regolith should help determine 
the meteoritic content of the soil. 

CONSTRAINTS FROM VIKING MEASUREMENTS 

The major meteoritic component in the Lunar mare soils is similar, in 
siderophile and volatile element composition, to the Cl carbonaceous 
chondrite meteorites. These meteorites contain an average of 3 to 5% carbon 
(Wasson, 1974), some in the form of organic matter. The cosmic dust 
particles collected from the Earth's stratosphere contain carbon at or 
above Cl concentrations (Blanford, et al., 1988), some of which may be in 
the form of polycyclic aromatic hydrocarbons (Allamandola, et al., 1987). 

If the Martian soil contained a substantial abundance of unmodified Cl 
carbonaceous chondrite material, the Viking gas chromatograph mass 
spectrometer should have detected its presence. 

Biemann, et al. (1977) found no detectable organic material in four 
Martian samples, one surface and one subsurface at each of the two Viking 
sites. Using the laboratory version of the Viking gas chromatorgaph mass 
spectrometer, they detect naphthalene at a level of about 1 ppm in Cl 
carbonaceous chondrite samples (Biemann, et al., 1977). They report a 
detection limit of 0.5 ppb for naphthalene at the Viking 1 site and 0.015 
ppb for Viking 2 (Biemann, et al., 1977). If the concentration of 
naphthalene were the same in the infalling meteoritic material as in their 
carbonaceous chondrite sample, then the corresponding upper limits on the 
meteoritic mass fraction in the analyzed Viking samples would be 0.05% at 
the Viking 1 site, and 0.0015% at the Viking 2 site. „ 

The organic content of the micrometeorites in the 10*° to 10 * z gram 
mass range has never been established, since the particles in this mass 
range collected on Earth are melted on atmospheric entry. Furthermore, as 
pointed out by Banin (1988) , the high redox potential of the Martian soil 
may have caused the decomposition of any organic matter from meteoritic 
infall. 
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MARTIAN AGGLUTINATES AND SOIL MATURITY 

If, as we calculate, micrometeorites are all slowed down by the Martian 
atmosphere, and assuming that most lunar agglutinates are made by 
micrometeorite impacts, no analogous Martian agglutinates would be expected 
(unless there were an era in which the atmosphere was considerably less 
dense than at present). However, many types of impact glasses would be 
expected from larger impacts , and some of these glasses may resemble lunar 
agglutinates in some respects. 

Gault and Baldwin (1970) have estimated a minimum impact crater size of 
50 meters, taking into account fragmentation and ablation of the incoming 
projectiles as well as atmospheric deceleration. The smallest craters 
noted in Viking orbiter images are about 100 meters in diameter (Blasius, 
1976), but smaller craters beyond the resolution limit of the photographs 
may still be present. Dycus (1969) predicts that projectiles as small as 
10 gm would still form craters. However, craters too small to be seen from 
the orbiter are not apparent in Viking lander images. Impact gardening 
associated with the 50 meter and larger craters predicted by Gault and 
Baldwin (1970) would determine regolith turnover rates and cause 
comminution of rocks into soils. The addition of micrometeorites would 
affect the petrology and chemistry of Martian soil. Weathering and 
sedimentary processes on Mars would also process the regolith components . 
The overall effect would be to make an exceedingly complex regolith. A new 
maturation scale will be necessary for Martian regolith. This scale will 
have to include terms which reflect (1) impact reworking, (2) addition of 
micrometeorites, and (3) Martian surface weathering and alteration. For 
example, if concentration mechanisms can be factored out, the abundance of 
micrometeorites (identified petrographically or chemically) in a soil layer 
might be directly related to its near -surface exposure time in a manner 
analogous to the abundance of agglutinates in lunar soils. In addition to 
soil evolution through maturation, physical mixing of soils of differing 
maturities should be common. 


CONCLUSIONS 
7 ^ 

Micrometeorites in the 10 to 10 gram mass range should be a major 
contributor to the meteoritic input on Mars. Unlike the Earth, where most 
particles above 10” b grams are melted or vaporized on atmospheric entry, a 
large fraction of these particles are expected to survive entry into the 
atmosphere of Mars unmelted. The addition of this meteoritic material to 
the Martian regolith could significantly perturb the chemical abundances in 
the soils, particularly the abundances of volatile and siderophile elements 
which are abundant in Cl meteorites but depleted in crustal materials, and 
of noble gases (and possibly hydrogen) which are implanted in 
micrometeorites during space exposure and carried into the soils with the 
particles . 

Uncertainties micrometeorite flux at Mars as well as the rate of 
production of soil through weathering processes on the planet give rise to 
large uncertainties in the meteoritic concentration in the Martian soils. 
However, our estimates are not incompatible with the suggestion by Boslough 
(1988) that the soils analyzed by Viking are 40% meteoritic. 

The first returned soil samples from Mars should provide the 
opportunity for recovery and analysis of unaltered micrometeorites larger 
than any sampled on earth, assessment of the magnitude of the meteoritic 
component, and and possibly an estimate of the rate of soil production on 
the planet. The larger micrometeorites which enter the atmosphere of Mars 
without melting may sample a different source population than is sampled by 
the smaller particles collected from the Earth's stratosphere. 
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Table I 


Mass Influx at Earth 


Mass Range 

Particle Diameter 
(microns ) + 

Mass Influx 

(grams) 

(kg/year) 

10*7 tQ i 0 *6 
10*® to 10* 5 

58 to 124 

2.0 x lof 

124 to 268 

4.2 x 10® 

10'? to 10*? 

268 to 576 

4.2 x 10® 

10* 4 to 10*;? 

576 to 1240 

2.7 x 10® 

10* J to I0' z 

1240 to 2600 

1.3 x 10 b 


Fraction of 
Total Mass* 

.12 

.26 

.26 

.17 

.08 


* Data, except for particle diameters, from Hughes (1978). , 

# Ratio of mass in this decade to total mass influx from 10* 1,5 to 10° 

grams . , 

+ Diameters calculated for spheres of density 1 gm/cm . 


Table II 

Particles Surviving Mars Atmospheric Entry Unmelted 


Size Range 

Upper Estimate 

Lower Estimate 

(micrometers) 

(particles/nr "year) 

(particles/m z * year) 

56 to 124 

1.6 x 10? 

2.1 x 10 1 

125 to 268 

3.2 x lOf 

4.3 

269 to 576 

2.7 x 10 1 

3.6 x 10 ? 

577 to 1240 

1.2 , 

1.6 x 10 m f 

1241 to 2600 

3.1 x 10* z 

4.1 x 10' 4 



Table III 


Number of Unmelted Micrometeorites Expected 


in 10 gram Average 

Soil Sample 

Size Range 
(micrometers) 

Upper Estimate 

Lower Estimate 

56 to 124 

16,000 

210 

125 to 268 

3,200 

43 

269 to 576 

270 

4 

577 to 1240 

12 

< 1 

1241 to 2600 

3 

< 1 
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ABSTRACT 


To assure the quality of potable water on the Space Station a number of 
chemical and physical tests must be conducted routinely. After reviewing 
the requirements for potable water, both direct and indirect analytical 
methods are evaluated that could make the required tests and inprovements 
are suggested be make them compatible with the Space Station operation. 

Of the possible analytical procedures that might be carried out on the 
SS, the following seem to be required to meet the water quality 
criteria: turbidity and color from photometric measurements, a conduc- 
tivity measurement, pH by an ion selective electrode, ion chromatography 
for inorganic ions with set MCIs, and a total organic carbon measurement. 
Also required are generic methods for total phenols (MCL 1 ppb) , total 
organic halogen compounds (MCL 10 ppb) , and organic acids (MCL 500 ppb) . 
Electrochemical methods of controlled potential amperametry and micro- 
coul ©metric titrimetry may meet this need. The other FW standards for 
organic alcohols (MCL 500 ppb) and TOC less nontoxicants (MCL 100 ppb) , 
are broadly defined and can only be met with either an approximate 
surrogate measure (example, UV absorption of aromatic compounds for 
"toxic TOC") , or methods that separate and measure specific compounds in 
those classes that are likely to be found in PW. Finally the problem of 
analyzing for specific toxicants that are yet to be defined will either 
require specific sensors (electrode or optrode) , or a separation and 
detection system. 

Each of the separation techniques evaluated appear to have special areas 
of applications for the SS, though supercritical fluid chromatography 
with CO2 could probably replace gas chromatography while extending the 
range of analyses to nonvolatile heat sensitive compounds. Most forms of 
HPLC are not very compatible with the SS environment, but ion chroma- 
tography is needed for most anion and cation analyses (anodic stripping 
voltammetry could measure the heavy metal ions) . It would be more 
suitable if it were converted to microbore column size. Capillary zone 
electrophoresis has the potential as a universal separation system, but 
it will require much more development in terms of buffer composition and 
sensitive micro volume detectors. 

The most universal detector for the separation techniques is the high 
resolution mass spectrum detector, however its use on the SS is limited 
because of size, power and weight requirements associated with its high 
vacuum system. The need to have one on the SS seems justified only if it 
would be used for a variety of experiments and not just to test FW. The 
electrochemistry/conductivity detector appears to be the best compromise 
for a detector for HPLC or CZE. 

Techniques are suggested for removing the interfering biocide I2/I” from 
water, and for carrying out chemical and electrochemical operations in 
microgravity where mixing and gas liquid separations are required. 
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SPACE STATION POTABLE WATER 


The potable water (FW) for use on the Space Station (SS) will come from 
the cabin air humidity condensate (HC) produced by the air revitalization 
system (ARS) . Ibis means the EW is in a direct recycle loop, a treatment 
process not previously implemented. To insure the EW remains safe for 
the duration of a 90 to 180 day mission, it needs to be continually 
monitored. This includes frequent or continuous testing of several 
important general or surrogate indicators that show the treatment system 
is operating properly. There are also a number of other specifications 
that the FW must continue to meet with regard to specific chemicals or 
classes of compounds. Table 1 gives the current criteria set for SS EW 
(1) . The table does not include the microbiological standards, as they 
are not of concern to this report. 

Unfortunately the majority of the organic carbon found in samples of HC 
from the two experiments that simulated the expected SS conditions (Space 
Lab and Bends) have not been correlated with specific organic compounds 
(2,3). Therefore, maximum contaminant levels (MCL) for possible toxic 
organic chemicals in EW can not be set until more specific information is 
obtained about the composition of the TOC in representative samples of HC 
from ground based tests. This, by itself, is a major analytical 
challenge and can influence the types of chemical analyses that are 
required on the SS. 

This report will attempt to evaluate the kinds of on-board analytical 
procedure that will be required and also look at alternate methods that 
accomplish the same ends and/or add valuable data for determining EW 
quality. Modifications of procedures needed to meet the required 
sensitivity and potential interferences will be suggested along with 
ideas for efficient operation within the limitations of the SS environ- 
ment. First the direct analytical methods which allow continuous on-line 
monitoring of the EW will be discussed, then the indirect methods which 
require discrete samples will be considered. These later methods require 
chemical reaction or separation and are the more sophisticated and 
complex analytical techniques. 

SPACE STATION OPERATIONAL LIMITATIONS 

There are four rather severe limitations on the analytical chemistry 
procedures that can be used on the SS. 

1. The biggest limitation is microgravity which does not allow the 
automatic phase separations of immiscible fluids. That is, gases and 
liquids do not separate without an external applied force. They can be 
separated by diffusion, but that is slow and inefficient. However, in 
zero G convection does not occur due to small density differences, this 
can be helpful for processes like electrochemistry or capillary zone 
electrophoresis (CZE) that will be discussed later. 
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TABLE 1.- SPACE STATION POTABIE WATER QUALITY REQUIREMENTS, MAXIMUM 
CONTAMINANT IEVELS (1) . 


Ehysical Parameter Limits 


Total Solids 

(mg/1) 

100 

pH 

6. 0-8.0 

Conductivity 

(Mmhos/cm) 

- 

Turbidity (NIU) 

1 

Color, true 

(Pt/Co units) 

15 

Dissolved Gas 
(free at 35 C) 

None 

Taste & Odor 

(TTN/TON) 

3 

Free Gas (§ STP) 

None 


Inorganics mg/1 (ppn) 


Ammonia 

0.5 

Cyanide 

0.2 

Mercury 

0.002 

Arsenic 

0.01 

Fluoride 

1.0 

Nickel 

0.05 

Barium 

1.0 

Iodine 


Nitrate 

10. 

Cadmium 

0.01 

Total 

15. 

Potassium 

340. 

Calcium 

30. 

Iron 

0.3 

Selenium 

0.01 

Chloride 

250. 

Lead 

0.05 

Silver 

0.05 

Chromium 

0.05 

Manganese 

0.05 

Sulfide 

0.05 

Copper 

1.0 

Magnesium 

50. 

Sulfate 

250. 


Zinc 5.0 


Organics Mg/1 (ppb) 

Aesthetics 

(mg/1) 

** Total Organic Carbon (TOC) 

500. 

Cations 

30 

TOC (Less Nontoxicants) 

100. 

Anions 

30 

Generic Classes 


002 

15 

Acids 

500. 



Alcohols 

500. 



Halogenated Hydrocarbons 

10. 



Fhenols 

1 . 



** Specific Toxicants 

TBD 




** Note, less than half the TOC from Space Lab humidity condensate 
samples have been associated with specific organic compounds. 


2. Mission personnel will not be specifically trained in chemistry and 
their time will be limited for chemical or instrumental manipulations. 
This requires careful consideration of analytical techniques in terms of 
in-flight serviceability, frequency of instrument reconfiguration, 
recalibration or component failure, and consumable servicing. 

3. Due to health and weight considerations, limits will be set on both 
the amount and toxic nature of consumable chemicals used in analytical 
procedures or instruments. There is special concern with toxic volatile 
waste chemicals that could get into the cabin atmosphere. 
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4. The SS will have limited amounts of power and space available for any 
operation, so these factors as well as weight will influence the choice 
of analytical instruments. Modularity and interchangeability of 
components between different instruments should be a goal in order to 
reduce the inventory of spare parts and improve serviceability. 

DIRECT METHODS 

Direct methods are able to continuously monitor the FW stream without 
personnel involvement. They have a specific response to the parameter of 
interest and should not respond to other substance. The methods we will 
consider are conductivity, photometry, and electrochemistry. These 
methods are roost useful for monitoring general or surrogate parameters 
relative to the function of the FW purification system, though some are 
applicable to specific substances. 

Conductivity 

This is the basic test for total ion concentration in the water and 
increasing conductivity indicates the ion exchange resins need renewing, 
or there has been a malfunction and excess ions are entering the FW 
system from some source. 

Usual equipment: Two platinum electrodes (with Pt black coating) with a 
fixed area and spacing, energized with a constant small amplitude sine 
wave voltage and read with a Wheatstone bridge circuit. The output is 
either in conductivity (units pS/ cm, previously /imhos/cm) or its 
reciprocal, as resistivity (ohms cm or megohm cm) . 

Suggested improvements for SS operation: The commercial equipment can be 
used "as is" with the only interference coming from gas bubbles in water. 
However, the system can be improved by using a bipolar pulse detector 
circuit (4) which allows simple rugged electrodes of materials like 
stainless steel to be built into the water line (example, ring electrodes 
separated with Teflon insulating rings) . This method is very sensitive 
and responds rapidly to conductivity changes. It is used in many modem 
instruments, often with micro cells. The electronics are simple and 
interface readily to computer control. In principle, one bipolar pulse 
circuit could be multiplexed to all the conductivity detectors and other 
resistance devices (like thermistors) on-board. If several units are 
used they would have the advantage of interchangeability. 

Photometry 

This is a group of methods including turbidimetry, spectrophotometry, and 
fiber optic sensors (optrodes) , all based upon sensing the amount of 
light that is absorbed or scattered at one or more wavelengths. 

Turbidimetry is a required water quality measurement related to the 
amount of light scattered from particulate matter suspended in water. If 
the scattered light is actually measured the method is called 
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nephelometry and is more sensitive for low light scattering. On the SS, 
increasing turbidity will generally indicate the presence of bacterial 
growth in the system or breakdown of the activated carbon or ion exchange 
resins. Detergent molecules, which form micelles, and gas bubbles will 
also scatter li^it. This permits the detection of unallowed gas bubbles 
in the FW as their size, and thus turbidity, is pressure sensitive. 

Spectrophotometry is required for the FW "color" requirement. It is also 
being considered as a surrogate measurement for total organic carbon 
(TOC) in FW. However, individual organic compounds vary considerably in 
the amount of light absorbed at a specific wavelength, with roost 
aliphatic compounds only weakly absorptive (e < 40) above 200 nm. The 
method, if carried out at several specific wavelengths, may give a 
reasonable measure of aromatic compounds, as they strongly absorb (e > 
8000) in the region of 200 to 350 nm. If one makes the assumption that 
aromatic containing chemicals are the most hazardous, this measurement 
would be closer to a surrogate response for the "TOC less nontoxicants". 

Usual equipment: Fhotametric processes require either laser or broad 
band, highly stable, light sources (with monochrcmeters or filters to 
select the desired wavelengths) and an optical train free from unwanted 
stray light and vibrations. The detectors need to be very sensitive with 
a wide dynamic range of response. This means they need very well 
regulated power supplies, often with either high voltage or current 
output capability. 

Suggested improvements for SS operation: If several optical measurements 
are to be made, it may be possible to combine them, using only one or two 
light sources (UV and Visible) and one monochromator, with fiber optics 
used to pick up specific wavelengths and distribute them to the saitple 
cells. Also, absorption measurements for contaminates of FW requires 10 
cm or longer path length cells. For a typical aromatic compound at 100 
ppb, with a molar absorbtivity of 8xl0 4 , a 10 cm cell would give only an 
absorbance of 0.01 units (only 2.3% of the li#it absorbed). Other 
problems for photometric measurements on the SS come from absorbance due 
to added biocide (forms of iyi") / this may need prior removal (see 
section on "Silver Bullet") , and need to correct for scattering due to 
gas bubbles or turbidity. Finally to circumvent these problems the use 
of "cptrodes", optical sensors equivalent to electrodes, may be used for 
some specific determinations. Cptrodes have a reagent immobilized on the 
end of a double fiber optic that gives a color or fluorescence with the 
species of interest. Light is sent down one fiber and returned up the 
other and does not enter the solution (5) . 

Electrochemistry 

Direct measurements with electrochemistry include potent iometry , using 
ion and molecular selective electrodes (ISE and MSE) and oxidation 
reduction potential (ORP) electrodes, and some applications of con- 
trolled potential airpercmetry (CPA) . 



An XSE for pH determination is required for I W, and the activity of 
several other specific ions could be measured, examples I and F - , 
however many other ISE suffer interference problems. MSEs could be used 
to measure 0 2 , and perhaps CO 2 and NH 3 . They also have the possibility 
for determining biochemical molecules through immobilized enzymes ( 6 ) . 
The problem with most of these ISE and MSE is their need for frequent 
recalibrations. While not required, an ORP electrode (made of Pt or Au) 
would be useful for monitoring the redox potential of EW, which would 
indicate the level of I 2 and/or 0 2 . A more negative trend in the 
potential may indicate anaerobic bacterial grcwth in the system. 

The use of CPA with large area inert electrodes could be valuable for 
monitoring the concentrations of several classes of organic compounds 
present in PW, Table 2. In aqueous solutions the range of potentials 
available for CPA, before H 2 or 0 2 formation occurs, is about + 1.0 to 
-1.0 volts. The method might be used as a surrogate method for phenols 
(MCL of 1 ppb) and halogenated hydrocarbons (MCL of 10 ppb) . with porous 
electrodes CPA could be used to remove the interfering 1 2/1 from EW 
before the trace organics are measured and to regenerate I 2 from I - . 


TABLE 2.- POTE N T I ALLY DETECTABLE ORGANIC FUNCTIONAL GROUPS BY CONTROLLED 
POTENTIAL AMPEROMETRY (7) . 


Oxidizable Groups 

Azines 

Hydrocarbons 

Amides 

Amines 

Quinolines 

Ehenols 


Potential 
(V vs SCE) 

> + 1.2 
> + 1.0 
> +0.5 
+1.0 > +0.5 
+ 0.6 > + 0.2 
+0.5 > 0.0 


Reducible Groups 

Halogens 
Nitro conpds 
Diazo ccnpds 
Ethers 
Esters 
Aldehydes 
Ketones 
Olefins 


Potential 
(V vs SCE) 
+ 0.2 > - 0.2 
-0.2 > -0.5 
- 0.2 > - 0.6 
-0.6 > -1.4 
- 0.8 > 

- 1.1 > 

-1.4 > 

-1.7 > 


(Note: This list represents general trends, as it is by no means complete 
nor totally accurate based upon my experience. RDG) 


Usual equipment: For potenticmetric measurements an indicator electrode 
and a reference electrode are needed, one reference electrode would work 
for several indicator electrodes if they are in close proximity. A very 
high impedance potentiometer, equivalent to a pH meter, is needed to 
sense the voltage without disturbing the measurement. For controlled 
potential amperometry usually separate reference, counter and working 
elec t rodes are needed, the latter two would be made of either Pt or Au. 
A three electrode potentiostat is usually needed to control and monitor 
the potential and current. 

Suggested improvements for SS operation: Commercial equipment can be 
used, one precaution is the need for gel filled reference electrodes to 
avoid bubble formation between the electrode and the salt bridge 
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junction. All the indicator electrodes could be multiplexed thru one 
output console if each electrode has its own very high impedance 
operational amplifier. The problems to overcome in potenticmetry are the 
stability of the electrodes and a means of cleaning and recalibrate them 
during a mission. If the electrodes are used in-line with the FW, a 
bypass loop to one or more special calibration solutions might be 
required. The CPA instrumentation is simple and should be interface thru 
a computer for digital control of potentials and currents. The 
electrodes may be subject to long term fouling if bacteria are present 
and may need periodic cleaning and exchange between missions. 

INDIRECT METHODS 

Indirect methods usually involve periodic sampling with chemical 
treatment, or sample separation followed by determination of the species 
of interest. The reaction processes to be discussed are TOC and two 
electrochemical techniques, anodic stripping voltammetry and micro- 
coulcmetric titrimetry (an insitu titration method) . Also disscussed are 
a semi-micro, zero G reaction system and a method for removing the 
interference of I^I". The separations methods of gas chromatography 
(GC) , supercritical fluid chromatography (SPC) , hic£i performance liquid 
chromatography (HPLC) , and CZE will be discussed briefly along with their 
detectors and required interfaces suitable for microgravity operation. 

Chemical reactions 

The TOC method oxidizes the organic carbon present in a FW sample and 
measures the C0 2 released. This is an important surrogate measure of the 
organic quality of FW and should be carried out routinely on the SS. 

Usual equipment: Organic matter in water samples (<10 ml) is oxidized in 
a reactor either by photochemical or thermal decorposition of persulfate, 
or by seme other oxidation catalyst. The 00 2 produced is measured by 
conductivity or IR, usually after traping the 00 2 in molecular sieves and 
purging it to the detector. The methods have problems with incomplete 
oxidation of refractory organic molecules. High temperature persulfate 
oxidation appears to give superior results and will be the method 
discussed. Inorganic carbon, H00 3 ", interferes and is removed by prior 
N 2 purging of the acidified sample. This causes seme volatile organics 
to be removed and they must be trapped and oxidized to 00 2 . These 
separate mesurements of 00 2 are called TIC and TVOC respectively. 

Suggested Improvements for SS operation: Since OC^ is sparged from the 
reactor with N 2 or air, the commercial method needs adaptation to work in 
microgravity. The method also consumes sodium persulfate and phosphoric 
acid solutions that should be considered hazardous in the SS environment. 
An "orbital reactor" (see indirect electrochemical methods) could be used 
to handle the gas-liquid separation problem and to conduct the 
persulfate high temperature oxidation. The acid and persulfate could be 
generated by standard electrochemically reactions within the reactor to 
avoid the need for consumables other than a sodium sulfate solution. 
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Microcoulometric titrimetry of the C0 2 collected by purge and trap 
methods would allow for smaller water samples, down to 0.1 ml, than the 
current IR detector system. The reaction could be carried out in a 0.5 
ml cell that would require less heat and smaller volume of consumables. 
These changes effectively convert the TOC determination to an electro- 
chemical method. 

Elecfa T*-»tt*inrs*i techniques often require the removal of interferences due 
to 1 2 / 1 ” and 0 2 at the concentrations expected in FW, as they give 
appreciable reduction currents at small electrode potentials. I would 
suggest a pair of porous silver/silver iodide electrodes separated with a 
Nafion ion exchange membrane that can reduce iodide to very low levels, 
base upon the applied potential and the Agl solubility product of 
1 . 5x1 6 at 25 °C. This "silver bullet" (SB) works on the principle 
that as the required water passes thru one electrode, 1^ will react spon- 
taneously with Ag to form Agl and I" and further anodic oxidation of Ag 
will form more Agl on the electrode. Meanwhile, as a small amount of 
water passes thru the other electrode, its coating of Agl is reduced to 
Ag and T~ is released (this concentrated iodide solution can be collected 
for re-use in the FW loop after oxidation to I^I") . The relative amount 
of I" in the two water streams is controlled by the applied voltage, 
with 350 mV giving a ratio of 10 6 . 

Note, the SB method can be readily scaled up to remove 1“ from FW for end 
uses. W at 1 liter per minute would need a maximum current of 190 mA 
(27 mA nominal, 3 ppm I t otal) . In this use the SB would have three 
connection positions, then as the Agl builds up on one electrode and is 
depleted on the other they can be interchanged so the opposite water 
streams flow thru the electrodes and their functions are exchanged. The 
third position would allow bypassing the electrodes to periodically 
disinfect the terminal portion of the water system with biocide. 

Removal of 0 2 from water can also be accomplished by electrochemical 
re du ction in similar flow cells with Pt electrodes, or by sparging the 
solution with pure nitrogen or argon. An "orbital reactor" (OR) is 
suggested as a microgravity method for the mixing and sparging of 
solutions required for chemical or electrochemical reactions. The OR is 
a platform holding a reaction vessel that is given a orbital motion while 
its orientation remains fixed, the only moving parts associated with the 
reactor are the connecting wire leads and tubing. These need to flex 
with the motion of the OR and their position and composition must be 
selected to reduce stress fatigue. The reaction cells can be of 
numerous sizes and shapes, but for electrochemistry a disk or washer 
shaped cell of 0.5 ml to 5 ml should work well. For example, a washer 
shaped cell 0.25 mm thick with radii of 1 and 2 cm, the volume would be 
0.5 ml and with an orbiting rate of 212 rpm the force would be 0.5 and 1 
G at the inner and outer surfaces respectively. For purging the OR cell, 
gas would be introduced at the periphery of the cell and it would move to 
the center in a spiral motion. The shear of moving water by the gas 
inlet would cause the bubbles to be small and well separated, giving very 
efficient sparging. The optimum amount of G force needs to be 
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determined, as little as 0.1 to 0.2 G may be sufficient for the these 
operations. The flat shape and large area of these OR cells allows the 
placement of a variety of sensors in the same cell, for example, conduc- 
tivity, large area, micro, reference, counter, and thin layer electrodes 
of several different materials, as well as optrodes and thermistors. 
Divided electrochemical cells for microcoulometric titrations, where the 
reagent is generated in situ, can be made by splitting the OR reactor 
horizontally with an ion exchange membrane. 

Anodic stripping volt a maetry is a common method for measuring trace (to 
sub ppb) levels of the heavy metal ions including Od, Cu, Pb, Hg, Ni, Ag, 
Fe, Zn, and possibly Cr, Mh, and As in water. This method would replace 
most of the water analyses done by ground based graphite furnace atomic 
absorption. This method concentrates the metals by plating them out on a 
inert electrode, held at a sufficiently negative potential, from an 
uniformly stirred solution for a specific period of time. After this 
time the stirring is stopped and when the solution cones to rest the 
metals are removed by a positive voltammetric sweep. The amount of 
charge required to remove each metal at its specific redox potential 
corresponds to the metals concentration in solution. Often the metals 
are plated into a co-deposited mercury film and the procedure may require 
the addition of supporting electrolyte. 

hhcrocoulanetric titrimetry of the acid and base content of EW can be 
carried out in a divided OR cell which allows the production of either H+ 
or OT, in the FW sample, depending upon the electrode polarity. By 
applying constant current pulses of short duration, specific amounts of 
acid or base are generated, and with the use of a sensitive bipolar pulse 
conductivity detector, the total acid and base content of the EW can be 
determined frcm titration curves. The detection limits should, in 
theory, be at the sub ppb level for acid and base compounds, as a 1 /xA 
pulse for 100 msec is equivalent to 0.1 microcoulcmb or 0.2 ppb for a 100 
molecular weight acid or base in 0.5 ml. This technique might be 
developed to give a surrogate measure of organic acids, bases and 
possibly phenols in EW. It might also be used for measuring 00 2 and NH 3 
in the EW and overcoming their analytical interference. As mentioned 
earlier the technique could be adapted to do the TOC analysis. The 
method could also be used in other SS chemical analyses where a needed 
titrarrt can be generated microcoulcmetrically. 

Separation Techniques 

There are four general types of high resolution separation techniques 
that may be considered for identifying and quantitating components 
(mainly organic) of FW and other samples on the SS. Three are based upon 
pressurized fluid elution: GC, generally using nitrogen gas, SFC, using 
C0 2 above its critical pressure (>1100 psi) , and several variations of 
HPLC, using aqueous/organic solvent mixtures or aqueous buffers (to 
pressures of 6000 psi) . The other one, CZE, uses very high electric 
field driven elution and ccmes in several variations. For analyzing 
trace organics in EW these techniques will probably require concentration 
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of samples using some arrangement of purge and trap. 

The GC technique is excellent for separating volatile heat stable organic 
compounds and interfaces well to a variety of detectors. It is 
considered the main tool for cabin atmosphere monitoring and could be 
useful for scare specific FW analyses. 

Usual equipment: The method requires a source of high pressure gas (100 
psi) a sample injection system (sampling loop and valve) , a temperature 
programmable oven (to 350 °C) , a long silica capillary column with a 
bonded stationary phase (several may be needed) and an interface to a 
detector. 

Suggested improvements for SS operation: Most commercial equipment would 
be compatible with the limitations of the SS environment. The small N 2 
requirement can be met by using cabin make-up gas if properly purified, 
or generated electrochemical ly from hydrazine sulfate (toxic) . The 
standard flame ionization detector (FID) should be readily interfaced to 
the column if a premix H 2 /O 2 micro burner is used to avoid microgravity 
effects. Otherwise a photot ionization detector (PID) might be adaquate. 
A fourier transform infrared detector (FI'IRD) should work directly, but 
the electrolytic conductivity detector (E1CD) that is excellent for 
halogen compounds would need modification of its gas/liquid reactor. The 
mass spectra detector (MSD) , which will be discussed later, would 
interface easily. 

The SFC technique is still developing and currently is able to separate 
heat sensitive compounds which are non-volatile, if they are not too 
polar, as well as the compounds handled by capillary GC. The technique 
falls between GC and HPLC in its applications and interfaces to most of 
the detectors used for GC. 

Usual equipment: The equipment is usually a modification to either GC or 
HPLC instruments, and requires high pressure pumps (often the syringe 
type) , usually a column oven for temperature up to 150 °C and an inter- 
face to a detector. 

Suggested improvements for SS operation: Of the three chromatographic 
techniques this is the most compatible with SS operation since it uses 
00 2 (which is recycled on the SS and is non-hazardous ) and has modest 
power requirements. An advantage of supercritical C0 2 is its ability to 
purge traps of all but their most polar organic compounds. The eluting 
ability of 00 2 can be modified by changing its density (pressure) and by 
saturating it with water or polar organics. SFC has about the same 
detector limitations as GC. 

The HPLC techniques have reached a high degree of sophistication and the 
ones of interest are normal and reverse phase partitioning for organics 
(NFP & RPP) , ion pairing for ionic and ionizable organics (IP) , and ion 
exchange, mainly for inorganic ions, called ion chromatography (IC) . In 
combination these techniques have the ability to separate most mixtures 
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of organic compounds including the most polar and non-volatile, though 
the resolution is not as high as with GC or SFC. 

Usual equipment: The methods all require special very high pressure 
pumps and controllers to push eluent solutions rapidly thru short (10 cm) 
small bore columns packed with special very small diameter particles 
designed or coated for the particular application. This means that 
either dedicated instruments are needed for each application or columns 
and eluent systems must be changed for different analyses. For each 
sample the NPP and RFP systems use about 10 to 15 ml of a complex eluent 
mixture, often containing acetonitrile or methanol as a major component. 
The IC requires several aqueous ionic buffer combinations and ion 
exchange columns to separate the different sets of ions effectively with 
a sub ppb level of sensitivity. The IP technique is a combination of the 
previous two and requires an organic solvent with ionic buffer and 
separation of organic ion pairs is usually on RPP columns. The detectors 
for these methods are generally ones that operate well with solutions, 
for example, conductivity, electrochemical (often CPA) , and various 
photometric methods. Recently, by means of a complex interfaces, several 
types of mass spectra detectors (MSD) have been successfully used with 
HPLC (see the section on detectors) . 

Suggested improvements for SS operation: The large amount of consumable 
buffers and their waste, some of which are toxic, may limit the use of 
HPLC in the SS. The amount of eluent can be reduced about 80 fold by 
using "microbore" columns (8) . However, the IC which is the most 
important for FW analysis apparently does not have the extra small 
particle ion exchange resins needed for microbore modification at this 
time. 

The CZE techniques are very new and hold much promise as a universal 
separation technique (9) . They use the principle of electroosmotic force 
to move buffer solutions down small bore (<100 jm ) while at the same 
time vising the electrophoretic mobility to separate charged species. 
Besides zone electrophoresis itself, there are two variant techniques, 
isotacophoresis that uses a discontinuous ionic buffer to focus and 
separate analyte ions, and electrokinetic chromatography that uses an 
ionic micelle buffer to separate neutral organic molecules. 

Usual equipment: No commercial equipment is currently available, but the 
apparatus is simple to build. The methods requires a 30 kV DC at 10 juA 
power supply, a 0.5 to 2 m lenght of fused silica capillary tubing, a 
very small volume of aqueous buffers, and an interface to a detector. 
Sample size for CZE is about 10 nl, and the resolution is about equi- 
valent to capillary GC. Because of the small sample size and column 
diameter it requires very sensitive and high resolution detectors. 

Suggested improvements for SS operation: The techniques are very 
compatible with the SS environment, with small weight, volume, power and 
consumables (non-hazardous) requirements. However, considerable 
development is needed to find the right buffer protocols to accomplish 
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the required separations. For convenience of SS operation, several 
columns and buffer setups may be required that can be switched between 
one power supply. It may be necessary to change capillary columns 
frequently due to plugging or other failure and replacements with easy 
means of installation is needed. Samples of EW will have to be con- 
centrated for use with these techniques. For the required sub m 1 sized 
samples, this might be accomplished without loss of too many volatile 
organics by diffusion of water out thru a membrane. 

Detectors 

A number of detectors are available for the separation methods just 
discussed, and some were commented upon in that section. The mass 
spectra detectror (MSD) in its several variations is close to the 
universal detector, giving nearly absolute identification of unknown 
components, requiring very small samples, having very high sensitivity, 
and interfacing easily to GC and SFC. The interface to HPLC and CZE is 
more difficult because the the eluent is a liquid. Several techniques 
have been developed, but there are problems with getting all analyte 
species ionized and into the MSD with equal efficiency. Electrochemical 
and conductivity detectors interface readily to the HPLC and CZE columns 
and are very sensitive and can be made with very small volume cells. 
There are several other detectors which may be useful but will not be 
discussed: photoionization, UV absorption, fluorescence and infrared. 

The MSD methods include the low resolution units (Ion Trap, Mass 
Sensitive Detector and related ion drift methods) that give only the 
parent mass/charge ion of unit resolution, and the high resolution mass 
spectra systems (quadripole units are the smallest and simplest) that 
give a complete molecular formula and often an unambiguous structure. 
Except for the ion drift method, the detectors all require a very good 
quality vacuum (10 -7 Torr) and need high velocity pumping for inter- 
facing to HPLC or CZE. The requirements of complexity, weight, volume, 
and energy consumption make the MSD the least compatible of the analyt- 
ical techniques for SS, however their advantage may be worth the effort 
to develop a workable system. The most likely candidates would be an ion 
drift, or ion trap detector. The ion trap is the smallest detector, 
excluding the turbamolecular pump, it is about 64 cubic inches, while 
the ion drift unit works with a plasma at atmospheric pressure. 

The MSD interface requires removing the volatile part of the eluent 
stream from the analyte species and converting them to "dry" ions while 
introducing them into the mass/charge detector system. Two interfaces 
currently in use are electrospray that uses high voltage to nebulize and 
ionize the droplets and theromospray that uses rapid heating of the 
effluent containing special buffers (usually NH 4 HCO 3 ) to produce ioniza- 
tion under these conditions. The removal of solvent is aided by a rapid 
flow of dry N 2 in front of the MSD inlet. There is still a need for 
research on these interface methods and a combination of the two may be 
needed for the special conditions on the SS. 
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Electrochemical and cxxiductivity detectors can both use the same 
microcell and electrode design to interface efficiently to HPLC and CZE 
systems. The conductivity detector would use the same bipolar pulse 
circuit that has been indicated for other conductivity measurements. 
The electrochemical detector would normally operate as a controlled 
potential amperemeter (CPA) , but could be used for rapid scan cyclic 
voltammetry for unknown identification. The micro volume cell that would 
fit directly on the end of either microbore HPLC or CZE columns could be 
constructed in sandwich fashion from thin Pt or An foils (0.01 mm) 
between thin Teflon films (0.03 mm) . A hole to match the inside diameter 
of the column (0.1 mm or less) would be punched or laser drilled through 
the layers to make a set of thin ring electrodes. The volume of the cell 
would be about 0.3 nl and would not degrade the resolution of a CZE 
column. Except for identifying unknowns, where accurate mass and 
structural information is needed, these detectors would work very well 
for most HPLC or CZE separation on the SS. 

RECOMMENDATIONS 

A variety of suggestions have been made throughout this report to improve 
the analytical techniques for SS operation. The most important 
recommendations are: 

(1) The silver/silver chloride electrode (SB) method of removing Io/I" 
biocide from the water, since it may interfere with analytical 
procedures for FW and also its end uses (page 13-10) . 

(2) The orbital reactor (OR) method of carrying out chemistry and 
electrochemistry in microgravity by using a disk shaped reactor on an 
orbital table to impart artificial G force to the contents, allowing 
solution mixing and separation of gases and liquids (page 13-10) . 

(3) A simple ultra low volume hi^ily sensitive electrochemical/- 
conductivity detector for use with a capillary zone electrophoresis 
apparatus (page 13-15) . 

It is also recommended, since several different conductivity and resis- 
tance measurements are made during the analysis of EW, that the bipolar 
pulse measuring circuit be used in all these applications for maximum 
compatibility and redundancy of equipment (p>age 13-5) . 
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ABSTRACT 


The Manipulator Emulator Testbed (MET) is to provide a facility 
capable of hosting the simulation of various manipulator configurations to 
support concept studies, evaluation, and other engineering development 
activities. Specifically, the testbed is intended to support development 
of the Space Station Remote Manipulator System (SSRMS) and related sys- 
tems. The MET is required to permit that components simulated in software 
may be replaced by future hardware components. 

The objective of this study is to evaluate the math models developed 
for the MET simulation of a manipulator's rigid body dynamics and the 
servo systems for each of the driven manipulator joints. Specifically, 
the math models are examined with regard to their amenability to pipeline 
and parallel processing processing. Based on this evaluation and the 
project objectives, a set of prioritized recommendations are offered for 
future work. 
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INTRODUCTION 


The Manipulator Emulator Testbed (MET) is to provide a facility in 
which different manipulator configurations may be simulated. The MET will 
be used as a tool to support Space Station manipulator development. It 
will be used to develop and evaluate concepts, support design and develop- 
ment studies, and evaluate hardware components and flight software 
modules. The MET is required to be designed such that initial software 
simulated components may be replaced by future hardware elements [1] - 
[3]. 


The MET is currently built around a network of IBM PC-AT computers. 
Each computer operates at an 8-MHz clock rate and is equipped with an 
Intel 80287 math co-processor, 640K bytes of memory, an extended graphics 
adapter (EGA) card, a color monitor, a 30-MB hard disk drive, and a 
high-capacity floppy disk drive. Each PC is equipped with a National 
Instruments GPIB-PCAA IEEE-488 interface card. The network includes an 
operator control station, data recording and display capability, and 
hardcopy output. A Multibus II "Network-in-a-Box" is scheduled to be 
added to the MET in August 1988. 

The software includes the operating system (iRMK), intercomputer 
communication software from National Instruments, external interface, 
executives, operator support, data recording, and math model modules. 

The math model application software includes (1) a multi-link manipulator 
rigid body dynamics model and (2) joint servo models for each of the 
driven manipulator joints. The MET software is intended to provide a 
simulation of proposed manipulator designs. The math model software is 
distributed across the PC's such that a hardware component may be 
substituted and consequently the math model removed. 

The objective of this study is to evaluate the math models developed 
for the MET simulation of a manipulator's rigid body dynamics and the 
servo systems for each of the driven manipulator joints. Specifically, 
the math models were examined with regard to their amenability to pipeline 
and parallel processing. Based on this evaluation and the project 
objectives, a set of prioritized recommendations are offered for future 
work. 


MODEL EVALUATION 


A recursive rigid body dynamics formulation was developed for 
real-time simulation on the MET by Nasser [4] . Nasser considers it to be 
a generalization of a method due to Armstrong [5] . The formulation is 
general in the sense that translational or rotational joint types can be 
modeled. Nasser also claims that the model assumes a topological tree 
configuration, but doesn't describe how this would be done. Procedures 
described in [6] - [7] are particularly suitable to extending Nasser's 
method to general topological trees. 
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The procedure for solving for the reaction force and torques at each 
joint and the joint accelerations given the external forces and moments on 
each link and the actuator forces and moments can be described by 
referring to Figure 1 [4]. 



Figure 1. Open Kinematic Chain 


First, set up the equations of dynamic equilibrium for link N and then 
solve for the reaction forces and moments and the joint acceleration in 
terms of the remaining variables . The remaining variables include the 
proximal links' joint displacements velocities and accelerations and the 
distal links' applied forces and moments, displacements, and velocities. 
Next, proceed to link N-l and substitute the expressions for the reactions 
that involve link N-l joint accelerations into the equations of dynamic 
equilibrium for link N-l. The reaction forces and moments exerted by link 
N-2 on link N-l and the joint acceleration of link N-l should then be 
solved for in terms of the remaining variables. This step is repeated 
until link 1 is reached. The joint acceleration for link 1 is then 
obtained in terms of the distal links' external forces and moments, 
actuator forces and moments and state variables, the state of link 1, and 
the base state. Then move distally to link 2 and substitute for the link 
1 joint acceleration to find the joint 2 acceleration. This step is 
repeated until link N is reached and all the joint accelerations are 
found. 

A free body diagram for any link i of the open kinematic chain is 
shown in Figure 2. 


14-4 




Figure 2 . Link Free Body Diagram 


F* is the link i inertia force, 

-i 

N* is the link i inertia torque, 

ffj is the control force exerted on link i at joint i-1, 
n^ is the control torque exerted on link i at joint i-1. 


f c 

— 1+1 


is the control force exerted on link i+1 at joint i. 


n c is the control torque exerted on link i+1 at joint i, 
-i+1 

f is the reaction force that link i+1 exerts on link 

— i+1 , i 


n 

— i+1 , i 
f 

“i-l,i 


is the reaction force that link i+1 exerts on link i, 
is the reaction moment that link i-1 exerts on link i 
is the reaction force that link i-1 exerts on link i, 
is the reaction torque that link i-1 exerts on link i. 


n is the reaction torque that link l-i exerts on link l, 

— i-1 ,i 

is the sum of the external forces acting on link i referred to 
the mass center of link i, 


N' is the sum of the external moments acting on link i when the 
i external forces are referred to the mass center of link i, 


r* is the position vector of the mass center of link i relative to 
the origin of the i-1 frame, 
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P c § is the position vector of the mass center of link i relative to 

* the origin of the i frame, 

P is the position vector of the origin of frame i relative to the 

* origin of frame i-1 

A number of researchers have investigated pipeline and parallel 
implementations of the equations of motion for various purposes. Luh and 
Lin [8] described a procedure for scheduling the sub tasks of a group of 
microprocessors computing the inverse dynamics using the Newton-Euler 
equations of motion. One microprocessor is assigned to each manipulator 
link. A variable branch-and-bound search finds an optional subtask- 
ordered schedule for each of the microprocessors. However, the total 
processing time of solving the minimum-time scheduling problem could not 
be easily reduced to a manageable level. 

Lathrop [9] proposed two parallel algorithms for solving the inverse 
dynamics problem using a group of special-purpose processors. One is a 
linear Newton-Euler algorithm, which is most closely related to the method 
proposed by Luh and Lin. The other is a logarithmic parallel Newton-Euler 
algorithm based on the "partial sum" technique, which achieves on 0(flog 2 Nl) 
total processing time. However, both methods have two main effects that 
deteriorate the performance of parallel computations. They both require 
potentially massive internal buffering to achieve pipelined computation 
between forward and backward recursions. They both involve complex inter- 
communication and bussing, which frequently cause data to be fetched, and 
as a result data for operand pairs are not properly aligned for parallel 
computations. 

Lee and Chang [10] proposed an algorithm for solving the inverse 
dynamics problem of an N-link manipulator using p processors in parallel 
on the Newton-Euler equations of motion. It was the "recursive doubling" 
technique to achieve a total processing time of OCk^fN/p] + k 2 flog 2 pl), 
where kj, and k 2 are constants. When p - N, the algorithm is 0(flog 2 N] ), 
the same as Lathrop. The p-fold parallel algorithm consists of p-parallel 
blocks with pipelined elements within each parallel block. The results 
from the p-parallel blocks form a homogeneous linear recurrence of size p. 
The parallel algorithm can also be implemented in a systolic pipelined 
architecture, requiring three floating-point operations for each complete 
set of joint torques. 

Binder and Herzog [11] described an algorithm based on the recursive 
Newton-Euler equations of motion where the parallel computations are 
distributed over multiple computing elements, one for each joint. 
Concurrency is achieved by substituting "predicted" values for the actual 
values of variables involved in the recursive equations. The authors 
simulated this method and compared it to other approaches such as 
simplification of the dynamic equations. 

Amin-Javaheri and Orin [12] proposed systolic architectures 
consisting of 1, N, and N(Nfl)/2 processors for computing the inertia 
matrix of an N degree of freedom manipulator. A VLSI-based Robotics 
Processor is being developed as the fundamental component of the 
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architecture. The algorithm used is based on recursive computation of the 
inertial parameters of sets of composite rigid bodies and is programmed to 
exploit any inherent parallelism. 

Lee and Chang [13] proposed two parallel algorithms for computing the 
forward dynamics for real-time simulation with N processors for an N 
degree of freedom manipulator. The first algorithm is based on Walker and 
Orin's [14] composite rigid-body method. It generates the inertia matrix 
using the parallel Newton-Euler algorithm, the parallel linear recurrence 
algorithm [10], and the modified row-sweep algorithm, and then inverts the 
inertia matrix to obtain the joint acceleration vector. The time complex- 
ity of this parallel algorithm is of the order 0(N 2 ) with 0(N) processors. 
Further reduction of the order of time complexity can be achieved by 
implementing the Cholesky factorization procedure on VLSI array processors 
to invert the symmetric, positive-definite, inertia matrix. The second 
parallel algorithm is based on Walker and Orin's [14] conjugate gradient 
method . 

Of the methods just described, only Nasser's method [4] was 
formulated to allow for hardware substitution in a simulation. Currently, 
it is implemented on the MET in a sequential algorithm with one IBM PC-AT. 
The Nasser method appears to be a variation of a method by Featherstone 
[15]. The A? matrix in Nasser's method seems to correspond to the 1^ 
matrix in Featherstone *s method. This would imply that Nasser's method 
has the same nonlinear recurrence that makes Featherstone 's method 
impossible to solve with known methods of parallel solution of recurrence 
problems [13]. Kung [16] has shown that any parallel algorithm using any 
number of processors cannot be essentially faster than the obvious 
sequential algorithm for any nonlinear polynomial recurrence problem. If 
not for this, the recursive doubling technique of Kogge and Stone [17] 
could have been used. The method of Nasser does not seem to be particu- 
larly adaptable to pipelining either. The calculation of the A* matrix 
alone requires about 42% of the total computations. The other processors 
in the pipeline would have to be idle, waiting for the A* to complete. 

The hardware that makes up the MET limits the use of pipelining. 
Instruction pipelining is possible on the Intel 80286 only by prefetching 
instructions. As an alternative, the MET would seem to be an ideal 
application for the new Intel 80960 architecture which allows parallel and 
out-of-order instruction execution of 32-bit operations at 20 MHz [18] . 

The Multibus II "Network-in-a-Box" promises to greatly speed up 
interprocessor communications with initially three Intel 80386 processors 
on board. The National Instruments IEEE-488 interface card and software 
are much too slow. 


RECOMMENDATION AND PRIORITIZING OF FUTURE WORK 

My recommendation of future work for the MET will be listed in order 
of decreasing priority. 
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1. The Nasser algorithm does not appear to be amenable to pipeline 
or parallel processing. One option would be to run it as a very fast 
sequential algorithm on a fast processor. Another option would be to 
modify the method of Binder and Herzog [11] and use prediction to estimate 
values, instead of waiting for a complete set of computed values to be 
ready. 

2. The simulation should include the rigid body dynamics of the base 
link. A mobile base, where motion relative to either a fixed base or base 
with rigid body dynamics should be included in the model. 

3. The servo systems are decoupled in their current design. This 
makes these computations naturally parallel. The rigid body model of 
Nasser could run on a fast processor and broadcast results to N 
processors, one for each of the N servo systems, operating in parallel. 

4. Flexible links should be added to the model. For open 
topological trees, this should be straight forward since the complete 
force system acting on each link is known. A beam element model of a 
flexible manipulator is described in Kelly and Huston [19]. 

5. Constrained motion should be modeled so that closed loops can be 
simulated. Rigid links should be simulated first and then flexible links. 
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